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Abstract 


This  contract  provided  support  for  the  procurement  of  equipment  to  build  and  operate  improved  near-field 
scanning  optical  microscopes  (NSOMs).  Improved  scanning  stages  were  constructed,  and  a  scheme  for 
high  speed  scanning  developed  based  on  distance  regulation  with  tuning  forks.  In  this  scheme  the  phase  of 
the  drive  to  the  fork  is  monitored  with  a  digital  phase-locked  loop,  which  allows  very  rapid  adjustment  of 
the  fork  position  in  response  to  damping  changes  that  occur  as  the  fork  and  optical  fiber  tip  are  scanned  in 
close  proximity  to  a  surface.  In  addition,  new  excitation  and  detection  schemes  were  developed,  which 
have  allowed  us  to  pioneer  the  use  of  nonlinear  NSOM  for  the  study  of  ferroelectric  and  ferromagnetic 
domains,  and  the  crystal  structure  of  surfaces 
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SUMMARY 


Our  research  on  near-field  scanning  optical  microscopy,  and  primarily  its  application  in  nonlinear  studies  of 
surfaces  was  greatly  enhanced  by  the  additional  scanning  and  detection  equipment  procured  from  contract 
funds.  We  have  used  nonlinear  near-field  imaging  to  study  various  interesting  characteristics  of 
technologically  important  anisotropic  materials  such  as  piezoceramics,  ferroelectrics,  and  ferromagnetics. 
These  studies  have  included  studies  of  both  thin  film  and  bulk  devices.  Near-field  imaging  with  bare 
tapered  fiber  optic  probes,  which  is  a  technique  pioneered  by  us,  provides  minimally  invasive  information 
about  domain  alignment,  poling  directions,  domain  structure,  and  domain  dynamics  in  both  passive  and 
active  devices  with  a  resolution  on  the  order  of  25-50nm.  Ferroelectric  materials  such  as  BaTi03  and 
Pb(ZrxTij.x)03  (PZT)  are  of  special  interest  because  these  materials  are  being  increasingly  used  in  memory 
systems.  We  have  shown  that  the  location  and  alignment  of  domains  in  these  materials  can  be  determined 
by  examining  the  polarization  dependence  of  second  harmonic  generation  (SHG)  in  the  near-field.  We  have 
improved  the  signal-to-noise  ratio  in  these  experiments  by  using  femtosecond  Ti:sapphire  laser  excitation, 
since  the  phenomena  we  are  studying  provide  only  weak  signals  unless  femtosecond  pulsed  lasers  of  high 
peak  intensity  are  used. 

We  have  also  carried  out  technological  improvements  on  scanning  techniques  in  NSOM  using  tuning  fork 
distance  control  and  rapid  digital  phase-locked  loop  techniques  to  acquire  tip  to  sample  separation  and 
permit  fast  scanning  and  a  shorter  residence  time  per  pixel 

INTRODUCTION 

Near-field  scanning  optical  microscopy  (NSOM)  has  become  a  relatively  mature  technique  during  the  last 
ten  years,  although  there  continue  to  be  difficulties  in  image  interpretation  when  this  technique  is  used  as  a 
pure  imaging  tool.  However,  the  use  of  near-field  probes  for  specific  excitation  or  interrogation  of  surfaces 
suffers  from  fewer  ambiguities.  This  is  especially  true  when  bare  tapered  fiber  probes  with  end  diameters  of 
50nm  or  less  are  used.  We  have  used  such  probes  for  direct  observation  of  near-field  emission  of 
developmental  diode  lasers,  for  studying  the  dynamics  of  surface  plasmons,  for  the  study  of  lithography  on 
a  lOOnm  scale,  for  studying  the  effects  of  surface  roughness  on  SHG  at  metal  films,  for  visualization  of 
ferromagnetic  domains  in  ferrromagnetic  materials,  and  for  examining  photon  bandgap  effects  in  periodic 
optical  structures.  Near-field  probes  with  a  metal  coating,  which  are  very  commonly  used  for  high 
resolution  suffer  from  artefacts,  especially  the  excitation  of  surface  plasmons  in  the  tip  coating  itself,  which 
makes  them  unsuitable  for  minimally  invasive  probing.  Bare  probes  provide  a  local  field  enhancement, 
although  this  is  not  as  great  as  can  be  achieved  with  sharp  metal  probes.  However,  this  can  still  give  a 
substantial  enhancement  in  local  SHG,  which  depends  on  the  fourth  power  of  the  applied  field. 

NANOMETER  SCALE  CRYSTAL  ANALYSIS  USING  NEAR-FIELD  OPTICAL  SECOND 
HARMONIC  MICROSCOPY:  APPLICATION  TO  THIN  FERROELECTRIC  FILMS 

In  this  section  we  provide  details  of  how  we  carry  out  nonlinear  NSOM  in  a  particular  case  to  illustrate  the 
scope  of  the  research  we  have  carried  out.  We  describe  how  near-field  optical  scanning  for  second 
harmonic  imaging  using  tapered  optical  fiber  tips  externally  illuminated  with  femtosecond  laser  pulses  is 
carried  out.  Enhancement  of  the  electric  field  at  the  tip  of  the  fiber  results  in  enhanced  second  harmonic 
generation  (SHG)  from  the  sample  region  near  the  tip.  This  SH  emission  is  collected  by  the  same  tapered 
fiber.  In  preliminary  experiments  polarization  properties  of  second  harmonic  generation  from  a  poled  single 
crystal  of  BaTi03  have  been  studied  experimentally  and  theoretically  in  the  near  and  far  field  zones  for 
different  orientations  of  poling  direction  with  respect  to  the  tip  and  the  incoming  fundamental  light.  Good 
agreement  between  theory  and  experiment  has  been  demonstrated.  Symmetry  properties  of  the  near-field 
SH  signal  allow  us  to  recover  the  local  poling  direction  of  individual  ferroelectric  domains  in 
polycrystalline  thin  ferroelectric  films  with  spatial  resolution  of  the  order  of  80  nm.  Near-field  strain- 
induced  SHG  from  thin  BaxSr1.xTi03  films  has  also  been  studied,  and  excellent  agreement  between  theory 
and  experiment  has  been  established.  Thus,  nonlinear  near-field  scanning  microscopy  has  been  shown  to 
be  a  novel  tool  for  nanometer  scale  crystal  analysis  of  polycrystalline  samples. 


NONLINEAR  GENERATION  AT  SURFACES. 


Optical  second  harmonic  generation  (SHG)  is  an  extremely  sensitive  technique  for  characterization  and 
investigation  of  a  wide  variety  of  materials,  surfaces  and  adsorbates  [1].  SHG  is  especially  useful  for 
recovery  of  the  symmetry  properties  of  the  sample  under  investigation.  Far-fleld  SH  microscopy  has 
already  demonstrated  its  great  potential  in  such  fields  as  magnetic  imaging  [2]  and  metal-oxide- 
semiconductor  device  characterization  [3].  However,  in  many  applications  a  spatial  resolution  below  the 
diffraction  limit  of  far-field  optics  is  required.  For  example,  better  than  100  nm  resolution  is  required  in 
characterization  of  thin  ferroelectric  films  such  as  BaTi03  and  Pb(ZrxTij_x)03  (or  PZT),  which  form  the 
basis  of  a  new  thin  film  technology  for  data  storage  [4].  The  ultimate  performance  (density  and  speed  of 
operation)  of  these  devices  depends  on  the  spatio-temporal  properties  of  individual  ferroelectric  domains. 
A  better  understanding  of  these  properties  requires  experiments  performed  with  nanometer  spatial  and  sub¬ 
nanosecond  temporal  resolution.  Near-field  optical  microscopy  may  be  the  technique  of  choice  for  such 
experiments.  Its  main  advantages  in  comparison  with  other  scanning  probe  techniques  are  the  possibility  of 
fast  time-resolved  measurements,  and  substantially  smaller  perturbation  of  the  sample  under  investigation 
caused  by  the  optical  probe.  As  a  result,  many  attempts  have  been  made  recently  [5-7]  to  implement 
different  near-field  optical  techniques  for  SH  imaging.  Unfortunately,  even  the  best  resolution  reported  so 
far,  140  nm  [7]  obtained  with  a  widely  adopted  aperture  near- field  probe  [8],  falls  short  of  the  desired 
range.  The  main  experimental  problem  in  near-field  SH  imaging  is  low  optical  signal,  since  SHG  is 
generally  a  weak  effect,  and  it  is  extremely  difficult  to  detect  using  an  aperture  fiber  probe  with  a  typical 
throughput  of  less  than  10'4  for  a  100  nm  aperture. 

An  apertureless  approach  to  near-field  microscopy  [9]  is  usually  based  on  local  field  perturbation  or  on 
the  local  field  enhancement  by  laser-illuminated  sharp  metal  tips.  In  the  latter  case  the  metal  tip  provides  a 
local  excitation  source  for  the  spectroscopic  response  of  the  sample  under  investigation.  Unfortunately,  this 
approach  is  difficult  to  apply  to  SH  microscopy,  since  a  laser  illuminated  metal  tip  would  be  a  source  of 
SHG  comparable  to  the  near-field  SH  signal  from  the  sample.  In  this  paper,  we  present  a  near-field  optical 
technique  for  second  harmonic  imaging  using  a  bare  tapered  optical  fiber  tip  externally  illuminated  with 
femtosecond  laser  pulses.  We  show  theoretically  and  experimentally  that  under  illumination  with  P- 
polarized  light  the  local  SHG  from  the  sample  region  under  the  tip  is  strongly  enhanced  even  in  the  case  of 
a  dielectric  tip.  At  the  same  time,  the  glass  tip  itself  does  not  generate  any  measurable  SH  light.  The  local 
SH  emission  from  the  sample  region  under  the  tip  is  collected  by  the  same  tapered  fiber,  which  has  a  field 
of  view  on  the  order  of  half  the  wavelength  of  SH  light  X2J2  in  diameter,  where  ^2co=400  nm  (this  is  the 
resolution  limit  of  a  linear  near-field  microscope  using  bare  tapered  fibers).  This  configuration 
automatically  solves  one  of  the  main  problems  in  apertureless  near-field  microscopy:  rejection  of  the 
background  signal  from  the  entire  illuminated  sample  area.  As  a  result,  spatial  resolution  on  the  order  of  80 
nm  has  been  achieved  in  imaging  of  thin  ferroelectric  films.  Symmetry  properties  of  the  near-field  SH 
signal  studied  experimentally  and  theoretically  have  allowed  us  to  recover  local  poling  direction  of 
individual  ferroelectric  domains  in  polycrystalline  thin  ferroelectric  films.  Thus,  the  developed  technique 
may  be  considered  as  a  novel  tool  for  nanometer  scale  crystal  analysis  of  polycrystalline  samples. 
Measurements  of  strain  induced  near-field  SHG  in  thin  BaxSr!_xTi03  films  have  also  allowed  us  to  establish 
near-field  second  harmonic  microscopy  as  a  unique  tool  for  nanometer  scale  strain  analysis  of  thin  films. 

EXPERIMENTAL  SETUP 

One  experimental  setup  we  have  used  is  shown  in  Fig.l.  Weakly  focused  light  (illuminated  spot  diameter 
on  the  order  of  50  pm)  from  a  Ti:sapphire  laser  system  consisting  of  an  oscillator  and  a  regenerative 
amplifier  operating  at  810  nm  (repetition  rate  up  to  250  kHz,  100-fs  pulse  duration,  and  up  to  10  pJ  pulse 
energy)  was  directed  at  an  angle  of  incidence  around  50°  onto  the  sample  surface.  Excitation  power  at  the 
sample  surface  was  kept  well  below  the  ablation  threshold.  The  local  SHG  was  collected  using  an  uncoated 
adiabatically  tapered  fiber  tip,  which  is  drawn  at  the  end  of  a  single  mode  fiber  by  the  standard  heating  and 
pulling  procedure.  Tip  curvature  was  measured  using  a  scanning  electron  microscope.  A  fiber  tip  with  an 
end  diameter  of  about  50  nm  was  scanned  over  the  sample  surface  with  a  constant  tip-surface  distance  of  a 
few  nanometers  using  shear  force  feedback  control  [8].  Therefore,  surface  topography  images  were 
obtained  while  simultaneously  recording  the  SH  near-field  images.  A  typical  set  of  such  images  of  a  0.5 
pm  thick  pulsed  laser  deposited  PZT  film  grown  on  a  Si  wafer  substrate  obtained  under  excitation  with  P- 
and  S-  polarized  fundamental  light  is  presented  in  Fig.2.  The  topographical  image  in  Fig.2(a)  shows  the 
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grain  structure  of  the  film.  Higher  than  usual  noise  in  the  image  is  caused  by  less  stable  behavior  of  the 
shear-force  feedback  control  when  the  tip  is  illuminated  with  femtosecond  laser  pulses.  This  may  be  caused 
by  local  evaporation  of  the  thin  water  layer  absorbed  on  the  sample  surface,  which  is  believed  to  contribute 
to  the  shear  force.  The  SH  image  in  Fig.2(b)  obtained  under  illumination  with  P-polarized  fundamental 
light  clearly  shows  fine  structure  of  SHG  within  individual  grains.  Fundamental  light  distribution  was  also 
measured  and  did  not  show  any  variations  on  this  scale. 

RESULTS  AND  DISCUSSIONS 

1.  Study  of  near-field  SHG  from  single  crystals  of  BaTi03. 

Optical  SHG  is  widely  used  to  characterize  the  average  structure  and  crystal  orientation  of  ferroelectric 
films  [10].  Near-field  second  harmonic  microscopy  is  the  logical  extension  of  these  far-field  techniques.  In 
order  to  understand  its  main  contrast  mechanisms  (and,  as  a  result,  understand  the  origin  of  the  fine 
structure  in  Fig.2)  we  performed  a  study  of  the  polarization  properties  of  the  near-  and  far-field  SH  signal 
using  a  poled  single  crystal  of  BaTi03  as  a  sample  (BaTi03  and  Pb(ZrxTii_x)03  have  similar  perovskite 
structure  and  belong  to  the  same  tetragonal  symmetry  class).  Phase-matched  SHG  is  prohibited  in  the  bulk 
of  BaTi03  because  of  its  strong  dispersion  [11].  As  a  result,  the  measured  SH  signal  originates  at  the 
surface,  which  makes  the  experimental  situation  look  very  similar  to  the  case  of  a  thin  film.  Near-field  SH 
signal  dependencies  on  the  polarization  of  the  fundamental  light  for  three  different  poling  (optical  axis) 
directions  of  BaTi03  crystal  are  shown  in  Fig.3.  The  SH  signal  exhibits  approximate  4-fold  symmetry  with 
respect  to  polarization  angle  rotation  in  the  case  (a)  of  poling  (optical  axis)  direction  located  in  the 
incidence  plane  of  the  fundamental  light  parallel  to  the  crystal  surface.  The  symmetry  becomes  2-fold  in  the 
case  (b)  when  the  optical  axis  is  perpendicular  to  the  incidence  plane.  In  case  (c)  (when  the  poling 
direction  is  perpendicular  to  the  crystal  surface)  SH  signal  is  much  weaker  than  in  cases  (a)  and  (b).  Big 
differences  between  these  dependencies  indicate  that  the  near-field  SHG  may  be  used  to  recover  local 
poling  direction  in  thin  ferroelectric  films.  In  order  to  do  this,  we  must  achieve  some  theoretical 
understanding  of  the  reasons  for  such  a  difference.  Another  important  question  is:  how  good  is  the  spatial 
resolution  of  the  near-field  SH  microscopy  in  the  experimental  geometry  discussed  above? 

Surprisingly,  most  of  the  features  of  the  near-field  SHG  may  be  understood  from  a  very  simple 
theoretical  model.  Let  us  consider  fundamental  light  illuminating  the  tip  of  the  microscope  located  near  the 
BaTi03  crystal  surface.  Three  different  cases  of  the  poling  vector  orientation  are  shown  in  Fig.4(a,b,c), 
where  a  is  the  incidence  angle,  and  the  z-axis  is  chosen  to  coincide  with  the  poling  direction  of  the  crystal. 
All  three  cases  may  be  considered  in  a  similar  way,  so  we  will  concentrate  on  case  (b)  and  only  present  the 
final  results  for  the  cases  (a)  and  (c).  Let  us  consider  the  fundamental  optical  field  distribution  in  the  tip- 
sample  region  at  distances  much  smaller  then  the  wavelength  X  of  the  fundamental  light.  In  this  region  the 
quasi-electrostatic  approximation  may  be  used.  The  tip  shape  may  be  approximated  as  an  ellipsoid  with  a 
very  high  aspect  ratio.  For  P-polarized  excitation  light  and  y-component  (the  component  perpendicular  to 
the  sample  surface)  of  the  optical  field  the  exact  analytic  solution  of  the  Laplace  equation  [12]  gives 
Etip(co)=eEyo(co)  where  s  is  the  dielectric  constant  of  the  tip,  Ey0  is  the  incoming  field  and  Etip  is  the  y- 
component  of  the  field  just  below  the  tip  apex.  As  a  result,  for  the  incidence  angles  a  around  n/2  the  local 
SHG  under  the  tip  is  enhanced  by  a  factor  of  n8  (since  I(2co)tip=  E2,ip(2co)=E  4tjp(«)=n8E4yo(co)),  where  n  is  the 
tip  refractive  index.  Assuming  n=l  .6  the  enhancement  factor  equals  about  40.  It  may  be  even  higher  (of  the 
order  of  (2s)4  )  in  the  vicinity  of  a  sample  surface  with  a  high  dielectric  constant,  such  as  a  metal  or  a 
ferroelectric,  due  to  the  image  potential.  Exact  tip  enhancement  factors  in  non-linear  near-field  microscopy 
were  calculated  by  Y.  Kawata  et  al  [13]  for  different  tip-sample  pairs.  For  example,  in  the  case  of  a 
diamond/air/diamond  tip-sample  gap  E4tipmay  be  enhanced  by  a  factor  of  600.  However,  in  our  discussion 
to  follow,  we  will  not  consider  any  particular  value  of  the  field  enhancement  factor.  Instead,  we  denote  the 
fundamental  field  enhancement  factor  by  y.  The  fundamental  optical  field  components  under  the  tip  apex  at 
a  polarization  angle  (j)  can  be  written  approximately  as: 

Ex  =  Ecos<j)cosa,  Ey  =  (y/sb)Ecos(|)sina,  Ez  =  Esin<|>,  ( 1 ) 

where  eb  is  the  dielectric  constant  of  BaTi03,  (J)=0  and  <J)=90°  correspond  to  P-polarized  and  S-polarized 
excitation  light  respectively.  Taking  into  account  the  non-zero  components  of  the  second  harmonic 
susceptibility  tensor  for  BaTi03:  x(2\yz  =  x(2)xxz  =  X(2)yzy  =  X(2)xzx  =17.7xl0'12  m/V,  x(2)zyy  =  X(2)zxx  =  - 
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18.8xlO'12  m/V,  and  x(2)zzz  =  -7.1xl0"12  m/V  (where  z  is  the  poling  direction)  [14],  the  second  harmonic 
field  components  at  the  tip  apex  may  be  written  as: 

D(2)x  =  x(2)XzxE2sin2(j)COsa,  D(2)y  =  (y/8b)x(2)xzxE2sin2(|)sina,  (2) 

D(2)z  =  (y2/ Eb2)x(2)zxxE2cos2(()sin2a  +  x{2)zXxE2cos2(()COS2a  +  x{2)zzzE2sin2(t) 

If  we  assume  that  the  microscope  tip  collects  only  dipole  radiation,  the  second  harmonic  optical  signal  will 
be  proportional  to 

I(2)  =  d^z2  +  D(2)x2  (3) 

(SH  dipoles  that  oscillate  along  the  y-direction  do  not  radiate  towards  the  tip).  Similar  calculations  of  the 
field  at  the  tip  apex  are  easy  to  perform  for  cases  (a)  and  (c).  Let  us  point  out  that  within  our  model  the 
angular  dependencies  of  the  field  E  for  the  points  in  space  around  the  tip  apex  should  be  approximately  the 
same  as  in  (1).  Only  the  enhancement  factor  y  is  different  for  these  points  (its  value  changes  from  a 
maximum  just  under  the  tip  apex  to  y=l  far  from  the  apex).  So,  in  the  final  result  one  must  replace  the 
factors  y4  and  y2  by  the  average  values  <y4>  and  <y2>  taken  over  all  the  area  around  the  tip  apex  where  the 
SH  signal  is  collected.  The  diameter  of  this  area  is  approximately  X2J3  -  X2J2  which  is  the  theoretical 
resolution  limit  of  near-field  optical  microscopy  in  the  collection  mode  performed  with  uncoated  fiber  tips 
[15]. 

The  final  results  of  our  calculations  are  presented  in  Fig. 5.  The  best  agreement  with  experimentally 
measured  symmetry  properties  of  the  near-field  SH  polarization  curves  was  achieved  in  the  y  =  5  -7  range 
for  the  field  enhancement  factor.  This  is  reasonably  close  to  the  expected  value  of  y  =  28^.  A  much  smaller 
value  of  the  SH  signal  in  case  (c)  is  a  consequence  of  the  fact  that  I(2)  contains  terms  proportional  to  <y4>  in 
the  cases  (a)  and  (b),  but  in  case  (c)  the  SH  signal  I(2)  is  proportional  to  <y2>  ; 

I(2)  =  D(2)x2  +  D(2)y2  =  (<y2>/8b2)x(2)xzx2E4(sin22(j)sin2a+cos4(t)sin22a)  (4) 

The  transition  from  near-field  mode  to  far-field  mode  operation  of  the  microscope  (which  corresponds  to 
an  increase  in  the  tip-sample  separation  from  a  few  nanometers  to  a  distance  of  a  few  wavelengths  of  light) 
can  be  described  mathematically  within  our  model  as  a  transition  from  y  =  5  -  7  to  y=l  (no  field 
enhancement  by  the  tip).  Results  of  numerical  calculations  for  such  a  transition  performed  using  (3)  are 
shown  in  Fig.6.  Experimental  data  shown  in  the  same  figure  are  in  good  qualitative  agreement  with  the 
calculations:  both  theoretical  and  experimental  data  have  the  same  symmetry  of  the  SH  polarization 
dependencies  in  the  near-field  and  far-field  zones.  This  indicates  that  our  simple  model  provides  an 
adequate  description  of  the  essential  physics  involved  in  near-field  observation  of  SHG  from  ferroelectric 
samples. 

2.  Measurements  of  local  poling  direction  in  thin  BaTi03  and  PZT  films. 

We  have  performed  similar  local  polarization  properties  measurements  in  different  areas  of  thin  BaTi03 
and  Pb(ZrxTi]_x)03  films.  An  example  of  such  measurements  is  presented  in  Fig. 7.  The  data  obtained  at 
points  A,  B  and  C  of  Fig. 7(b)  show  close  resemblance  to  the  curves  (a),  (b)  and  (c)  respectively  in  Fig. 3. 
Optical  properties  of  Pb(ZrxTii_x)03  and  BaTi03  are  similar  to  each  other,  with  Pb(ZrxTit_x)03  (or  PZT) 
having  somewhat  higher  values  of  refractive  index  and  nonlinear  susceptibility:  no=2.610  and  ne=2.611  at 
800  nm;  n0=2.966  and  ne=2.967  at  400  nm;  x(2)yyZ  =  X(2)xxz  =  X(2)yzy=  X<2)xzx  =37.9x1  O'12  m/V,  x<2)zyy=  X(2) zxx  = 
-  42.8xl0'12  m/V,  and  y}2)za=  -8.5xl0'12  m/V  at  1060  nm  [16,17],  Thus,  our  theoretical  analysis  performed 
for  BaTi03  should  also  be  applicable  to  PZT  films.  As  a  result,  the  localized  areas  of  dark  and  bright  SHG 
in  Fig. 7(b)  should  be  identified  as  regions  with  different  poling  (optical  axis)  directions.  Similarly,  the 
areas  of  bright  SHG  from  Fig.2(b)  localized  within  bigger  crystalline  grains  (visible  in  topography),  and 
separated  by  thin  dark  boundaries,  should  be  identified  as  individual  ferroelectric  domains  separated  by 
domain  boundaries.  A  similar  multidomain  structure  of  grains  in  thin  pulse  laser  deposited  PZT  films  was 
observed  in  [18].  Very  high  spatial  resolution  of  the  near-field  SH  image  is  evident  from  Fig.2(d)  where  the 
apparent  thickness  of  a  domain  boundary  is  equal  to  80  nm.  Such  a  high  spatial  resolution  may  seem 
surprising  since  theoretical  estimates  of  the  resolution  obtained  with  uncoated  fiber  tips  in  linear  near-field 
measurements  varies  within  the  X/2  -  X!3  range  [15].  The  following  simple  estimate  shows  that  in  the  SH 
mode  the  resolution  obtained  with  adiabatically  tapered  uncoated  tips  can  be  an  order  of  magnitude  better. 
From  the  theoretical  discussion  above  and  the  experimental  data  shown  in  Fig.6  it  is  clear  that  the 
fundamental  field  component  parallel  to  the  microscope  tip  (and  perpendicular  to  the  sample  surface)  is 
enhanced  by  a  factor  y  of  the  order  of  y  =  5  -  7.  This  leads  to  the  local  SH  signal  enhancement  at  the  tip 
apex  by  a  factor  of  up  to  y4  for  P-polarized  excitation  light.  The  fiber  tip  of  the  microscope  collects  SHG 
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from  a  sample  area  limited  by  a  circle  with  a  diameter  on  the  order  ofX2J2,  where  \2(0  is  the  wavelength  of 
SH  light.  But  the  sample  area  just  below  the  tip  dominates  the  resulting  SH  signal.  As  a  result,  the 
resolution  of  the  microscope  in  the  SH  mode  is  determined  by  the  fiber  tip  diameter  d,  which  was  on  the 
order  of  50  nm  in  our  experiments.  Even  better  resolution  may  be  expected  using  thinner  tips  or  tips  with 
higher  refractive  indices.  The  tip  contribution  and  the  contribution  from  all  the  other  sample  surface 
become  equal  at  a  tip  diameter  determined  by  y4d2=(X2o/2)2  or  d=X2J(2y2)  which  may  be  as  small  as  10-30 
nm.  In  fact,  more  severe  limitations  on  the  resolution  of  the  microscope  are  caused  by  the  small  number  of 
SH  photons  collected  by  the  microscope  tip.  In  our  setup  the  SH  signal  was  measured  with  a 
photomultiplier  and  a  gated  photon  counter.  The  gate  was  open  only  for  the  duration  of  a  single  photon 
pulse  from  the  photomultiplier  in  sync  with  the  femtosecond  pulses  of  the  laser.  The  SH  signal  at  every 
point  of  the  image  was  averaged  over  100-300  ms.  The  characteristic  SH  photon  counting  rate  was  on  the 
order  of  one  photon  count  per  100-300  laser  pulses.  Thus,  the  image  acquisition  time  was  on  the  order  of 
20-30  minutes.  Further  increase  of  the  microscope  resolution  is  limited  by  the  time  of  measurement  rather 
than  by  achievement  of  a  better  tip  geometry. 

SH  imaging  performed  with  S-polarized  excitation  light  provides  another  verification  of  the  proposed 
mechanism  for  high  spatial  resolution.  Since  fundamental  field  components  perpendicular  to  the  tip  axes  do 
not  experience  any  substantial  enhancement,  one  should  expect  much  smaller  SH  signal  collected  in  this 
case.  SH  images  in  Fig.2(b)  and  Fig.2(c)  obtained  from  the  same  region  of  the  PZT  film  under  P-polarized 
and  S-polarized  excitation  light  clearly  show  the  expected  behavior. 

Another  proof  of  high  spatial  resolution  achieved  in  our  measurements  of  local  poling  direction  has  been 
obtained  by  observation  of  sudden  symmetry  changes  in  SH  polarization  dependencies  performed  with  45 
nm  spatial  steps  (Fig.8).  After  obtaining  the  topographical  and  SH  images  shown  in  Fig.8(a)  and  Fig.8(b) 
respectively,  we  have  performed  point  by  point  measurements  of  local  SH  polarization  dependencies  over 
the  region  which  show  clear  change  in  SH  contrast.  These  measurements  were  done  by  moving  the  tip  in  45 
nm  steps  from  left  to  right  through  the  points  on  the  film  surface  indicated  in  both  the  topographical  and  SH 
images.  The  data  obtained  are  shown  in  Fig.8(c),  where  the  curves  are  numbered  in  succession.  A  clear 
transition  from  4-fold  symmetry  (curves  1-3)  to  2-fold  symmetry  (curves  8  and  9)  has  been  observed  with  a 
narrow  intermediate  region  (curves  4-7)  which  does  not  show  any  polarization  dependence  of  SHG. 

3.  Strain  induced  near-field  SHG  in  thin  BaxSr!_xTi03  films. 

We  also  applied  our  near-field  setup  to  visualize  strain  induced  SHG  in  thin  BaxSri_xTi03  (or  BST)  films 
grown  on  LaO  substrates.  At  room  temperature  bulk  BST  has  the  centrosymmetric  cubic  m3m  crystal 
symmetry,  for  which  bulk  SHG  should  be  absent.  The  Curie  temperature  of  BST  for  a  phase  transition  into 
a  ferroelectric  state  is  about  240K.  As  has  been  shown  recently  [19],  thin  BST  films  may  generate  SH  light 
because  of  strain  due  to  the  lattice  mismatch.  The  nonlinear  susceptibility  values  of  these  films  measured  in 
[19]  are  sufficiently  different  from  those  of  BaTi03  and  PZT:  X(2)yyz  =  X(2)xxz  =  X(2)yzy  =  X(2)xzx  =2.2xl0'12 
m/V,  X(2zyy=  X(2)zxx=  “  4.7xl0"12  m/V,  and  x(2)zzz=  -ll.OxlO'12  m/V  (where  z  is  the  strain  direction).  As  a 
result,  theoretical  near-field  SH  signal  polarization  dependencies  calculated  for  different  strain  directions 
and  shown  in  Fig.9  differ  substantially  from  the  dependencies  in  Figs.  3  and  5.  For  instance,  curve  (b) 
obtained  for  a  strain  direction  perpendicular  to  the  incidence  plane  of  the  fundamental  light  shows  SH 
maxima  for  S-polarized  excitation  light.  This  behavior  was  never  observed  with  BaTi03  and  PZT  samples, 
but  it  was  routinely  observed  (see  Fig.  10(c))  in  experiments  conducted  with  BST  films.  Curve  B  in 
Fig.  10(c)  represents  typical  polarization  behavior  of  the  SH  signal  measured  in  different  areas  of  the  SH 
image  in  Fig.  10(b).  After  the  sample  had  been  rotated  by  90°,  the  typical  polarization  behavior  changed 
(curve  A  in  Fig.  10(c)).  Curves  A  and  B  look  similar  to  theoretical  curves  in  Fig.9  (a)  and  (b)  respectively, 
which  is  another  indication  that  our  simple  theoretical  model  provides  an  adequate  description  of  near-field 
SH  measurements. 

Near-field  SH  signal  observed  from  BST  films  was  an  order  of  magnitude  weaker  than  the  signal 
observed  with  PZT  films,  which  required  substantially  longer  image  acquisition  times.  Nevertheless,  SH 
imaging  was  possible  and  an  example  of  such  images  obtained  with  a  0.4  pm  thick  BST  film  are  presented 
in  Fig.l0(a,b).  The  SH  image  (Fig.  10(b))  is  supposed  to  show  the  distribution  of  mechanical  strain  in  the 
film.  Features  as  small  as  200  nm  are  clearly  detectable  in  this  image.  They  are  probably  caused  by  the 
terrace  structure  of  the  LaO  substrate,  which  was  clearly  visible  and  has  the  same  spatial  orientation  as  the 
orientation  of  stripes  in  the  SH  image. 
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In  conclusion,  an  apertureless  near-field  optical  technique  for  second  harmonic  imaging  with  optical 
resolution  on  the  order  of  80  nm  has  been  developed.  Polarization  properties  of  second  harmonic 
generation  from  a  poled  single  crystal  of  BaTi03  have  been  studied  experimentally  and  theoretically  in  the 
near  and  far  field  zones  for  different  orientations  of  poling  direction  with  respect  to  the  tip  and  the 
incoming  fundamental  light.  Good  agreement  between  theory  and  experiment  is  demonstrated.  Symmetry 
properties  of  the  near-field  SH  signal  allow  us  to  recover  local  poling  direction  of  individual  ferroelectric 
domains  in  polycrystalline  thin  ferroelectric  films  with  spatial  resolution  on  the  order  of  80  nm.  Near-field 
strain-induced  SHG  from  thin  BaxSr]_xTi03  films  has  also  been  studied,  and  excellent  agreement  between 
theory  and  experiment  has  been  established.  Thus,  the  described  technique  provides  a  novel  tool  for 
nanometer  scale  crystal  and  strain  analysis  of  polycrystalline  samples. 
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Figure  1.  Schematic  view  of  our  near-field  SH  microscopy  setup.  P-  and  S-  polarization  directions  of 
the  fundamental  light  are  indicated. 
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(C)  (d) 

Figure  2.  Simultaneously  measured  shear-force  topographical  (a)  and  SH  near-field  optical  images  (b,c) 
of  a  0.5  pm  thick  PZT  film  surface  obtained  respectively  with  P-  and  S-polarized  excitation  light.  The  size 
of  the  images  is  3  pm  by  3  pm.  Height  variation  within  the  topographical  image  is  82  nm.  (d)  Cross-section 
of  the  near-field  image  along  the  direction  shown  in  (b).  The  apparent  domain  boundary  width  shown  by 
the  markers  equals  80  nm. 


Figure  3.  Near- field  SH  signal  dependencies  on  the  polarization  of  fundamental  light  for  three  different 
poling  directions  of  BaTi03  crystal  (shown  by  the  arrows)  with  respect  to  the  incoming  fundamental  light 
(shown  by  the  straight  lines).  P-  and  S-  polarization  directions  of  the  fundamental  light  are  indicated. 
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Figure  4.  Orientation  of  the  BaTi03  crystal  with  respect  to  the  illuminating  fundamental  light  for  curves 
(a),  (b),  and  (c)  in  Fig.3. 
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Figure  5.  Theoretical  near-field  SH  signal  dependencies  on  the  polarization  angle  of  fundamental  light 
obtained  in  the  geometry  (a),  (b)  and  (c)  from  Fig.4  for  the  field  enhancement  factor  y  =  7  and  a  =  45°. 

Figure  6.  Theoretical  (open  circles  and  triangles)  and  experimentally  measured  (filled  circles  and 
triangles)  data  for  the  transition  from  the  near-field  to  the  far-field  behavior  of  the  SH  signal  polarization 
dependency  for  geometry  (b)  from  Fig.4. 


90  100  110  120  130  140  150  160  170  180  190  200 

Polarization  angle  (degrees) 


Polarization  Angle  (degrees) 


(c) 


Figure  7.  Simultaneously  measured  (a)  topographical  and  SH  near-field  optical  image  (b)  of  a  PZT 
film.  3  pm  by  3  pm  images  were  obtained  using  P-  polarized  excitation  light.  Height  variation  within  the 
topographical  image  is  90  nm.  Near-field  SH  signal  dependencies  on  the  polarization  of  the  fundamental 
light  taken  at  three  different  points  A,  B  and  C  of  the  images  are  shown  in  (c). 
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(c) 


Figure  8.  Simultaneously  measured  (a)  topographical  and  SH  near-field  optical  image  (b)  of  a  PZT  film. 
3  pm  by  3  pm  images  were  obtained  using  P-  polarized  excitation  light.  Height  variation  within  the 
topographical  image  is  75  nm.  Near-field  SH  signal  dependencies  on  the  polarization  of  the  fundamental 
light  taken  at  9  points  indicated  in  the  images  are  shown  in  (c).  The  spacing  between  the  points  is  45  nm. 
The  points  are  numbered  in  succession  from  left  to  right. 
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Figure  9.  Theoretical  near-field  SH  signal  dependencies  on  the  polarization  angle  of  fundamental  light 
obtained  for  strained  BST  films  in  the  geometry  similar  to  (a),  (b)  and  (c)  in  Fig. 4,  where  the  arrows 
represent  the  direction  of  strain  in  the  film.  The  field  enhancement  factor  is  assumed  to  be  y  =  4  and  the 
incidence  angle  a  -  50°. 


Figure  10.  Simultaneously  measured  (a)  topographical  and  SH  near-field  optical  image  (b)  of  a  BST  film. 
2.2  pm  by  0.9  pm  images  were  obtained  using  S-  polarized  excitation  light.  Height  variation  within  the 
topographical  image  is  50  nm.  (c)  Typical  near-field  SH  signal  dependencies  on  the  polarization  angle  of 
the  fundamental  light  measured  over  a  thin  BST  film.  Curve  B  was  measured  in  the  middle  point  of  the  SH 
image  (b).  Curve  A  was  measured  over  the  same  area  of  the  film  after  the  sample  had  been  rotated  by  90°. 
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Copies  of  some  key  journal  articles  published  during  the  contract  period. 


A  phase-locked  shear-force  microscope  for  distance  regulation  in  near¬ 
field  optical  microscopy 

Walid  A.  Atia  and  Christopher  C.  Davisa> 

Electrical  Engineering  Department ,  University  of  Maryland ,  College  Park,  Maryland  20742 
(Received  16  April  1996;  accepted  for  publication  19  November  1996) 

A  nonoptical  phase-locked  shear-force  microscope  utilizing  a  quartz  crystal  tuning  fork  acting  as  a 
voltage-controlled  oscillator  in  a  phase-locked  loop  has  been  implemented.  A  tapered  optical  fiber 
is  rigidly  mounted  on  one  of  the  prongs  of  the  fork  to  serve  as  both  a  shear-force  pickup  and  a 
near-field  optical  probe.  The  crystal  is  driven  at  its  resonance  frequency  through  positive  feedback 
of  the  monitored  current  through  the  crystal.  This  signal  is  used  as  the  voltage-controlled  oscillator 
in  a  phase-locked  loop.  The  scheme  allows  for  scan  speeds  far  beyond  the  Q-limited  amplitude 
sensor  bandwidth  and  exhibits  excellent  sensitivity  for  a  high-Q  resonator.  Furthermore,  given  the 
small  vibration  amplitude  of  the  tip  (<0.5  nm)  and  the  distance  over  which  it  is  reduced  (>  6  nm), 
it  is  unlikely  that  the  tip  is  making  direct  contact  with  the  sample  surface  as  has  been  suggested  for 
the  optical  shear-force  detection  scheme.  ©  1997  American  Institute  of  Physics. 
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Shear-force  microscopy1,2  has  been  the  method  of  choice 
for  distance  regulation  in  near-field  scanning  optical 
microscopy,3  allowing  one  to  simultaneously  maintain  the 
aperture  within  the  near-field  of  the  sample  and  obtain  a 
topographical  map  of  the  sample  surface.  To  date,  the  most 
popular  implementation  of  shear-force  microscopy  involves 
vibrating  an  optical  fiber  probe  at  its  resonance  frequency  in 
a  direction  lateral  to  the  sample  surface,  and  optically  moni¬ 
toring  the  vibration  amplitude.  This  amplitude  decreases  as 
the  tip  approaches  some  tens  of  nanometers  from  the  sample 
surface.  The  mechanism  for  this  amplitude  reduction  is  still 
not  well  understood.  Recently,  however,  Gregor  et  alA  have 
proposed  that  the  damping  is  simply  due  to  the  tip  edge 
lightly  contacting  the  sample  in  a  manner  similar  to  tapping¬ 
mode  atomic  force  microscopy.5  They  propose  that  as  the  tip 
makes  contact  with  the  surface,  it  exhibits  a  new  compliance 
term  because  the  very  end  of  the  tip  is  pressing  against  the 
sample  surface.  As  this  excess  compliance  results  only  when 
the  tip  is  in  contact  with  the  sample,  which  occurs  over  a 
small  fraction  of  the  pendulum-like  trajectory  of  the  fiber  tip, 
the  differential  equation  that  describes  the  tip-sample  inter¬ 
action  is  nonlinear  and  depends  on  the  tip-sample  geometry. 
An  alternative,  nonoptical  distance  regulation  scheme  de¬ 
scribed  by  Karrai  and  Grober6  involves  the  use  of  a  quartz 
crystal  tuning  fork  to  sense  the  force  exerted  on  a  fiber  tip 
glued  to  one  of  the  arms  of  the  fork.  Because  the  typical 
vibration  amplitude  of  the  fork  is  very  small  (<0.5  nm6) 
compared  to  that  required  for  optical  shear-force  detection 
(typically  tens  of  nanometers,4)  the  tip-sample  angle  would 
have  to  be  unrealistically  large  for  the  measured  crystal  ap¬ 
proach  curves6  to  be  consistent  with  the  surface  contact 
model.4  Specifically,  the  surface  contact  model  would  sug¬ 
gest  that  o'  =  tan_1(0.5  nm Jd),  where  a  is  the  tip-sample 
angle  and  d  is  the  distance  from  which  the  tip  amplitude 
starts  to  reduce  to  where  it  is  reduced  to  zero  by  complete  tip 
contact,  a  distanced  ^  10  nm  measured  in  Ref.  6  would  yield 
an  unrealistically  small  tip-sample  angle  a< 3°  (i.e.,  the  tip 
would  have  to  be  nearly  parallel  to  the  sample  surface).  A 
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more  likely  explanation  is  that  the  tuning  fork  scheme  in¬ 
volves  a  lateral  force  exerted  on  the  tip  before  tip  contact, 
possibly  caused  by  interaction  with  a  surface  contamination 
layer.  Gregor  et  aiA  acknowledged  that  over  “larger  separa¬ 
tions,”  an  interaction  before  surface  contact  may  occur,  and 
the  tuning  fork  scheme,  with  a  force  sensitivity  of  less  than  1 
nN,6  is  an  ideal  detector  to  probe  such  interactions. 

One  limitation  of  the  tuning  fork  scheme  is  its  scan 
speed.  This  results  from  the  extremely  high  Q  of  the  tuning 
fork,  which  allows  it  to  be  used  as  a  sensitive  force  detector 
despite  its  stiff  spring  constant  (typically  20  000  N/m),  but 
limits  the  transient  amplitude  response  to  a  time  constant  of 
T^lQyJl/ajQ,  where  Q  and  o>0  are  the  measured  quality 
factor  (at  full  width  half  maximum)  and  resonance  frequency 
of  the  crystal  with  the  fiber  attached,  respectively.  For  typical 
values  of  2"  1000  and  o>0~  (2tt)33  000  rad/s,  this  yields  a 
time  constant  of  r—  17  ms,  which  implies  a  rather  large  set¬ 
tling  time  (to  1%)  of  5r=85  ms.  In  a  conventional  servo 
loop,  one  cannot  increase  the  bandwidth  much  beyond  that 
of  the  slowest  system  element  (either  actuator  or  sensor). 
Furthermore,  if  the  fork  were  operated  in  a  vacuum,  then 
Q  would  typically  increase  by  a  factor  of  10,  yielding  an 
unreasonably  large  settling  time  that  would  make  the  maxi¬ 
mum  scanning  speed  unacceptably  slow  for  most 
applications.7 

In  order  to  overcome  this  limitation,  we  have  imple¬ 
mented  a  frequency  modulation/demodulation  scheme  simi¬ 
lar  to  that  used  for  high-0  cantilevers  in  atomic  force 
microscopy.7  Conceptually,  if  the  shear  force  mechanism  is 
not  purely  dissipative,  as  suggested  previously,6  but  also  in¬ 
cludes  a  conservative  component  that  changes  the  effective 
compliance  of  the  cantilever,  then  the  resonance  frequency 
of  the  cantilever  will  change  instantaneously  independent  of 
the  cantilever  Q  as  the  tip  approaches  the  sample.  A  higher 
Q  would  simply  increase  the  signal-to-noise  ratio  ( S/N )  for 
observation  of  the  resonance  frequency.7  A  significant  com¬ 
pliance  change  would  be  expected  if  the  shear  force  mecha¬ 
nism  results  from  the  very  end  of  the  tip  pressing  against  the 
sample,4  or  against  the  sample’s  contamination  layer.  In  or¬ 
der  to  facilitate  frequency  modulation,  the  tuning  fork  was 
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Distance  Approached  (nm) 

FIG.  1 .  Simultaneous  amplitude  and  frequency  deviation  of  the  quartz  crys¬ 
tal  oscillator  upon  approach  to  the  sample.  Within  6  nm  from  the  interaction 
area,  the  amplitude  decays  by  over  20%  while  the  frequency  increases  by 
2.3  Hz.  The  center  frequency  /0=  33,  089  kHz  with  £  =  2100. 

driven  at  its  resonance  frequency  by  sampling  the  current 
through  it  and  positively  feeding  back  this  signal  into  the 
crystal,  as  is  commonly  done  in  conventional  crystal  oscilla¬ 
tor  circuits.8  As  the  resonance  frequency  of  the  cantilever  is 
(ol= k/m,  where  k  is  its  effective  spring  constant  and  m  is 
the  effective  mass  of  the  cantilever,  any  variations  of  k  or 
m  will  result  in  an  instantaneous  resonance  frequency 
change.  If  we  consider  the  amplitude  reduction  upon  ap¬ 
proach  to  result  from  a  change  in  viscosity,6  then  the  reso¬ 
nance  frequency,  within  the  framework  of  a  simple  harmonic 
oscillator  model,  would  decrease  slightly  given  the  initial 
large  Q.  However,  as  plotted  in  Fig.  1,  we  have  experimen¬ 
tally  observed  an  increase  in  the  resonance  frequency  of  the 
oscillator  upon  approach  of  over  2  Hz,  suggesting  a  phenom¬ 
enological  repulsive  lateral  force  gradient  dFIdx^l  N/m  for 
a  center  frequency  f0  =  33  kHz9  (the  tip  was  vibrating  with 
an  amplitude  of  0.47  nm  during  the  approach).  This  finding, 
together  with  the  unrealistic  tip-sample  geometry  needed  for 
tapping-induced  damping4  suggests  that  both  the  simple  har¬ 
monic  oscillator  model  and  the  direct  contact  nonlinear 
bending  force  model  fail  to  adequately  describe  the  shear 
force  mechanism. 

Because  the  frequency  shift  upon  approach  is  small 
(<0.01%),  a  simple  frequency  demodulator,  for  example  a 
phase-locked  loop  integrated  circuit  or  digital  frequency 
counter,  could  not  be  used  and  provide  the  required  band¬ 
width.  Therefore,  we  employ  the  phase  difference  between  a 
reference  signal  and  the  oscillator  signal  from  the  monitored 
current  of  the  tuning  fork  to  demodulate  the  frequency  shift 
in  a  phase-locked  loop  (PLL)  configuration.10  A  similar 
scheme  for  an  atomic  force  microscope  cantilever  was  dem¬ 
onstrated  by  Kikukawa  et  alu  Since  frequency  is  the  deriva¬ 
tive  of  phase,  a  variation  of  1  Hz  between  the  reference  and 
oscillator  frequency  manifests  itself  as  a  1  Hz  sinusoidal  sig¬ 
nal  at  the  output  of  the  phase  detector.  This  signal  readily 
serves  as  an  error  signal  in  a  PLL  arrangement  where  the 
voltage-controlled  oscillator  (VCO)  incorporates  the  piezo 
controlling  the  tip-sample  distance  of  the  fork.  A  schematic 
diagram  of  the  shear  force  PLL  is  shown  in  Fig.  2.  The  error 
signal  is  the  dc  output  of  the  phase  detector.  It  is  a  null  when 
the  tip-sample  distance  causes  the  fork  oscillation  frequency 
to  be  exactly  equal  to  that  of  the  reference  signal  and  simul- 
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FIG.  2.  Schematic  diagram  of  the  phase-locked  shear-force  microscope.  The 
current  through  the  quartz  crystal  is  converted  to  a  voltage  and  is  then 
positively  fed  back  to  the  crystal  to  form  an  oscillator  that  tracks  the  reso¬ 
nance  frequency  of  the  fiber-crystal  system.  This  frequency  shift  is  demodu¬ 
lated  by  means  of  a  phase  detector  whose  reference  input  is  set  by  the  user 
to  the  frequency  deviation  for  a  given  tip-sample  height.  The  loop  filter’s 
output  is  then  sent  into  the  z-piezo  where  a  zero  error  signal  condition  is 
obtained  when  the  tip-sample  separation  causes  the  oscillator  frequency  to 
equal  the  reference  frequency. 

taneously  90°  out  of  phase.  The  experimentally  measured 
closed-loop  frequency  response  for  a  crystal  with  <2-2100 
along  with  a  theoretical  curve  based  on  a  modeling  of  the 
phase-locked  loop  and  piezo  stage  is  shown  in  Fig.  3,  yield¬ 
ing  a  bandwidth  of  over  70  Hz  (defined  at  which  the  output 
is  attenuated  to  a  factor  of  0.707  times  the  input  relative  to 
the  zero-frequency  gain12).  This  defined  bandwidth  is  also 
inversely  proportional  to  the  loop’s  settling  time  for  the  tran¬ 
sient  to  decay  to  its  steady  state  value  to  within  1%,  thus 
facilitating  direct  comparison  with  the  amplitude  settling 
time  of  the  tuning  fork.  For  the  crystal  used  to  produce  Fig. 
3,  the  theoretical  settling  time  of  5r=  10<2\/3/tt0  =  175  ms 
(for  Q  =  2100,  o)0=  277- 33089  rad/s),  which  has  an  ampli¬ 
tude  response  bandwidth  of  5.7  Hz.  Thus,  the  experimentally 
attained  bandwidth  of  70  Hz  is  over  an  order  of  magnitude 
faster  than  that  obtainable  by  amplitude  detection.  Fig.  4 
shows  a  scan  of  a  500  nm  period  diffraction  grating  taken 


Frequency  (Hz) 

FIG.  3.  The  measured  closed-loop  frequency  response  plotted  with  a  theo¬ 
retical  curve  based  on  a  control  system  analysis  that  models  the  piezo  re¬ 
sponse,  loop  filter,  and  crystal  VCO.  The  bandwidth  is  over  70  Hz  which 
exceeds  the  limit  of  5.7  Hz  imposed  by  amplitude  detection  for  the  mea¬ 
sured  £  =  2100  by  over  an  order  of  magnitude.  The  noise  spike  near  50  Hz 
is  due  to  the  feedback  compensated  piezo  stage. 
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FIG.  4.  Shear- force  image  of  a  500  nm  period  diffraction  grating  obtained 
with  a  frequency  deviation  of  1.5  Hz.  The  tip  amplitude  was  0.47  nm  (in- 
terferometrically  measured),  the  dither  direction  is  perpendicular  to  the  grat¬ 
ing,  there  are  200  by  200  pixels,  and  the  scan  rate  is  (1/70  Hz)200=2.86 
s/line.  The  noise  in  the  image  is  less  than  3  nm  and  is  caused  primarily  by 
noise  in  the  feedback  compensated  piezo  stage. 

with  a  frequency  deviation  of  1.5  Hz.  The  noise  in  the  image 
was  determined  to  be  less  than  3  nm,  primarily  due  to  noise 
in  the  feedback-compensated  piezo  stage  (a  Melles-Griot 
Nanoblock).  The  300  Hz  bandwidth  of  this  stage  limits  the 
fastest  possible  scanning  speed. 

Besides  faster  response,  there  are  additional  benefits  of 
the  phase-locked  shear-force  microscope.  There  is  no  need 
for  an  external  dither,  resulting  in  a  more  compact  and  stable 
tip  head.  The  crystal  is  made  to  automatically  operate  at 
resonance  through  self-oscillation,  so  it  is  essentially  self¬ 
tuning  and  requires  no  time-consuming  frequency  sweeps  to 


find  the  resonance  frequency.  Because  the  scheme  is  based 
on  frequency  demodulation,  it  is  immune  from  the  usual  dc 
drifts  and  noise  associated  with  amplitude  detection.  Further¬ 
more,  the  scheme  is  nonoptical,  eliminating  the  problem  of 
the  large  background  radiation  associated  with  optical  shear- 
force  control,  and  obviating  the  need  for  alignment  of  an 
external  laser  beam  and  detector.  We  have  implemented  the 
technique  using  both  a  bare  glass  fiber  and  an  aluminum 
coated  fiber  with  no  observable  difference  in  performance. 
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A  direct-write  ablation  technique  has  been  implemented  in  a  photon  scanning  tunneling  microscope  setup. 
This  combination  allows  us  to  study  surface-plasmon  (SP)  scattering  by  in  situ  created  individual  surface 
defects,  while  the  sizes  and  shapes  of  the  defects  are  varied  continuously.  It  is  found  that  within  a  certain  range 
of  its  size,  a  “hill”  on  an  otherwise  flat  surface  can  be  the  source  of  a  very  narrow  plasmon  beam.  This  effect 
is  explained  using  the  Huygens-Fresnel  principle.  Shadowing  and  refraction  of  the  SP  field  by  smaller  defects 
has  also  been  observed.  In  order  to  explain  these  results  we  introduce  an  effective  SP  refractive  index  for  two 
classes  of  surface  defects:  shallow  topographical  defects  and  areas  covered  with  absorbed  molecular  layers. 
This  concept  allows  us  to  achieve  a  qualitative  understanding  of  plasmon  scattering  in  many  practical  cases. 
Some  simple  optical  elements  for  the  control  of  SP  propagation  are  suggested  and  demonstrated.  Our  obser¬ 
vations  suggest  numerous  practical  applications  in  multichannel  chemical  sensing,  biosensing,  and  integrated 
optics.  [SOI 63- 1829(97)05928-6] 


L  INTRODUCTION 

The  surface  plasmon  (SP)  is  a  fundamental  electromag¬ 
netic  excitation  mode  of  a  metal-dielectric  interface.1,2  The 
properties  of  SP’s  have  been  the  subject  of  extensive  studies 
because  of  their  extreme  sensitivity  to  surface  roughness  and 
the  dielectric  properties  of  surface  absorbed  layers.  This  sen¬ 
sitivity  results  from  the  exponential  dependence  of  the  SP 
electromagnetic  field  on  the  distance  from  the  interface  in 
each  of  the  bounding  media.  The  SP  is  free  to  propagate 
along  the  metal  surface  with  a  mean  free  path  determined  by 
the  quality  of  the  surface  and  Ohmic  losses.  If  the  surface  is 
perfectly  flat  then  Ohmic  loss  is  the  only  channel  for  decay, 
since  the  momentum  conservation  law  prohibits  SP  conver¬ 
sion  into  photons.  In  the  visible  range  the  SP  mean  free  path 
on  a  silver  surface  is  limited  to  about  50  /zm  but  can  extend 
to  many  mm  in  the  infrared  region.3  Similar  waves  exist  on 
the  surfaces  of  a  number  of  dielectrics  and  semiconductors. 
Since  the  properties  of  SP’s  depend  strongly  on  the  proper¬ 
ties  of  the  interface  along  which  they  propagate,  they  have 
found  applications  in  the  surface-roughness  analysis,1  gas 
detection,  and  biosensing.4 

Until  recently  there  was  no  tool  for  the  direct  study  of  SP 
propagation,  since  the  SP  field  decay  length  perpendicular  to 
the  interface  is  on  the  order  of  100  nm  in  the  visible  range. 
Indirect  information  about  SP  properties  had  been  extracted 
from  far-field  measurements  of  light  scattered  into  a  free 
space  by  SP  interaction  with  surface  defects  such  as  in 
surface-plasmon  microscopy.5  The  development  of  scanning 
probe  techniques  has  opened  the  way  to  a  study  of  SP  prop¬ 
erties  at  the  surface  along  which  the  SP  propagates  with  reso¬ 
lution  in  the  nanometer  range.  The  scanning  tunneling  mi¬ 
croscopy  (STM)  was  the  first  such  technique.  In  STM 
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experiments  an  additional  tunneling  current  has  been  ob¬ 
served  in  the  presence  of  the  SP  wave.6  A  closely  related 
technique,  in  which  the  SP  field  scattered  by  a  STM  tip  is 
detected  in  the  far-field  zone,  has  been  introduced  by  Specht 
et  al?  However,  a  metallic  STM  tip  positioned  within  a  frac¬ 
tion  of  a  nanometer  from  the  surface  causes  a  very  strong 
perturbation  of  the  initial  SP  field  distribution.  The  results  of 
these  experiments  reflect  the  properties  of  localized  plasmon 
modes8  that  develop  between  the  STM  tip  and  the  sample 
rather  than  the  unperturbed  SP  field  distribution.  Similarly  an 
atomic  force  microscope  (AFM)  tip  has  been  used  to  scatter 
the  SP  field.9  Although  less  intrusive  then  STM,  this  tech¬ 
nique  also  relies  on  the  perturbation  of  the  SP  field  as  its 
sensitive  mechanism.  Kim  et  al  even  refer  to  the  potential 
benefits  of  a  metal-coated  AFM  tip  for  obtaining  an  improve¬ 
ment  in  the  sensitivity  of  this  technique.9 

The  development  of  near-field  optical  microscopy10  has 
created  different  ways  for  probing  the  SP  field  distribution. 
For  example,  in  a  tunneling  near-field  optical  microscopy11 
light  from  a  metal-coated  optical  fiber  tip  is  used  for  SP 
excitation.  However,  the  most  successful  experimental  ge¬ 
ometry  for  measuring  the  true  unperturbed  SP  field  distribu¬ 
tion  appears  to  be  the  photon  scanning  tunneling  microscopy 
(PSTM).12  In  the  PSTM  geometry  a  fiber  tip  is  used  to  detect 
the  evanescent  field  of  the  light  being  totally  internally  re¬ 
flected  at  the  sample  surface.  In  Ref.  13,  PSTM  was  used  to 
measure  the  exponential  decay  of  the  SP  field  away  from  a 
metal-air  interface.  SP  propagation  along  the  interface  has 
also  been  studied14  and  the  SP  propagation  length  measured. 
SP’s  on  both  an  external  metal-air  interface  and  an  internal 
interface  between  a  metal  film  and  a  glass  substrate  have 
been  excited15  and  the  SP  field  distribution  for  both  modes 
has  been  mapped.  Extension  of  SP  studies  to  the  control  of 
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FIG.  I .  Schematic  view  of  the  experimental  setup.  Dispersion 
curves  for  surface  plasmons  and  photons  in  air  and  in  glass  are 
shown  in  the  inset. 

SP’s  by  two-dimensional  optics  is  anticipated.  This  could 
create  numerous  applications  in  integrated  optics  and  multi¬ 
channel  chemical  sensing  and  biosensing.  However,  until  re¬ 
cently  there  has  been  no  attempt  to  create  a  real  two- 
dimensional  optical  element  for  SP’s. 

Recently  we  studied  the  scattering  of  SP’s  by  in  situ  cre¬ 
ated  individual  surface  defects.16  Some  prototype  two- 
dimensional  optical  elements  able  to  control  SP  propagation 
were  demonstrated.  Using  a  combination  of  a  PSTM  with  a 
direct-write  lithography  technique17  we  can  create  submi¬ 
crometer  size  patterns  on  the  surface  of  a  metal  film  and 
measure  the  SP  field  redistribution  caused  by  these  patterns. 
An  uncoated  tapered  optical  fiber  allows  us  to  deliver  sub- 
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FIG.  2.  Dependence  of  the  optical  signal  on  the  tip-sample  dis¬ 
tance. 


FIG.  3.  Topography  (a)  and  optical  near-field  distribution  (b) 
around  a  crater  created  by  a  UV  laser  pulse  at  full  energy. 

stantial  light  power  into  a  diffraction-limited  area  of  the 
sample  near  the  tip.  We  are  able  to  ablate  submicrometer 
craters  on  virtually  any  opaque  surface  without  damaging  a 
transparent  fiber  tip.  On  the  other  hand,  the  resolution  of  our 
microscope  in  the  PSTM  mode  depends  strongly  on  the  di¬ 
ameter  of  the  uncoated  fiber  tip:  it  is  not  dominated  by  dif¬ 
fraction.  An  uncoated  fiber  tip  introduces  a  relatively  small 
perturbation  into  the  SP  field  distribution,  since  the  differ¬ 
ence  e.ip-  1  between  the  dielectric  constant  of  the  tip  and  a 
vacuum  is  much  smaller  than  the  difference  between  the  di¬ 
electric  constant  of  the  tip  and  a  metal  (in  the  case  of  silver 
eAg=  ~  17.9+0.7/  at  X  =  633  nm).  A  comparison  of  theoret¬ 
ical  and  experimentally  measured  SP  field  distribution  (dis¬ 
cussed  later  in  this  paper)  confirms  this  conclusion.  All  these 
facts  make  our  experimental  technique  a  very  reliable  tool 
for  a  study  of  the  optics  of  SP’s. 

Below  we  discuss  in  detail  our  experimental  setup  and 
report  experimental  observations.  In  the  discussion  section 
we  present  some  model  calculations  of  the  SP  field  distribu¬ 
tion  around  simple  defects  and  compare  them  with  our  ob¬ 
servations.  We  present  the  design  of  some  optical  elements 
that  could  control  SP  propagation.  Finally,  some  practical 
applications  of  these  elements  are  suggested. 

II.  EXPERIMENTAL  SETUP 

In  our  experiment  SP’s  are  excited  in  an  attenuated  total 
reflection  (ATR)  arrangement18  on  the  surface  of  an  80-nm- 
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FIG.  4.  Topography  (a)  and  optical  near-field  distribution  (b)  around  two  plasmon  “flashlight”  structures,  (c)  and  (d)  are  far-field  images 
taken  at  distances  1  and  4  yum  from  the  sample. 


thick  silver  film  evaporated  onto  the  hypotenuse  face  of  a 
right-angle  glass  prism  (Fig.  1).  An  unfocused  P-polarized 
HeNe  laser  beam  is  used  for  excitation  of  the  SP’s.  Wave- 
vector  matching  to  the  SP  mode  at  the  silver-air  interface  is 
achieved  when  light  is  internally  incident  on  the  silver  film  at 
an  angle  a=arcsin{[eAg/(egeAg  +  ep]l/2}  =  42°,  where  eAg 
and  eg  are  the  dielectric  constants  of  silver  and  the  glass 
prism  (this  value  stems  from  the  SP  dispersion  relation  in  the 
thick-film  limit14).  The  thickness  of  the  film  is  a  little  bit  less 
than  the  decay  length  of  the  SP  field  into  silver,  so  the  cou¬ 
pling  of  fields  between  the  top  and  bottom  interfaces  is  suf¬ 
ficient  for  effective  SP  excitation.  The  SP  local  field  is 
probed  with  an  uncoated,  adiabatically  tapered  fiber  tip  that 
is  drawn  at  the  end  of  200-yLtm  UV  fiber  by  heating  it  with  a 
CO  2  laser  beam  in  a  micropipette  puller.  The  fiber  tip  can  be 
scanned  over  the  sample  surface  with  a  constant  tip-surface 
distance  (~5  nm)  using  shear  force  feedback  as  described 
elsewhere.17  Therefore,  surface  topography  can  be  imaged 
with  a  resolution  on  the  nanometer  scale,  while  simulta¬ 
neously  allowing  the  recording  of  a  near-field  optical  image 
corresponding  to  the  distribution  of  the  SP  local-field  inten¬ 
sity  near  the  film  surface.  This  configuration  of  a  scanning 
probe  microscope  is  usually  called  a  PSTM.12 

For  the  purpose  of  lithography  the  free  end  of  the  fiber 
can  be  connected  to  the  output  of  a  248-nm  excimer  laser. 


The  laser  power  is  easily  attenuated  in  order  to  change  the 
size  of  defects  created  in  the  film.  The  use  of  uncoated  UV 
grade  tapered  fiber  and  nanosecond  pulsed  248-nm  excimer 
laser  light  allows  us  to  deliver  an  estimated  200  GW/m2  UV 
light  power  into  a  tip-sample  region  smaller  then  200X200 
nm2  during  a  single  laser  pulse.  This  power  is  sufficient  for 
local  surface  ablation  of  gold  or  silver  and  local  melting  of 
silicon  without  melting  of  the  transparent  UV  fiber  tip.17  The 
topographic  changes  caused  by  UV  irradiation  are  immedi¬ 
ately  recorded  by  the  shear  force  microscope.  The  developed 
technique  is  to  some  extent  similar  to  STM  lithography  but  is 
not  restricted  by  the  conductance  of  the  sample  surface. 

III.  EXPERIMENTAL  RESULTS 

The  quality  of  the  silver  films  used  in  our  experiments 
was  good  enough  to  provide  a  quite  homogeneous  optical 
signal  distribution  over  the  unmodified  surface  areas.  As  a 
first  step,  we  measured  the  optical  signal  dependence  on  the 
distance  between  the  tip  and  the  silver  film  (Fig.  2).  The 
signal  exhibits  an  exponential  decay  with  distance  increase: 
Lu(I)=  ~2ad  superimposed  on  a  small  background.  This 
background  signal  is  due  to  the  scattering  of  SP’s  into  pho¬ 
tons  over  the  entire  illuminated  film  area.  The  value  of  the 
exponent  {a—  1/416  nm)  is  consistent  with  the  angle  of  SP 
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FIG.  5.  More  complicated  defects  emitting  plasmon  beams:  to¬ 
pography  (a)  and  optical  near-field  distribution  (b)  around  an  indi¬ 
vidual  defect;  (c)  is  a  near-field  distribution  around  a  big  defect. 

excitation.  According  to  theory,  the  SP  field  should  decay 
into  air  with  a  decay  constant  kz  =  {k2x~  orlc2)m-  1/412 
nm,  where  kx  is  the  wave  vector  parallel  to  the  surface. 

Below  we  present  images  of  the  optical  near-field  around 
in  situ  created  defects  having  different  shapes  and  sizes. 
These  images  result  from  the  scattering  of  the  SP  field  by  the 
defects.  They  also  contain  an  evanescent  component  of  the 
totally  internally  reflected  HeNe  laser  beam  used  for  the  ex¬ 
citation  of  SP’s.  But  this  evanescent  component  does  not 
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FIG.  6.  Cross  section  of  the  optical  signal  in  the  plasmon  beam 
in  Fig.  5(c). 

propagate  along  the  surface,  so  one  can  distinguish  its  con¬ 
tribution  from  the  contribution  of  SP  field  scattering.  It  is 
important  to  mention  that  the  contrast  of  each  image  is  de¬ 
termined  internally  by  scaling  the  difference  between  the 
maximum  and  minimum  of  the  detected  signals  to  the  maxi¬ 
mum  number  of  gray  levels  in  the  image.  A  comparison 
between  images  therefore  needs  to  be  performed  with  care. 


FIG.  7.  Topography  (a)  and  optical  near-field  distribution  '(b) 
around  a  surface  defect  emitting  a  weaker  reflected  plasmon  beam. 
A  stronger  plasmon  beam  is  emitted  in  the  forward  direction. 
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FIG.  8.  Scattering  of  a  plasmon  beam  emitted  from  a  surface 
defect  in  the  top  left  comer  of  the  image  by  another  surface  defect. 
As  usual,  images  (a)  and  (b)  are  the  topography  and  the  near-field 
distribution,  respectively. 

At  maximum  laser  power  1  nJ  per  pulse  measured  in 
the  far  field  of  the  tip)  a  single  nanosecond  UV  pulse  deliv¬ 
ered  to  the  silver  film  surface  by  the  uncoated  tapered  fiber 
was  able  to  ablate  a  hole  through  the  80-nm-thick  film.  The 
topography  of  such  a  defect  [Fig.  3(a)]  looks  like  a  large 
crater  with  an  approximately  0.5-,um  hole  in  the  middle.  The 
optical  near-field  distribution  [Fig.  3(b)]  has  a  strong  maxi¬ 
mum  around  the  hole  because  the  glass  prism  surface  is  al¬ 
most  clear  of  metal  in  the  hole  and  we  see  the  near  field  of 
the  totally  internally  reflected  HeNe  laser  beam  used  for  the 
excitation  of  SP’s.  By  attenuating  the  UV  laser  pulse  power 
we  are  able  to  make  craters  of  shallower  depth  and  smaller 
diameter.  At  low  power  craters  may  even  disappear  and  we 
observe  a  “hill”  at  the  irradiated  site.  Such  hills  created  at 
approximately  f  of  the  initial  power  are  shown  in  Fig.  4. 
Strong  maxima  of  the  optical  near  field  are  again  seen 
around  these  hills.  The  surprising  feature  of  the  near-field 
image  in  Fig.  4(b)  is  the  appearance  of  quite  strong  and 
narrow  beams  (their  width  is  around  1  /mm)  shining  from  the 
hills  in  the  direction  determined  by  the  HeNe  laser  beam 
used  for  SP  excitation.  These  “beams”  are  localized  within 
the  surface.  The  near-field  intensity  seems  not  to  change  sig¬ 
nificantly  along  the  l-/mm  length  of  the  left  “beam”  visible 
in  the  image.  On  the  other  hand,  a  far-field  image  taken  7 


FIG.  9.  Topography  (a)  and  optical  near-field  distribution  (b) 
around  smaller  defects  created  at  around  {  of  the  initial  power.  The 
direction  of  plasmon  propagation  is  shown  by  the  arrows.  Shadows 
are  seen  behind  the  defects. 

/im  away  from  the  surface  is  completely  dark  on  the  same 
half-tone  scale.  A  number  of  images  taken  at  intermediate 
distances  show  the  spreading  and  exponential  attenuation  of 
the  beam  “brightness”  with  increasing  distance  from  the 
surface.  Some  of  them  are  shown  in  Figs.  4(c)-4(d).  (Note, 
however,  that  the  contrast  for  both  images  is  internal  and  not 
relative  to  each  other.)  Evidently,  in  the  near-field  image 
[Fig.  4(b)]  we  see  surface-plasmon  beams  shining  from  the 
peaks  of  the  optical  near-field  around  the  hills.  We  like  to 
call  these  structures  plasmon  “flashlights.”  Our  experience 
shows  that  the  appearance  of  a  plasmon  beam  is  a  rather 
common  features  of  an  optical -field  distribution  near  surface 
defects  in  the  ATR  geometry.  In  Fig.  5  one  can  see  some 
other  images  showing  plasmon  beams.  If  the  defect  has  a 
complicated  shape  as  in  Fig.  5(c)  the  plasmon  beam  can 
shine  from  some  region  of  the  defect.  The  cross  section  (Fig. 
6)  of  the  plasmon  beam  in  this  image  shows  the  lateral  ex¬ 
ponential  decay  of  the  SP  field  intensity  with  the  distance 
from  the  source  of  the  beam.  The  value  of  the  SP  propaga- 
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FIG.  10.  Topography  (a)  and  optical  near-field  distribution  (b) 
around  three  surface  defects. 


lion  length  equals  9.1  fi m  in  this  case.  This  value  is  close  to 
the  I3.2-yiun  value  reported  in  Ref.  14  yet  is  much  smaller 
than  the  theoretically  estimated  propagation  length  42.4 
Aim.14  A  possible  explanation  of  this  discrepancy  is  that  the 
image  in  Fig.  5(c)  was  recorded  using  a  rather  old  silver  film 
sample  whose  surface  looked  slightly  oxidized.  The  value  of 
the  propagation  length  derived  from  the  image  in  Fig.  4(b) 
taken  using  a  freshly  prepared  silver  film  equals  30  yam.  This 
number  is  much  closer  to  the  theoretical  value. 

In  some  images  of  surface  defects  emitting  SP  beams  one 
can  see  not  only  a  forward  beam  but  also  a  weaker  beam 
emitted  in  the  backward  direction  (Fig.  7).  This  means  that  a 
“plasmon  mirror”  can  be  created  as  another  prototype  sur¬ 
face  optical  element.  A  plasmon  beam  from  one  defect  can 
be  directed  onto  another  one,  as  is  the  case  in  Fig.  8.  One  can 
see  a  large  defect  in  the  top  left  corner  of  the  image  that 
emits  a  wide  plasmon  beam.  This  beam  is  scattered  by  an¬ 
other  defect  in  the  middle  of  the  field  of  view.  This  experi¬ 
ment  shows  that  one  can  deal  with  plasmon  “flashlights”  in 
the  manner  that  is  very  similar  to  the  operation  of  usual 


three-dimensional  light  sources. 

Attenuation  of  the  UV  power  during  the  lithography  pro¬ 
cess  allows  us  to  create  smaller  surface  defects  and  study 
plasmon  scattering  by  these  defects.  The  results  of  these 
studies  are  shown  in  Figs.  9  and  10.  The  UV  power  used  to 
create  these  smaller  defects  was  approximately  {  of  the  ini¬ 
tial  power.  We  cannot  see  bright  near-field  maxima  around 
these  defects  any  more.  Instead,  we  see  prominent  shadows 
behind  the  defects  (the  direction  of  plasmon  propagation  is 
shown  by  the  arrows).  Two  characteristic  features  of  these 
figures  emerge.  First,  the  sources  of  the  shadows  appear  to  be 
larger  in  the  near-held  images  than  in  the  topographic  ones 
(a  comparison  of  these  images  with  the  near-field  image  in 
Fig.  4(b)  shows  that  this  effect  does  not  stem  from  the  reso¬ 
lution  of  the  near- field  images).  Second,  for  different  defects 
the  shadows  have  different,  but  still  large,  angular  widths. 
The  approximately  90° -angular  widths  cannot  be  accounted 
for  by  diffraction.  The  diffraction  angle  can  be  easily  esti¬ 
mated  from  the  uncertainty  principle: 
/3~£Ap|/£p|~\pi/(2  77Y/),  where  \p]  =  590  nm  is  a  plasmon 
wavelength  and  d  is  the  size  of  a  scatterer.  A  micrometer 
size  scatterer  gives  diffraction  angle  values  on  the  order  of 
5°.  This  value  is  consistent  with  the  rather  small  visible 
width  of  the  SP  beams  from  the  plasmon  “flashlights”  in 
Fig.  4(b).  Another  interesting  feature  is  an  apparent  differ¬ 
ence  in  the  ability  to  scatter  SP  field  at  hills  and  dips  of 
similar  geometrical  size.  Looking  at  Fig.  10  one  can  notice 
three  defects  of  about  the  same  lateral  size.  Two  of  them  are 
shallow  hills.  The  third  is  a  dip  of  about  the  same  vertical 
size.  The  hills  seem  not  to  affect  the  SP  propagation  at  all. 
On  the  other  hand,  the  dip  gives  a  distinct  angular  shadow.  In 
the  next  part  of  the  paper  we  propose  an  explanation  for 
these  phenomena.  We  have  used  these  phenomena  to  create 
another  simple  two-dimensional  optical  element  which  we 
call  a  “beam  stop.”16  A  number  of  defects  were  created 
along  a  single  line  on  the  surface  of  a  silver  film  (Fig.  4  in 
Ref.  16).  The  near-field  image  shows  some  increase  in  the 
optical  field  around  the  defects  and  a  prominent  dark  shadow 
behind  the  line  of  defects.  Thus,  an  effective  optical  shield 
has  been  created  that  does  not  allow  plasmons  to  pass  by. 


IV.  DISCUSSION 

The  scattering  of  SP’s  by  surface  defects  has  been  con¬ 
sidered  in  numerous  theoretical  papers.  A  recent  paper  by 
Pincemin  et  al. 19  contains  a  good  overview  of  the  research  in 
this  field.  Generally  this  problem  has  no  analytical  solution. 
Numerically  it  has  been  solved  for  a  limited  number  of  ex¬ 
perimental  geometries.  No  comparison  between  such  calcu¬ 
lations  and  experimental  near-field  data  has  been  carried  out. 
Unfortunately,  we  have  not  found  in  the  literature  any  result 
of  numerical  calculations  that  could  be  directly  applied  to 
our  observations.  Nevertheless  we  propose  two  simple  mod¬ 
els  that  allow  us  to  understand  our  results  and  compare  them 
with  theoretical  predictions. 

The  plasmon  flashlight  phenomenon  can  be  understood 
directly  from  the  Huygens-Fresnel  principle.  Let  us  consider 
an  array  of  ten  secondary  sources  of  SP  circular  waves  (Fig. 
11)  equally  spaced  along  a  2-/zm  line.  The  resulting  SP  field 
distribution  is  given  by  the  expression 
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FIG.  I  1.  Field  distribution  from  a  2-/xm-long  array  of  point  sources  in  two-dimensional  (a),  (b)  and  three-dimensional  (c),(d)  cases.  A 
pronounced  SP  beam  appears  for  a  two-dimensional  field  distribution.  Horizontal  size  is  10X  10  fx m.  Halftone  views  and  line  scan  views  are 
presented. 


eMrn  +  *n)lr\nJ  5  (j) 

where  rn  is  the  distance  between  the  point  of  observation  and 
the  nth  secondary  source.  This  SP  field  distribution  is  shown 
in  Fig.  1 1  (a),  1 1  (b)  for  the  case  where  the  phases  of  second¬ 
ary  sources  change  linearly  with  n.  A  plasmon  beam  similar 
to  the  one  in  Fig.  5(b)  is  clearly  visible  in  the  half-tone 
representation.  It  appears  that  the  primary  reason  for  the 
beam  existence  is  a  weak  1/2  distance  dependence  of  the 
fields  from  individual  point  sources  of  circular  SP  waves. 
For  a  comparison,  the  field  distribution  from  the  same  source 
calculated  using  the  bulk  r~]  distance  dependence  is  shown 
in  Figs.  11(c),  11(d).  No  prominent  beam  appears  in  this 
case.  We  believe  that  the  strong  maxima  of  the  optical  near 
field  observed  in  Figs.  4(b)  and  5(b)  can  be  represented  by 
the  arrays  of  secondary  sources  of  SP  circular  waves.  The 
phase  shift  between  the  secondary  waves  is  due  to  the  angu¬ 
lar  illumination  of  the  defects  by  the  He-Ne  excitation  laser. 
It  results  in  the  plasmon  beam  shining  in  the  direction  of  the 
incoming  beam  of  HeNe  laser.  The  enhancement  of  the  near 
field  around  the  defects  should  be  calculated  from  the  exact 
theory.  The  main  factors  contributing  to  the  enhancement  are 
a  ‘‘lightning  rod  effect,”  which  leads  to  field  enhancement 


near  the  prominent  hills  on  the  silver  film,  and  larger  photon 
tunneling  through  the  thinner  regions  of  the  film. 

Figure  12  represents  our  modeling  of  data  shown  in  Fig. 
4.  Two  arrays  of  secondary  sources  similar  to  the  one  dis¬ 
cussed  in  the  previous  paragraph  are  positioned  at  different 
distances  from  each  other.  Figure  1 2(c)  demonstrates  a  very 
close  resemblance  to  Fig.  4(b).  It  is  really  a  surprising  result 
that  two  plasmon  flashlight  structures  can  be  positioned 
within  a  few  micrometers  from  each  other  giving  rise  to  dis¬ 
tinctly  separate  plasmon  beams.  One  can  envision  numerous 
applications  of  these  structures  in  sensor  arrays.  Each  plas¬ 
mon  beam  could  probe  the  surface  chemistry  in  a  very  local¬ 
ized  region.  Based  on  the  success  of  our  plasmon  beam  mod¬ 
eling  we  suggest  some  optical  elements  that  could  control  SP 
propagation.  One  could  think  of  focusing  mirrors  or  lenses 
consisting  of  surface  defects  similar  to  ones  in  Fig.  4.  Some 
examples  of  these  optical  elements  are  shown  in  Fig.  13.  A 
curved  defect  described  (in  polar  coordinates)  by  the  expres¬ 
sion  R(a)  =  2F/(  1  +cosa)  would  be  a  perfect  focusing  mir¬ 
ror  with  focal  length  F:  R(a)  is  a  solution  of  the  equation 
kR(a)  +  <£(a)  =  const  [where  <f>(  a)  =  kR(a)cos(a)]  for  the 
positions  of  the  point  sources  that  would  direct  reflected  sec¬ 
ondary  SP  waves  to  be  exactly  in  phase  at  some  particular 
focal  point.  A  curved  defect  as  shown  in  Fig.  13(b),  which 
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FIG.  12.  Numerical  modeling  of  two  plasmon  flashlights  at  different  distances  from  each  other.  Horizontal  size  is  10X  10  /im. 


consists  of  two  parts  of  this  focusing  mirror  positioned  next 
to  each  other  should  have  two  focal  points.  These  calcula¬ 
tions  are  justified  by  the  observation  of  a  “plasmon  mirror” 
in  Fig.  7.  Another  possibility  is  to  create  a  focusing  lens 
using  a  zone  plate  principle.  Figure  13(c)  represents  the  field 
distribution  of  such  an  element.  The  defects  are  positioned 
on  a  circle  of  radius  R  in  order  to  compensate  for  the  dis¬ 
tance  dependence  of  the  field  in  a  circular  wave.  The  angular 
position  of  the  // th  point  SP  source  is  given  by  the  expression 
kR-kRcosla„)  =  2im.  All  these  optical  elements  could  find 
numerous  applications  in  integrated  optics. 

Let  us  now  discuss  our  data  concerning  the  SP  scattering 
by  smaller  defects.  The  simplest  explanation  of  the  observed 
discrepancy  between  the  visible  sizes  of  the  shadow  sources 
in  the  near-field  and  topographic  maps  could  be  that  the  ac¬ 
tual  size  of  the  scatterers  is  not  the  one  that  is  seen  in  the 
surface  topographic  map.  The  bulk  structure  of  the  silver 
film  can  be  deformed  around  the  surface  defect.  Also  the 
chemical  composition  of  the  surface  oxide  layer  (and/or  ab¬ 
sorbed  layer)  could  he  changed  under  UV  irradiation  on  a 
much  larger  scale  than  is  apparent  in  the  visible  size  of  the 
defect.  Precise  control  of  these  parameters  would  require  ex¬ 
periments  under  ullrniiigh-vacutim  conditions.  But,  there  is 
strong  theoretical  evidence  that  the  visible  size  of  even  a 
point  scatterer  in  the  near-field  map  may  he  on  the  order  of 
some  micrometers.  Let  us  imagine  a  point  defect  on  an  oth¬ 
erwise  ideally  flat  metal  film  surface.  This  defect  is  irradiated 
by  a  plane  surface-plasmon  wave.  The  defect  breaks  transla¬ 
tion  symmetry  of  the  film  surface.  The  SP  wave  can  transfer 
part  of  its  momentum  to  the  defect  and  dccn>  into  a  photon. 


Consider  the  cross  section  of  the  point  defect  for  such  a 
process.  This  point  defect  can  be  Fourier  decomposed  into  a 
set  of  diffraction  gratings  with  periods  from  zero  to  infinity 
(in  reality,  the  upper  limit  is  the  inverse  atomic  size).  Since 
SP  scattering  is  a  linear  process  we  can  treat  the  SP  interac¬ 
tion  with  each  diffraction  grating  separately.  Each  grating 
can  change  the  momentum  of  a  plasmon  by  a  multiple  of  its 
inverse  period.  Most  probably  the  plasmon  momentum  will 
change  by  only  one  inverse  grating  period.  If  the  inverse 
grating  period  is  smaller  than  the  difference  between  the  mo¬ 
mentum  of  plasmon  and  the  momentum  of  photon,  then  plas¬ 
mon  decay  will  not  happen.  Hence,  only  grating  with  suffi¬ 
ciently  large  inverse  periods  will  contribute  to  the  cross 
section  of  SP  decay  into  photons.  In  order  to  find  this  cross 
section  we  should  cutoff  the  Fourier  spectrum  of  the  defect  at 
wave  vectors  less  then  the  difference  between  the  plasmon 
and  photon  wave  vectors  (&pi— £ph)  and  perform  an  inverse 
Fourier  transformation.  Such  a  procedure  gives  an  effective 
cross  section  s~  V(kpl-kph).  This  cross  section  is  s—1.6 
yum  at  HeNe  laser  wavelengths.  This  value  is  consistent  with 
the  effective  sizes  of  the  shadow  sources  in  Fig.  9(b). 

Another  effect  that  might  produce  a  shadowlike  phenom¬ 
enon  around  a  surface  defect  is  attenuated  energy  transfer 
from  the  He-Ne  laser  excitation  beam  into  the  SP’s.  Reso¬ 
nant  excitation  of  SP’s  is  happening  because  of  the  phase 
matching  between  the  SP  and  the  HeNe  beam  in  the  glass 
prism.  The  phase-matching  length  L  is  determined  by  the 
angular  width  of  the  SP  resonance  (see  Fig.  1)  which  is  usu¬ 
ally  on  the  order  of  da-  1  °:  SkplL~  Sakp]L - 1 .  This  gives 
a  value  of  the  phase-matching  length  on  the  order  of  L~5 
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FIG.  13.  Numerical  modeling  of  plasmon  mirrors  (a),(b)  and  a  plasmon  lens  (c).  Horizontal  size  is  10X  10  /nm. 


/xm.  A  surface  defect  could  worsen  this  phase  matching  and 
cause  smaller  energy  transfer  in  to  the  SP’s.  A  smaller  inten¬ 
sity  of  the  SP  field  would  look  like  a  shadow  behind  the 
defect.  The  two  effects  discussed  above  lead  to  essentially 
the  same  effective  size  of  the  shadow  source.  We  believe 
both  effects  contribute  to  the  observed  field  distribution.  Un¬ 
fortunately,  neither  of  these  effects  could  be  responsible  for 
the  angular  shape  of  the  shadows  and  an  apparent  difference 
between  SP  scattering  by  shallow  hills  and  dips. 

The  most  probable  reason  for  shadows  to  have  an  angular 
shape  is  refraction.  As  we  have  discussed  earlier,  the  SP 
interaction  with  a  surface  defect  can  be  decomposed  into 
interactions  with  separate  gratings  that  constitute  the  Fourier 
spectrum  of  a  defect.  Interaction  with  gratings  that  have  pe¬ 
riods  much  smaller  than  2irl(kpl-kph)  leads  to  SP  decay 
into  photons.  Interaction  with  the  rest  of  the  gratings  can  be 
treated  like  two-dimensional  optics:  a  plasmon  will  remain  a 
plasmon  after  the  interaction.  Thus,  we  can  talk  about  SP 
refraction. 

It  is  relatively  easy  to  introduce  an  effective  refractive 
index  for  a  shallow  dip  or  hill  that  scatters  a  SP  wave.  The 


dispersion  relation  for  the  SP  in  the  three-media  air-metal- 
glass  system  can  be  written  as  follows:20 


(ej  !kx  4-  e2/k2)(€2/k2+  e2lk3) 

"P (^i !k\ “ €2 lk2)( €fk2  —  £3 !k^)e  (2) 

where  kn  =  (k2- a)2en  !c2) 1/2,  and  n=  1,2,3  denotes  glass, 
silver  film,  and  air,  respectively.  After  some  transformations 
the  dispersion  law  can  be  rewritten  in  the  more  explicit  form: 


u)=zkxc[{€2+\)l€2\m{\+Ae  2*2^2),  (3) 


where 


A  —  2e2 


[(-61f2+f2-6|)1/2+/61] 

{(*2~1  )[(-*,  *2  +e2-el)m+i€  1]} 


— 0.11 


(4) 
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effective  SP  refractive 
index  changes: 


FIG.  14.  Summary  of  different  physical  mechanisms  affecting 
surface-plasmon  field  distribution. 

(here  we  assume  that  the  thickness  of  the  silver  film  is  large 
enough  that  e~2kld2<  1).  As  we  can  see,  the  phase  velocity 
of  SP  depends  on  the  thickness  of  the  silver  film.  The  SP’s 
are  faster  in  the  areas  where  the  thickness  of  a  film  is 
smaller.  The  effective  refractive  index  of  a  shallow  defect 
can  be  written  as  follows: 

n=VQ/V=l+A(e~2k^o~e~2k^),  (5) 

where  V  is  a  SP  phase  velocity  at  thickness  d  and  V0  is  a  SP 
phase  velocity  in  unmodified  areas  of  a  silver  film  with  a 
thickness  d0.  Numerical  analysis  of  this  expression  shows 
that  in  our  experiments  (with  d0=80  nm)  the  effective  re¬ 
fractive  index  n  varies  within  the  range  0.89- 1.0001.  For 
shallow  hills  n>  1  and  the  changes  in  n  are  negligible,  while 
shallow  dips  have  n<  1.  This  results  agrees  qualitatively 
with  our  observations  in  Fig.  10.  Shallow  hills  should  have 
much  less  of  an  effect  on  the  SP  propagation  than  shallow 
dips.  Evidently,  shallow  dips  act  as  negative  lenses.  This 
results  in  the  angular  shape  of  shadows. 

Similar  effective  refractive  indices  could  be  assigned  to 
areas  of  the  silver  film  where  the  composition  of  an  absorbed 
layer  is  modified.  Sometimes  the  presence  of  an  absorbed 
layer  could  lead  to  appearance  of  additional  SP  modes.21  But 
in  many  cases  the  presence  of  a  few  molecular  monolayers 
on  the  surface  of  a  metal  film  leads  only  to  a  shift  in  the  SP 
phase  velocity.22  In  terms  of  effective  refractive  index  this 
shift  can  be  expressed  as  follows: 

n=1  +  (€a-l)/[ea(-eAg-l)''2]X2iTda/X,  (6) 

where  ea  is  a  dielectric  constant  and  da  is  the  thickness  of  an 
absorbed  layer  (here  we  assume  that  efl<^eAg).  According  to 
Ref.  23  four  monolayers  of  cadmium  arachidate  gives  a  shift 
in  the  SP  phase  velocity  which  corresponds  to  n  =  1.033.  The 
potential  consequence  of  this  result  for  our  experiment  is  that 
possible  modification  of  an  absorbed  layer  by  an  UV  light 
pulse  during  defect  creation  could  cause  refraction  of  SP’s 
around  the  defect.  On  the  other  hand,  this  effect  can  be  used 
to  create  simple  optical  elements  for  SP  field  redistribution 
such  as  lenses  or  prisms.  A  triangular  region  of  the  silver 
film  covered  with  a  layer  of  molecules  such  as  cadmium 
arachidate  would  act  like  a  prism  on  a  propagating  SP  beam. 
One  could  calculate  the  dispersion  of  this  prism  by  taking 
into  account  the  dispersion  of  eAg  and  ea .  This  has  the  same 


order  of  magnitude  as  the  dispersion  of  any  bulk  prism.  The 
dispersion  of  such  optical  elements  can  be  used  in  multichan¬ 
nel  chemical  sensors  and  biosensors. 

Figure  14  summarizes  our  discussion  of  different  effects 
that  influence  the  SP  field  distribution  near  surface  defects.  A 
change  of  metal  film  thickness  and  the  local  surface  tilt  cause 
the  modification  of  the  phase  matching  between  the  He-Ne 
laser  light  and  SP’s.  Some  sharp  surface  “hills”  and  holes 
through  the  film  can  cause  strong-field  enhancement  in  their 
vicinity.  These  effects  influence  the  electromagnetic  field 
distribution  in  close  proximity  to  these  defects.  One  needs  to 
use  an  exact  theory  to  describe  these  “near-field”  effects. 

The  SP  “far-field”  distribution  is  much  easier  to  under¬ 
stand  using  such  simple  concepts  as  the  Huygens-Fresnel 
principle  and  effective  SP  refractive  index.  One  can  replace 
strong  maxima  of  the  near  field  around  the  hills  and  the  holes 
with  secondary  SP  sources.  Adiabatic  changes  in  film  thick¬ 
ness  and  local  variations  in  the  absorbed  molecular  layer  can 
be  described  by  an  effective  refractive  indices.  Then  the  na¬ 
ture  of  SP  propagation  can  be  treated  as  a  two-dimensional 
optics  problem. 


V.  CONCLUSIONS 

In  summary,  we  have  constructed  a  combined  PSTM/ 
direct- wire  lithography  setup  that  has  allowed  us  to  study  the 
scattering  of  surface  plasmons  by  individual  defects,  while 
continuously  adjusting  their  shapes  and  sizes.  The  use  of  an 
uncoated  fiber  tip  as  a  local  probe  causes  much  less  pertur¬ 
bation  of  the  SP  field  as  compared  to  other  scanning  probe 
techniques  such  as  STM,  AFM,  or  near-field  optical  micros¬ 
copy  with  metal-coated  tips. 

Surface  defects  that  emit  micrometer  wide  SP  beams 
(plasmon  flashlights)  have  been  observed.  This  phenomenon 
has  been  explained  using  the  Huygens-Fresnel  principle. 
Scattering  and  refraction  of  SP  field  by  smaller  surface  de¬ 
fects  has  also  been  studied.  The  concept  of  effective  refrac¬ 
tive  index  has  been  introduced  for  two  classes  of  surface 
defects:  shallow  hills  or  dips,  and  areas  covered  with  a  layer 
of  absorbed  molecules.  This  concept  allows  us  to  achieve  a 
qualitative  understanding  of  SP  field  redistribution  by  de¬ 
fects,  which  is  otherwise  a  complicated  theoretical  problem 
without  an  analytical  solution. 

Some  simple  optical  elements  able  to  govern  SP  field 
propagation  have  been  suggested.  We  believe  our  technique 
has  the  potential  to  create  analogs  of  any  three-dimensional 
optical  devices  in  two  dimensions.  Similar  techniques  can  be 
used  not  only  in  the  two-dimensional  optics  of  surface  plas¬ 
mons  on  metal  surfaces  but  also  in  the  optics  of  any  surface 
waves  on  semiconductors  and  dielectrics.  Plasmon  flash¬ 
lights  in  combination  with  other  SP  optical  elements  could 
have  numerous  potential  applications  in  multichannel  chemi¬ 
cal  sensing,  biosensing,  and  integrated  optics. 
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Local  second-harmonic  generation  (SHG)  from  a  rough  metal  surface  has  been  measured  using  a  near-field 
optical  microscope.  The  dependence  of  the  SH  signal  on  the  tip-surface  distance  for  excitation  with  S -  and 
P -polarized  light  shows  the  presence  of  evanescent  SH  field  components.  With  excitation  by  S-polarized  light 
the  SHG  over  the  rough  surface  is  found  to  be  strongly  related  to  surface  topography.  For  P -polarized  light 
excitation,  local  enhancement  of  SH  emission  have  been  observed  that  may  be  attributed  to  SH  field  localiza¬ 
tion.  [S0163- 1829(97)07040-9] 


I.  INTRODUCTION 

In  recent  years  there  has  been  increasing  interest  in  linear 
and  nonlinear  optics  of  disordered  media,  particularly  rough 
surfaces.1  Because  second-harmonic  generation  (SHG)  is 
known  to  be  surface  sensitive  on  an  atomic  scale2  (especially 
when  the  nonlinearity  itself  is  due  to  the  presence  of  the 
surface,  as  is  the  case  for  isotropic  materials*),  studies  of  the 
nonlinear  optical  response  are  particularly  suitable  for  im¬ 
proving  our  knowledge  of  the  relationship  between  optical 
properties  and  morphology  (roughness,  defects,  impurities, 
adsorbates)  of  metallic  and  semiconductor  surfaces.  ’ 5 

Second-harmonic  generation  from  very  flat  surfaces  has 
been  extensively  studied.2,6  However,  real  surfaces  always 
have  different  degrees  of  roughness.  The  inclusion  of  rough¬ 
ness  leads  to  several  phenomena  affecting  SHG  (such  as  re¬ 
laxation  of  polarization  selection  rules  and  an  enhancement 
of  the  electromagnetic  field  by  the  roughness).  Less  obvious 
is  the  strong  influence  of  weak-localization  effects  on  the 
SHG  from  rough  surfaces7  that  result  in  significant  changes 
of  the  angular  spectrum  of  SHG  observed  in  the  far-field 
region.8  Even  more  pronounced  changes  should  be  expected 
in  the  near-field  region  close  to  the  surface. 

Until  now,  SHG  from  rough  surfaces  has  been  studied 
only  in  the  far-field  region  without  exact  knowledge  of  the 
surface  topography.  Information  provided  by  the  usual  (far- 
field)  optical  studies  of  surfaces  is  derived  from  data  aver¬ 
aged  over  a  surface  region  limited,  in  the  best  case,  by  dif¬ 
fraction  of  the  probe  light.  However,  the  usual  sizes  of  these 
regions  is  much  larger:  from  several  square  micrometers  to 
several  square  millimeters.  Thus,  optical  data  obtained  from 
linear  and  nonlinear  spectroscopy  are  valid  for  very  large 
defect  ensembles.  At  the  same  time,  the  lateral  distribution  of 
the  electromagnetic  field  over  a  surface  is  not  uniform  and 
depends  itself  on  the  surface  defect  structure.  The  local-field 
intensity  can  vary  by  several  orders  of  magnitude  on  a  scale 
less  than  half  a  wavelength  along  the  surface.9,10  Therefore, 
investigation  of  the  averaged  optical  response  in  many  cases 
does  not  result  in  an  understanding  of  the  underlying  physics 
(especially  in  nonlinear  spectroscopy  where  optical  response 
depends  on  the  driving  field  in  a  nonlinear  manner). 

The  recent  development  of  near-field  optical  microscopy 
(NFOM)  has  opened  the  possibility  for  study  of  numerous 
optical  phenomena  with  a  resolution  well  below  the  diffrac- 
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tion  limit.11  Most  near-field  optical  studies  have  concentrated 
on  the  investigation  of  the  linear  optical  response  of  materi¬ 
als.  The  possibility  of  studying  nonlinear  optical  processes 
(especially,  wave  mixing)  with  near-field  optical  character¬ 
ization  methods  has  been  discussed  recently,12,13  but  only 
phase-conjugation  experiments  were  experimentally  realized 
in  the  near-field  region.14  This  slow  progress  in  nonlinear 
near-field  optics  stems  from  numerous  experimental  difficul¬ 
ties  caused  by  the  small  optical  throughput  of  conventional 
metal-coated  near-field  apertures.  This  small  throughput  en¬ 
ters  quadratically  into  the  measured  SH  optical  signal.  An¬ 
other  disadvantage  of  metal-coated  fiber  tips  in  nonlinear 
optical  measurements  is  a  strong  perturbation  introduced  by 
a  tip  into  the  local  surface  field  distribution.  For  example,  in 
SHG  experiments  SH  light  intensity  is  proportional  to  the 
fourth  power  of  the  local  field.  Hence,  by  using  a  metal- 
coated  fiber  tip  one  can  study  only  the  properties  of  a  mi¬ 
croresonator  formed  between  the  tip  and  the  sample,  but  not 
the  true  local  SH  field  distribution.  Uncoated  fiber  tips  (used 
in  Ref.  14)  introduce  much  smaller  perturbation  into  the  lo¬ 
cal  field  near  the  sample,  since  the  dielectric  constant  6  of  a 
fiber  is  much  less  than  e  of  a  metal.  However,  the  decrease  in 
optical  confinement  does  lead  to  somewhat  lower  resolution. 

The  combination  of  second-harmonic  generation  tech¬ 
niques  with  near-field  optical  microscopy  has  significant  po¬ 
tential  for  probing  the  nonlinear  optical  response  of  a  surface 
locally  with  subwavelength  lateral  resolution  while  simulta¬ 
neously  measuring  surface  topography  using  shear-force 
feedback.  The  data  obtained  in  such  a  way  may  allow  com¬ 
parison  of  experiment  and  theory  in  sufficient  detail  to  un¬ 
derstand  the  essential  features  of  SHG  at  rough  surfaces  and 
provide  better  understanding  of  the  underlying  microscopic 
electrodynamics.  Since  SHG  is  sensitive  to  any  asymmetry 
of  the  sample,  it  is  sensitive  to  surface  magnetic  or  electric 
dipole  moments.  Near-field  SHG  microscopy  can  be  used  to 
image  magnetic  or  electric  domains  and  is  an  alternative  to 
magnetic  and  electric  force  microscopies  in  this  regard. 
Near-field  SHG  microscopy  of  a  rough  metal  surface  is  an 
important  first  step  in  this  direction. 

In  this  paper  we  report  the  observation  of  second- 
harmonic  light  generated  at  a  rough  silver  surface  using  a 
near-field  optical  microscope.  The  SH  signal  dependence  on 
the  tip-surface  distance  has  been  measured  for  S  and  P  po¬ 
larization  of  the  excitation  light.  The  spatial  distribution  of 
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FIG.  1.  Schematic  view  of  the  near-field  optical  microscope  for 
second-harmonic  light  studies. 


SHG  over  the  surface  has  been  imaged.  Such  images  dem¬ 
onstrate  that  different  mechanisms  affect  surface  SHG  at 
rough  metal  surfaces  for  S-  and  P -polarized  excitation  light. 
Our  data  clearly  show  near-field  features  inaccessible  in  the 
usual  far-field  measurements. 


II.  EXPERIMENTAL  SETUP 


(b) 


Our  experimental  setup  for  near-field  microscopy  of 
second-harmonic  generation  is  shown  in  Fig.  1.  A  400-nm- 
thick  silver  film  was  prepared  by  thermal  evaporation  on  a 
glass  substrate  and  was  used  as  the  sample  in  our  near-field 
optical  microscope.  The  local  SH  field  distribution  has  been 
probed  with  an  uncoated  adiabatically  tapered  fiber  tip, 
which  is  drawn  at  the  end  of  a  single  mode  fiber  by  heating 
it  with  a  C02  laser  in  a  micropipette  puller.  The  fiber  tip  can 
be  scanned  over  the  sample  surface  with  a  constant  tip- 
surface  distance  of  a  few  nanometers  using  shear  force  feed¬ 
back  control.  Therefore,  surface  topography  can  be  imaged 
with  a  resolution  on  a  nanometer  scale,  while  simultaneously 
recording  the  SH  near-field  image.  The  SHG  has  been  ex¬ 
cited  at  an  angle  of  incidence  of  60°  with  a  Nd:  YAG  (yttrium 
aluminum  garnet)  laser  operating  at  1064  nm  (repetition  rate 
10  Hz,  pulse  duration  20  ns,  pulse  energy  10-15  mJ).  The 
excitation  power  at  the  sample  surface  is  estimated  to  be 
about  1  MW/cm2.  The  SH  signal  has  been  measured  with  a 
photomultiplier  and  gated  electronics.  The  SH  signal  at  every 
point  of  the  image  has  beert  averaged  over  70  or  100  laser 
pulses  for  P-polarized  and  S-polarized  light  excitation,  re¬ 
spectively.  The  characteristic  SH  photon  counting  rate  was 
on  the  order  of  one  photon  count  per  3-5  laser  pulses.  The 
main  source  of  noise  in  our  data  is  statistical  fluctuation 
(Nm)  of  the  number  of  photon  counts  detected.  The  mea¬ 
sured  SH  intensity  depends  quadratically  on  the  fundamental 
light  power.  The  spatial  resolution  of  the  microscope  in  the 
SH  light  collection  mode  has  been  determined  to  be  better 
than  150  nm  (this  is  the  size  of  the  smallest  features  visible 
in  the  SH  images  of  piezoelectric  ceramics  and  the  sharpness 
of  magnetic  domain  walls  observed  in  SH  images  of  a  Ni 
surface16).  The  fiber  tip  itself  did  not  contribute  to  the  SH 


FIG.  2.  (a)  Variation  of  SH  signal  on  the  tip-sample  distance  for 
S-  and  P -polarized  excitation.  These  variations  have  been  measured 
in  the  center  of  the  silver  film  region  imaged  in  (b).  The  noise  level 
in  these  data  is  determined  by  the  statistic  fluctuations  (Ny2). 
Dashed  lines  are  simple  power-law  curves  drawn  as  a  guide  to  the 
eye. 

signal:  we  measured  zero  SH  signal  when  the  sample  was  far 
from  the  tip. 

III.  RESULTS  AND  DISCUSSION 

To  determine  the  contribution  of  near-field  processes  to 
the  SH  generation,  we  have  studied  the  dependence  of  the 
second-harmonic  signal  on  the  tip-sample  distance  for  5-  and 
P -polarized  excitation  light  [Fig.  2(a)].  These  measurements 
were  made  near  the  center  of  the  silver  film  surface  shown  in 
Fig.  2(b).  A  rather  strong  decrease  in  the  signal  occurs  about 
500  nm  from  the  surface.  This  behavior  is  especially  pro¬ 
nounced  for  S-polarized  excitation.  This  initial  drop  in  the 
signal  is  followed  by  a  number  of  oscillations  for  both  po¬ 
larizations.  Note  that  each  maximum  for  P-polarized  excita¬ 
tion  corresponds  to  a  minimum  for  S-polarized  excitation 
and  vice  versa.  Similar  oscillations  have  been  observed  at 
every  point  on  the  sample  surface.  We  believe  that  this  be¬ 
havior  results  from  the  standing  waves  of  the  1064  nm  fun¬ 
damental  light  formed  between  the  glass  tip  and  silver  film 
surface.  The  phase  shift  of  rr  between  the  distance  dependen¬ 
cies  can  be  explained  by  the  phase  shift  of  tt  between  P-  and 
^-polarized  fundamental  waves  reflected  between  tip  and 
metal  at  angles  larger  than  the  Brewster  angle. 

The  ratio  of  the  SH  signals  induced  by  P-  and  S-polarized 
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light  changes  substantially  in  the  near-field  region.  This  ratio 
hills  from  3.1  at  a  distance  of  3600  nm  from  the  surface  to 
1 .5  near  the  surface.  This  is  a  nontrivial  near-field  phenom¬ 
enon  .  In  the  dipole  approximation  the  nonlinear  polarization 
P(2w)  induced  by  incident  laser  light  E(<o)  is  equal  to 

Pi(2<o)=XfjkEj(<o)Ek(a>)+---.  (1) 

Mere  xljk  is  an  electric  dipole  related  source  of  the  nonlinear 
susceptibility.  Symmetry  considerations  show  that  this  con¬ 
tribution  is  absent  in  centrosymmetric  media  and  in  this  case 
the  only  remaining  source  of  nonlinearity  is  a  surface  inter¬ 
lace  where  inversion  symmetry  is  broken.  SHG  from  a  per¬ 
fectly  flat  metal  surface  can  be  excited  mainly  with 
/'-polarized  light  while  the  much  smaller  5-polarized-light- 
induced  SHG  is  a  measure  of  the  roughness  of  the  surface.17 
In  the  near-field  region  there  is  a  substantial  contribution 
Irom  evanescent  fundamental  and  second-harmonic  fields. 
Phase-matching  condition  and  polarization  selection  rules 
are  meaningless  in  the  near-field  region  of  a  rough  surface. 
Both  are  defined  with  respect  to  the  plane  of  incidence  by  the 
direction  of  propagation  and  the  normal  to  the  surface.  There 
is  no  direction  of  propagation  for  an  evanescent  wave.  The 
normal  direction  depends  on  the  scale  of  observation  and 
does  not  coincide  with  the  average  normal  to  the  sample 
surface.  As  a  consequence  there  is  a  substantial  increase  in 
5-polarized-light-induced  evanescent  SH  field  in  the  vicinity 
of  a  rough  surface. 

Simultaneously  measured  topographical  and  near-field 
second-harmonic  images  of  the  silver  film  are  shown  in  Figs. 

3  and  4.  For  these  figures  characteristic  data  acquisition  time 
was  4-6  h.  The  power  of  the  1064  nm  light  has  been  selected 
in  a  trade  off  between  signal  magnification  and  the  rate  of 
thermal  drift  of  the  image.  The  noise  in  the  topographical 
image  is  due  to  the  laser  pulses  striking  the  sample  surface. 
The  noise  in  the  SH  image  results  from  statistical  variations 
in  the  number  of  photons  counted  in  individual  pixels. 

Near-field  images  in  Figs.  3(b)  and  3(d)  have  been  mea¬ 
sured  using  5-polarized  excitation  light.  It  appears  that  the 
main  contrast  mechanism  for  the  5-polarized-light-induced 
near-field  SH  image  is  topography  variation.  The  surface  to¬ 
pography  in  Fig.  3(a)  can  be  described  as  essentially  flat 
regions  at  the  top  and  at  the  bottom  of  the  image  separated 
by  a  narrow  groove.  The  SH  image  [Fig.  3(b)]  looks  almost 
like  a  “negative”  of  the  topographical  image.  The  SH  signal 
measured  in  the  flat  regions  is  much  smaller  than  the  signal 
measured  near  the  groove.  In  Fig.  3(c)  one  can  see  generally 
flat  regions  at  the  left  and  at  the  right  sides  of  the  image 
divided  by  a  step.  Again,  the  SH  signal  measured  near  the 
step  is  much  larger  than  the  signal  measured  in  the  fiat  re¬ 
gions  [Fig.  3(d)].  In  both  cases  the  enhancement  of  SH  gen¬ 
eration  induced  by  5-polarized  excitation  occurs  near  the 
places  on  the  surface  where  inhomogeneities  are  present. 
This  effect  has  the  same  origin  as  the  growth  of  the  ratio  of 
5-polarized-  to  P-polarized-light-induced  SH  signals  near 
the  rough  surface  discussed  above. 

Generally,  P-polarized-light-induced  SH  images  are 
much  more  complicated.  The  topographical  image  in  Fig. 
4(a)  shows  a  rather  usual  picture  of  different  crystalline 
grains  randomly  distributed  over  the  silver  film  surface.  In 
the  SH  image  these  grains  have  different  brightness:  the 
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FIG.  3.  Topography  (a),(c)  and  SH  light  distribution  (b),(d) 
measured  for  5-polarized  excitation  light. 

grain  in  the  bottom  right  comer  looks  much  darker  than  the 
grain  in  the  top  right  comer.  This  contrast  may  be  due  to  the 
different  crystalline  orientation  of  different  grains.  SH  gen¬ 
eration  is  known  to  be  sensitive  to  the  atomic  structure  of  the 
surface,2  since  in  an  isotropic  metal  such  as  silver  the  SH 
light  is  generated  within  the  top  few  atomic  surface  layers.  In 
addition,  our  experiments  have  been  conducted  under  ambi¬ 
ent  conditions.  Although  we  have  used  freshly  prepared  sil¬ 
ver  films,  surface  oxidization  could  have  affected  our  results. 
The  rate  of  surface  oxidization  depends  on  crystalline  orien¬ 
tation.  Thus,  we  may  see  a  contrast  that  is  related  to  surface 
chemistry  of  the  silver  film.  It  will  be  necessary  to  conduct 
further  experiments  in  ultrahigh  vacuum  to  clarify  this  issue. 
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FIG.  4.  Topography  (a)  and  SH  light  distribution  (b)  measured 
for  P -polarized  excitation  light.  The  cross  section  of  the  SH  image 
along  the  line  determined  by  the  arrows  is  shown  in  the  inset. 

An  interesting  feature  of  Fig.  4(b)  is  the  presence  of  small 
bright  spots  in  the  top  part  of  the  image.  The  cross  section 
through  two  of  these  spots  is  shown  in  the  inset.  The  width 
of  the  spot  shown  by  the  markers  is  equal  to  240  nm,  or 
about  half  of  the  wavelength  of  the  SH  light.  Bright  spots 
with  similar  sizes  have  been  detected  in  the  surface  near-field 
distribution  with  surface  plasmon-polariton  (SPP) 
excitation.10  These  spots  have  been  identified  as  the  result  of 
optical  field  localization. 


In  fact,  SH  generation  at  a  rough  surface  and  SPP  excita¬ 
tion  are  closely  related  phenomena  since  any  surface  defect 
is  a  source  of  a  surface  wave.15  Localization  effects  in  the 
SH  generation  from  rough  metal  surfaces  have  been  pre¬ 
dicted  by  McGum,  Leskova,  and  Agranovich.7  It  has  been 
shown  that  multiple  scattering  leading  to  light  localization 
gives  a  contribution  to  the  SH  generation.  Both  fundamental 
and  second-harmonic  wave  localization  contributes  to  the  ef¬ 
fect  giving  rise  to  components  of  SH  light  propagating  per¬ 
pendicularly  to  the  average  sample  surface  and  in  the  reverse 
direction  to  the  fundamental  wave  propagation,  respectively. 
The  latter  process  results  in  a  local  surface  enhancement  of 
the  SH  field.  The  spots  observed  in  Fig.  4(b)  may  correspond 
to  this  localization.  Previously,  the  SH  field  localization  has 
been  observed  indirectly  by  measuring  weak  changes  in  the 
far-field  angular  distribution  of  the  SH  lights.8 

IV.  CONCLUSION 

In  summary,  we  have  presented  images  of  SH  light  emis¬ 
sion  taken  in  near-field  proximity  to  a  metal  surface  together 
with  corresponding  images  of  surface  topography.  Differ¬ 
ences  in  the  mechanisms  of  SH  generation  for  different  po¬ 
larizations  of  the  excitation  light  have  been  demonstrated. 
For  S -polarized  excitation  the  observed  local  SH  generation 
is  related  mainly  to  the  topology  of  the  surface.  For 
P -polarized  light  excitation,  local  enhancement  of  SH  emis¬ 
sion  have  been  observed  that  may  be  attributed  to  SH  field 
localization.  We  believe  that  these  results  are  important  for 
further  development  of  SH  near-field  microscopy,  which 
may  become  a  useful  tool  in  magnetic  and  electric  domain 
imaging  and  has  the  potential  for  accessing  the  microscopic 
electrodynamics  of  rough  surfaces. 

ACKNOWLEDGMENTS 

We  would  like  to  acknowledge  helpful  discussions  with 
Chi  Lee  and  Julius  Goldhar.  This  research  was  supported  by 
the  Department  of  Defense  through  Contracts  Nos.  DAAH 
04-96-10279  and  MDA  904-94-C-6203. 


]K.  A.  O’Donnel,  in  Current  Trends  in  Optics,  edited  by  A.  Con¬ 
sorting  ICO  Series  (Academic  Press,  London,  1995). 

2J.  F.  McGilp,  Prog.  Surf.  Sci.  49,  1  (1995). 

3N.  Blombergen,  R.  K.  Chang,  S.  S.  Jha,  and  C.  H.  Lee,  Phys.  Rev. 
174,  813  (1968). 

4Y.  R.  Shen,  J.  Vac.  Sci.  Technol.  B  3,  1464  (1985). 

5T.  F.  Heinz,  M.  M.  T.  Loy,  and  W.  A.  Thompson,  Phys.  Rev. 
Lett.  54,  63  (1985). 

6 Nonlinear  Surface  Electromagnetic  Phenomena,  edited  by  H.-E. 
Ponath  and  G.  I.  Stegeman  (North-Holland,  Amsterdam,  1991). 

7  A.  R.  McGum,  T.  A.  Leskova,  and  V.  M.  Agranovich,  Phys.  Rev. 

B  44,  11  441  (1991). 

8  Y.  Wang  and  H.  J.  Simon,  Phys.  Rev.  B  47,  13  695  (1993). 

9D.  P.  Tsai  et  al ,  Phys.  Rev.  Lett.  72,  414  (1994). 

,0S.  I.  Bozhevolnyi,  I.  Smolyaninov,  and  A.  Zayats,  Phys.  Rev.  B 


51,  17  916  (1995);  S.  I.  Bozhevolnyi,  B.  Vohnsen,  1. 1.  Smolya¬ 
ninov,  and  A.  V.  Zayats,  Opt.  Commun.  117,  417  (1995). 

1 1  Near  Field  Optics,  edited  by  D.  W.  Pohl  and  D.  Couijon  (Kluwer, 
Dordrecht,  1993). 

12X.  Zhao  and  R.  Kopelman,  Ultramicroscopy  61,  69  (1995). 

13J.  M.  Vigoureux,  C.  Girard,  and  F.  Depasse,  J.  Mod.  Opt.  41,  49 
(1994). 

14 S.  I.  Bozhevolnyi,  O.  Keller,  and  I.  Smolyaninov,  Opt.  Lett.  19, 
1603  (1994). 

15 1.  I.  Smolyaninov,  D.  L.  Mazzoni,  and  C.  C.  Davis,  Phys.  Rev. 
Lett.  77,  3877  (1996). 

16 1. 1.  Smolyaninov,  A.  V.  Zayats,  and  C.  C.  Davis,  Opt.  Lett,  (to  be 
published). 

17  K.  J.  Song,  D.  Heskett,  H.  L.  Dai,  A.  Liebsch,  and  E.  W.  Plum¬ 
mer,  Phys.  Rev.  Lett.  61,  1380  (1988). 


1592  OPTICS  LETTERS  /  Vol.  22,  No.  21  /  November  1,  1997 


Near-field  second-harmonic  imaging  of  ferromagnetic 
and  ferroelectric  materials 


Igor  I.  Smolyaninov,  Anatoly  V.  Zayats,  and  Christopher  C.  Davis 

Department  of  Electrical  Engineering,  University  of  Maryland,  College  Park,  Maryland  20740 

Received  June  6,  1997 

A  novel  scanning  probe  technique  that  is  able  to  image  surface  magnetic  and  electric  properties  has  been 
developed.  It  is  based  on  the  near-field  microscopy  of  surface  second-harmonic  generation.  We  have 
demonstrated  the  capability  of  the  technique  by  imaging  the  domain  structure  of  a  Ni  single  crystal  and  a 
piezoelectric  ceramic.  ©  1997  Optical  Society  of  America 


Ferromagnetic  and  ferroelectric  materials  have  nu¬ 
merous  applications  in  modern  technology.  The  per¬ 
formance  of  these  materials  is  determined  by  their 
properties  at  the  microscopic  level,  particularly  by  the 
sizes  of  magnetic  (or  electric)  domains  and  domain 
boundaries.  For  example,  one  bit  of  information  in  a 
modern  magnetic  disk  is  recorded  in  a  1-^tm  size  do¬ 
main.  Ferroelectric  materials  are  widely  used  in  com¬ 
puter  DRAlM  chips.  Since  there  is  a  need  to  increase 
data  density  in  magnetic  and  electric  storage  devices, 
it  is  important  to  develop  high-resolution  tools  to  char¬ 
acterize  their  structure. 

There  are  three  major  techniques  for  magnetic  struc¬ 
ture  observations;  they  rely  on  the  use  of  polarized  elec¬ 
trons,  the  magneto-optical  Kerr  effect,  and  magnetic 
force  microscopy.1"3  The  recent  development  of  near¬ 
field  microscopy  allows  the  diffraction-limited  reso¬ 
lution  of  conventional  microscopy  to  be  overcome.4 
This  achievement  allows  one  to  use  many  potential 
advantages  of  optical  methods  in  high-resolution  mag¬ 
netic  material  imaging,  such  as  time  resolution,  spec¬ 
tral  sensitivity,  and  noninvasiveness  with  respect  to 
magnetic  structure.  Measurements  in  high  magnetic 
fields  can  be  performed  since  the  probe  used  is  non¬ 
magnetic. 

Recently,  magneto-optical  Kerr  near-field  transmis¬ 
sion  microscopy  has  been  developed  for  transparent 
samples.5  Implementation  of  a  similar  technique  in 
reflection  (to  study  opaque  samples)  is  rather  diffi¬ 
cult  because  reflected  light  often  becomes  depolarized 
when  the  near-field  microscopy  tip  is  in  close  proximity 
to  a  conducting  surface.6  Very  recently  such  measure¬ 
ments  have  been  conducted,7  although  the  achieved 
signal-to-noise  ratio  was  1:1. 

Magnetization-induced  second-harmonic  generation 
(SHG)  is  a  nonlinear  version  of  the  magneto-optical 
Kerr  technique8,9  that  has  been  shown  to  have  a  very 
large  magneto-optical  response.10  SHG  has  been  used 
successfully  in  far-field  imaging  of  magnetic-domain 
structures.11  We  report  implementation  of  this  tech¬ 
nique  in  near-field  optical  imaging  of  magnetic  materi¬ 
als  with  spatial  resolution  of  better  than  150  nm.  This 
implementation  is  based  on  recently  developed  near¬ 
field  SHG  microscopy.12 

Our  experimental  setup  for  near-field  microscopy 
of  SHG  shown  in  Fig.  1.  The  local  second-harmonic 
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(SH)  field  distribution  was  probed  with  an  uncoated 
adiabatically  tapered  fiber  tip,  which  was  drawn  at 
the  end  of  a  single-mode  fiber  by  heating  with  a  C02 
laser  in  a  micropipette  puller.  The  fiber  tip  can  be 
scanned  over  the  sample  surface  with  a  constant  tip- 
surface  distance  of  a  few  nanometers  by  shear-force 
feedback  control.  Therefore,  surface  topography  can 
be  imaged  with  resolution  on  a  nanometer  scale,  while 
the  SH  near-field  image  is  simultaneously  recorded. 
SHG  was  excited  at  an  angle  of  incidence  of  60° 
with  a  Nd:YAG  laser  operating  at  1064  nm  (repetition 
rate,  10  Hz;  pulse  duration,  20  ns;  pulse  energy,  10- 
15  mJ).  The  SH  signal  was  measured  with  a  photo¬ 
multiplier  and  gated  electronics.  The  signal  depends 
quadratically  on  the  fundamental  light  power.  The 
SH  signal  at  every  point  of  the  image  was  averaged 
over  a  few  tens  of  laser  pulses,  so  image-acquisition 
time  was  several  hours.  This  drawback  can  be  over¬ 
come  by  use  of  a  Ti:sapphire  laser  with  a  high  repeti¬ 
tion  rate. 

In  the  dipole  approximation  the  nonlinear  polariza¬ 
tion  P(2co)  of  a  magnetic  medium  induced  by  incident 


magnetic  sample 

Fig.  1.  Experimental  setup  for  near-field  SH  microscopy. 
An  auxiliary  laser  is  used  for  shear-force  feedback  control. 
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laser  light  E(w)  is  equal  to 

Pit 2«)  =  xfjkEj(co)  +  xfaEjMEMM,  +  ....  (1) 

Here  xtjk  an  electric  dipole  contribution  to  the 
nonlinear  susceptibility,  xTjki is  magnetization-induced 
term,  and  M  is  the  magnetization  of  the  medium  (the 
second  term  describes  direct  magnetization-induced 
SH  polarization  and  the  influence  of  the  magnetization 
on  the  symmetry  of  the  electric  dipole  contribution). 
Symmetry  considerations  show  that  both  terms  are  ab¬ 
sent  in  centrosymmetric  media,  and  the  only  remain¬ 
ing  source  of  nonlinearity  is  a  surface  interface  where 
inversion  symmetry  is  broken.  SHG  arises  from  the 
top  few  atomic  surface  layers.  Magnetization-induced 
contrast  depends  on  the  polarization  state  of  the  ex¬ 
citation  light  and  the  orientation  of  the  magnetiza¬ 
tion  vector  in  the  sample.  At  zero  magnetization,  SHG 
from  a  perfectly  flat  metal  surface  can  be  excited 
mainly  with  p-polarized  light,  while  the  much  smaller 
s-polarized-light-induced  SHG  is  generally  a  measure 
of  the  roughness  of  the  surface13  (since  metals  are  not 
exactly  described  by  a  jellium  model  there  is  also  some 
contribution  that  is  due  to  the  crystal  structure  of 
metals).  Under  p-polarized  excitation,  the  motion  of 
electrons  has  a  component  perpendicular  to  the  sur¬ 
face.  This  perpendicular  motion  of  electrons  in  the  top 
atomic  layers  is  strongly  affected  by  the  presence  of 
the  surface.  The  consequent  asymmetric  motion  gives 
rise  to  SH  generation.  Under  s-polarized  excitation, 
the  driving  field  is  parallel  to  the  surface,  and  SHG  is 
much  weaker  in  this  case. 

If  the  magnetization  is  nonzero  and  has  a  component 
parallel  to  the  surface  and  perpendicular  to  the  electric 
field  of  the  s-polarized  excitation  light  (Fig.  1),  the  mo¬ 
tion  of  electrons  near  the  surface  becomes  asymmetric 
because  of  the  Lorentz  force  acting  perpendicular  to  the 
surface.  This  gives  rise  to  SH  generation.  This  SHG 
is  insensitive  to  the  sign  of  magnetization.  The  com¬ 
ponent  of  the  electric  field  of  the  p-polarized  excitation 
light  parallel  to  the  surface  gives  similar  SHG.  In  this 
case  magnetization-induced  SHG  adds  to  the  nonmag¬ 
netic  contribution  from  the  component  of  the  electric 
field  perpendicular  to  the  surface.  Equation  (1)  shows 
that  p-polarized-light-induced  SHG  is  sensitive  to  the 
sign  of  magnetization,  since  nonlinear  polarization  is 
linear  in  the  magnetization  vector. 

In  a  jellium  model,  magnetization  perpendicular  to 
an  ideally  flat  surface  does  not  contribute  to  the  SHG 
for  both  polarizations  of  the  excitation  light.  However, 
this  magnetization  can  be  allowed  by  the  symmetry 
of  the  surface  electronic  band  structure  and  has  been 
observed  recently.14  It  can  also  be  observed  from  a 
rough  surface  since  the  normal  direction  to  the  surface 
is  not  well  defined.12  We  observed  this  behavior  in  our 
experiment  with  a  strong  permanent  SmCo  magnet. 
The  SH  intensity  observed  from  the  side  surface  of 
the  magnet  was  almost  an  order  of  magnitude  stronger 
than  from  the  pole  surface. 

Taking  into  account  the  considerations  above,  we 
chose  to  use  s-polarized  excitation  light  in  magnetic- 
domain  imaging.  We  studied  magnetic  domains  at  the 
[100]  surface  of  a  Ni  single  crystal.  Before  the  ex¬ 
periment  the  sample  had  not  been  specially  magne¬ 


tized.  Images  of  surface  topography  and  spatial  dis¬ 
tribution  of  SHG  are  shown  in  Figs.  2(a)  and  2(b). 
Bright  and  dark  regions  with  characteristic  rectan¬ 
gular  shapes  are  seen  in  the  SH  image.  These  fea¬ 
tures  do  not  correlate  with  surface  topography,  which 
is  rather  flat.  The  rectangular  features  resemble  the 
characteristic  shapes  of  the  domains  of  closure  in  fer¬ 
romagnetics.15  We  conclude  that  the  regions  with  uni¬ 
form  brightness  are  magnetic  domains  with  different 
orientations  of  magnetization.  The  sizes  of  the  do¬ 
mains  observed  are  of  the  order  of  a  few  micrometers, 
which  is  the  usual  size  of  magnetic  domains  in  Ni.  Af¬ 
ter  the  sample  is  exposed  to  a  strong  magnetic  field 
perpendicular  to  its  surface,  the  domain  structure  is 
changed,  as  can  be  seen  in  the  SH  image  in  Fig.  2(c). 
SH  light  intensity  decreased  since  the  induced  magne¬ 
tization  was  perpendicular  to  the  surface. 

The  cross  section  of  the  SH  image  in  Fig.  2(b) 
perpendicular  to  a  domain  wall  is  shown  in  the  inset  of 
Fig.  2.  The  width  of  a  transition  region  between  a  pair 
of  domains  is  equal  to  150  nm.  This  number  is  close  to 
the  known  thickness  of  a  Bloch  wall  (magnetic-domain 
boundary),15  that  is,  —100  nm,  and  gives  us  an  estimate 
of  the  resolution  of  our  near-field  SH  microscope. 

A  ferroelectric  material  that  we  studied  is  poled  lead 
zirconate  titanate  piezoceramic.  Crucial  parameters 
of  piezoceramic  performance  are  hysteresis,  linearity, 
and  creep,  which  depend  on  the  quality  of  poling. 
A  typical  image  of  the  SH  intensity  distribution  is 
shown  in  Fig.  3.  It  was  recorded  with  p-polarized 
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Fig.  2.  (a)  Topography  and  (b),  (c)  SH  light  distribution 
measured  for  the  Ni  single  crystal  (b)  before  and  (c)  after 
exposure  to  a  strong  magnetic  field.  The  inset  shows  the 
cross  section  of  the  SH  image  in  (b)  perpendicular  to  a 
domain  wall. 
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Fig.  3.  SH  intensity  distribution  over  the  piezoceramic 
surface.  This  image  is  a  superposition  of  the  surface 
topography  shown  as  a  quasi-three-dimensional  view  and 
the  half-tone  SH  image. 

fundamental  light.  The  poling  vector  was  parallel 
to  the  sample  surface.  The  polycrystalline  structure 
of  the  ceramic  is  clearly  seen  in  the  topography. 
Different  crystallites  show  different  brightness  of  SH 
emission  that  may  correspond  to  variations  in  poling. 
One  of  the  crystallite  boundaries  (shown  by  the  arrows) 
also  shows  a  difference  in  SH  brightness  from  the 
rest  of  the  material.  We  believe  that  such  defects 
contribute  to  the  effects  of  hysteresis  and  creep  of 
piezoceramics. 

In  summary,  we  have  developed  a  novel  scanning 
probe  technique  that  is  able  to  image  surface  magnetic 
and  electric  properties.  It  is  based  on  near-field  op¬ 
tical  imaging  of  the  surface.  We  have  demonstrated 


resolution  of  better  than  150  nm  by  using  uncoated  ta¬ 
pered  fiber  tips. 
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We  present  a  focused  ion-beam  (FIB)  fabrication  method  for  very  clean  and  well-defined 
subwavelength  fiber  probes  with  metallic  apertures  of  a  desired  diameter  for  use  in  near-field 
scanning  optical  microscopy.  Such  probes  exhibit  improved  features  compared  to  probes  coated 
with  metal  by  the  conventional  angled  evaporation  technique.  Examples  of  FIB  fabricated  fiber 
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In  the  past  decade,  near-field  scanning  optical  micros¬ 
copy  (NSOM)  has  successfully  demonstrated  its  capability  of 
subwavelength  optical  resolution  by  breaking  the  diffraction 
limit  of  conventional  optical  microscopy.1’2  The  optical  near 
field  of  a  source  exists  within  a  distance  of  less  than  one 
optical  wavelength  from  it.  By  following  Synge’s3  original 
suggestion,  involving  the  use  of  a  tiny  aperture  (diameter 
<^\)  in  a  metal  screen  that  is  placed  in  near-field  proximity 
to  a  sample,  optical  microscopy  can  achieve  a  resolving 
power  below  the  classical  Abbe  limit.  Practical  implementa¬ 
tion  of  a  tiny  aperture  in  a  conductive  screen  is  accomplished 
by  tapering  down  an  optical  fiber  to  subwavelength  diameter 
using  the  heating-pulling  technique4  followed  by  metallic 
overcoating  of  the  sides  of  the  fiber  to  define  an  aperture. 

The  major  drawback  of  pulled  near-field  tips  is  their  low 
optical  throughput,  which  is  a  key  parameter  influencing 
their  better  utilization  in  applications  such  as  optical 
spectroscopy,5  nanofabrication,6  and  high-density  data 
storage.7  To  overcome  this  problem  many  researchers8- 1 1 
have  attempted  to  fabricate  NSOM  probes  of  a  desired  shape 
and  with  higher  throughput.  For  example,  Zeisel  et  ai9  used 
a  chemical  etching  technique  similar  to  that  developed  ear¬ 
lier  by  Pangaribuan  et  al .8  for  a  photon  scanning  tunneling 
microscopy.  The  resultant  fiber  tips  had  a  much  shorter  taper 
region  and  higher  cone  angle;  they  provided  a  three  orders  of 
magnitude  higher  throughput  than  normally  obtained  with 
pulled  fiber  probes.  In  another  experiment  Islam  and 
co-workers11  fabricated  hybrid  pull-etch  tapered  probes  with 
throughput  four  orders  of  magnitude  higher  than  that  of 
pulled  fibers. 

The  process  of  applying  a  metallic  overcoat  to  the  tip  is 
an  equally  important  aspect  of  fabricating  high-resolution 
NSOM  probes.4  Through  the  addition  of  an  opaque  metallic 
coating  such  as  aluminum  that  covers  the  sidewalls  of  the 
tapered  fiber,  a  tiny  aperture  at  the  end  of  the  fiber  is  created. 
The  usual  process  for  coating  the  tip  is  through  vacuum 
deposition  of  aluminum  by  thermal  evaporation  while  the  tip 
is  rotated  and  held  at  an  angle  to  the  source.  The  exposed 
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glass  fiber  end  resulting  in  this  manner  is  generally  located 
inside  the  metallic  aperture,  is  partially  obstructed  by  rough 
aluminum  grains  situated  near  the  aperture  boundary,  and 
often  contains  aluminum  protrusions,  as  evidenced  by  nu¬ 
merous  scanning  electron  microscope  (SEM)  images  of  fiber 
tips.12,13  This  may  contribute  to  the  poor  polarization  capa¬ 
bilities  of  the  majority  of  fiber  tips  prepared  by  angled 
evaporation.14  A  larger  than  optimal  separation  of  fiber  tip 
and  sample  and  the  absence  of  a  clean  opening  result  in 
lower  throughput  of  the  probe  and  poor  optical  resolution  of 
the  NSOM.  Furthermore,  the  angle-evaporation  procedure  is 
very  time  consuming  and  may  not  result  in  the  fabrication  of 
a  batch  of  reproducible  and  well-defined  apertures  of  desired 
diameter.  In  addition,  other  problems  associated  with  angle 
evaporation  are  the  fabrication  of  a  special  holder  to  accom¬ 
modate  many  fibers  with  adjustable  tilt  angle,  the  loading  of 
the  probes,  rotation  of  the  tips  about  their  axis  during  the 
coating  operation,  and  keeping  the  fibers  untangled  through¬ 
out  the  coating  procedure. 

It  is  the  aim  of  this  work  to  develop  a  method  for  fabri¬ 
cating  clean,  well-defined,  and  highly  reproducible  subwave¬ 
length  apertures  for  NSOM  applications  using  a  focused  ion- 
beam  (FIB)  milling  technique.  Fiber  probes  were  produced 
through  the  use  of  a  heating-pulling  process  and  by  adjust¬ 
ing  the  parameters  described  by  Valaskovic  et  al.4  for  ob¬ 
taining  probes  of  optimum  shapes  and  throughputs.  The 
pulled  fiber  probes  were  then  uniformly  coated  all  over  with 
a  100-150  nm  aluminum  film  in  batches  of  10-35  fibers 


FIG.  1.  Schematic  diagram  of  the  FIB  milling  procedure  to  produce  a  near¬ 
field  aperture.  The  rectangular  shaded  region  represents  the  area  in  which 
the  ion-beam  raster  scans  over  the  tip  from  right  to  left  in  order  to  remove 
material  cleanly. 


0003-695 1/98/72(24)/31 33/3/$  1 5.00 


3133 


©  1998  American  Institute  of  Physics 


3134 


Appl.  Phys.  Lett.,  Vol.  72,  No,  24,  15  Juno  1998 


Pilevar  et  al. 


(C) 

FIG.  2.  SEM  photographs  of  the  near-field  optical  fiber  probe  (a)  before  FIB 
milling  (bar=2/xm),  (b)  after  FIB  fabrication  of  a  100  nm  diam  aperture 
(bar=l  /am),  and  (c)  after  FIB  fabrication  of  a  300  nm  diam  aperture 
(bar=I  /im). 


with  a  sputtering  unit  (Edwards  Xenosput  XE200)  usually 
used  to  metal  coat  specimens  prior  to  inspection  by  scanning 
electron  microscopy.  The  coating  is  smooth  and  exhibits  no 
leakage  or  pin  holes  upon  injection  of  laser  light  into  such 
fiber  probes. 

In  order  to  use  the  tips  in  a  near-field  microscope,  a 
subwavelength  opening  at  the  apex  of  the  probe  must  be 
formed.  We  have  used  FIB  milling  to  fabricate  an  aperture  of 
desired  diameter  at  the  very  end  of  the  metal  coated  fiber. 
The  FIB  system  used  was  a  FEI  Dualbeam  620  FIB/SEM 
workstation  (30  keV  30  keV  Ga+  ions,  minimum  spot  size  5 
nm).  The  electron  beam  of  the  system  imaged  the  profile  and 
diameter  of  each  individual  tip  in  the  batch  containing  5-10 
fibers  positioned  inside  the  chamber.  To  remove  the  metal 
coating  and  thus  fabricate  the  aperture,  the  ion  beam  was 
raster  scanned  in  a  rectangular  pattern  at  the  apex  of  the 
fiber,  as  represented  by  the  shaded  region  in  Fig.  1.  In  con¬ 
trast  to  simply  cutting  off  the  end  of  the  fiber  in  one  sweep 


(c) 

FIG.  3.  Simultaneous  (a)  shear-force  topographical  and  (b)  near-field  optical 
images  of  a  10X  10  jxm2  region  of  a  test  sample  using  the  100  nm  aperture 
diameter  FIB  fabricated  fiber  probe.  An  expanded  view  of  a  region  (3.25 
X3.25  Atm2)  of  the  100  nm  lines  is  shown  in  (c). 


with  the  FIB,  which  would  result  in  redeposition  of  material 
onto  the  aperture,  this  procedure  completely  removes  the  ma¬ 
terial  at  the  end  of  the  fiber.15  By  varying  the  size  of  the 
rectangular  scanned  region  and  with  knowledge  of  the  coat¬ 
ing  thickness  and  fiber  end  diameter,  we  can  precisely  make 
a  clean-cut  aperture  of  the  desired  diameter  and  optical 
throughput  suitable  for  a  variety  of  NSOM  applications. 

Figure  2(a)  shows  a  SEM  image  of  a  tapered  fiber  probe 
coated  with  150  nm  of  aluminum.  The  aperture  is  then 
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formed  with  the  aforementioned  FIB  milling  procedure  cre¬ 
ating  an  extremely  well-defined  100  nm  diam  opening,  as 
shown  in  Fig.  2(h).  Figure  2(c)  shows  another  aluminized 
fiber  lip  milled  by  the  FIB  to  make  a  300  nm  diam  aperture. 
The  entire  FIB  milling  time  for  each  tip  was  between  2  and 
t  min.  As  revealed  in  Fig.  2.  the  fabricated  near- field  aper¬ 
tures  are,  to  the  authors'  best  knowledge,  the  besl-delined 
and  clean  cut  ones  ever  fabricated.  In  contrast  to  angled 
metal  evaporation,  the  FIB  milling  procedure  removes  any 
rough  aluminum  grains  at  the  tip  end  .surrounding,  the  aper¬ 
ture.  The  advantages  of  this  are  twofold:  the  optical  through 
put  is  improved  and  the  aperture  can  give  the  true  near-field 
optical  contrast  of  a  sample  without  the  influence  of  tip  pro¬ 
trusions. 

The  limitations  of  our  FIB  tip  aperture  fabrication  tech¬ 
nique  involve  the  coating  uniformity,  tip  geometry,  tip 
charging  under  ion  imaging,  and  the  ion-beam  profile.  An  rf 
sputterer  was  used  to  achieve  a  smooth  and  uniform  coating 
over  the  tip  so  that  for  a  known  coating  thickness,  one  can 
cut  the  aperture  to  a  desired  diameter.  Clearly,  nonuniformity 
of  the  coating  can  lead  to  errors  in  the  desired  aperture  di¬ 
ameter.  Tip  geometry  is  also  a  factor  since  a  high  aspect  ratio 
tip  will  allow  for  more  precise  control  of  the  final  aperture 
diameter  (this  is  especially  important  for  fabrication  of  aper¬ 
tures  smaller  than  100  nm).  A  tradeoff  must  be  struck,  how¬ 
ever,  because  high  aspect  ratio  tips  have  less  throughput  be¬ 
cause  the  evanescent  wave  propagates  through  a  longer 
region.  Our  FIB  has  a  Gaussian  profile  with  a  minimum  5 
nm  spot  size,  allowing  for  precise  aperture  milling  of  even 
ultrahigh  aspect  ratio  tips.  A  final  important  consideration  is 
the  grounding  of  the  tip  to  avoid  charging  while  imaging. 
This  was  accomplished  using  a  small  amount  of  silver  paint, 
which  was  then  easily  removed  with  acetone. 

FIB  fabricated  fiber  probes  were  tested  in  the  near-field 
system  developed  in  our  laboratory.  Near-field  distance  regu¬ 
lation  is  maintained  through  the  use  of  shear  force  phase- 
locked  tuning  fork  feedback.16  The  sample  imaged  was  a 
series  of  lines  of  various  widths  (1  yum-100  nm)  milled  by 
the  FIB  through  a  50  nm  aluminum  film  coated  on  a  glass 
substrate.  The  sample  was  imaged  using  a  FIB  milled  fiber 
tip  of  100  nm  aperture  diameter  in  our  NSOM  setup  operat¬ 
ing  in  collection  mode  using  a  He-Ne  laser  source  (\=633 
nm).  In  Figs.  3(a)  and  3(b),  we  show  the  topographical  and 
near-field  optical  images  of  a  lOXlOyttm2  region  of  the 
sample,  respectively.  The  smallest  line  imaged  at  the  left¬ 
most  area  is  approximately  100  nm  wide.  A  magnified  near¬ 
field  image  of  the  sample  showing  optically  resolved  lines  of 
less  than  200  nm  is  shown  in  Fig.  3(c).  Because  the  probe 
end  is  flat  and  free  of  aluminum  protrusions,  the  optical  con¬ 
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trast  is  much  less  affected  by  topographically  induced  inten¬ 
sity  changes.17  Thus,  the  optical  images  in  Figs.  3(b)  and 
3(c)  yield  optical  resolution  of  about  100  nm  (when  100  nm 
width  of  the  test  line  itself  is  taken  into  account),  as  expected 
for  the  fabricated  100  nm  aperture.  The  signal-to-noise  ratio 
in  the  optical  image  is  rather  low  because  the  sample  itself 
defines  an  aperture  smaller  than  the  wavelength  with  a  thick¬ 
ness  of  50  nm.  Thus,  the  optical  image  is  taken  in  a  situation 
that  is  equivalent  to  a  double  pass  through  the  NSOM  tip. 

In  conclusion,  we  have  implemented  a  FIB  technique  to 
fabricate  well-defined  near-field  apertures  of  desired  diam¬ 
eters  that  exhibit  improved  features  over  the  conventional 
angled  metal  evaporation  method.  Currently  under  investiga¬ 
tion  is  the  use  of  FIB  milling  applied  to  chemically  etched 
liber  probes,  which  has  the  potential  for  creating  a  probe 
with  a  much  higher  optical  throughput  and  smaller  aperture. 
This  combined  technique  could  boost  the  optical  signal  by 
orders  of  magnitude  and  improve  NSOM  resolution  signifi¬ 
cantly. 

This  research  was  supported  by  (he  Department  of  De¬ 
fense  through  Contract  No.  MDA  904-94-C-6203.  The  au¬ 
thors  would  like  to  thank  Dr.  John  Orloff  for  helpful  sugges¬ 
tions  and  Victor  E.  Yun  for  help  with  the  sputter  coating. 


'D.  W.  Pohl,  W.  Denk.  ami  M.  lain/..  Appl.  Phys.  Lett.  44.  651  (1984). 
2E.  Betzig,  J.  K.  Trautman.  T.  D.  Harris.  J.  S.  Weiner,  ami  R.  I..  Kosielak, 
Science  251,  1468  (1992). 

3E.  H.  Synge,  Philos.  Mag.  6.  356  (1928). 

4G.  A.  Valaskovic,  M.  Holton,  and  G.  H.  Morrison.  Appl.  Opt.  34.  125 
(1995). 

5W.  P.  Ambrose,  P.  M.  Goodwin,  J.  C.  Marlin,  and  R,  A.  Keller,  Science 
265,  364  (1994). 

6 A.  Shchemelinin,  M.  Rudman,  K.  Liberman,  and  A.  Lewis.  Rev.  Sci. 
Instrum.  64,  3538  (1993). 

7E.  Betzig,  J.  K.  Trautman,  R.  Wolfe,  E.  M.  Gyorgy.  and  P.  L.  Finn.  Appl. 
Phys.  Lett.  61,  142  (1992). 

8T.  Pangaribuan,  K.  Yamada,  S.  Jiang.  H.  Ohsawa.  and  M.  Ohisu.  Jpn.  J. 
Appl.  Phys.,  Part  1  31,  1302  (1992)." 

9D.  Zeisel,  S.  Nettesheim,  B.  Dutoit,  and  R.  Zenobi.  Appl.  Phys.  Lett.  68, 
2491  (1996). 

10R.  U.  Maheswari,  S.  Mononobe,  and  M.  Ohtsu,  Appl.  Opt.  35,  6740 
(1996). 

11 M.  N.  Islam,  X.  K.  Zhao,  A.  A.  Said,  and  C.  F.  Vail.  Conference  on  Lasers 
and  Electro-Optics  (Optical  Society  of  America,  Washington,  DC.  1997). 
Vol.  11,  p.  70. 

12 E.  Betzig  and  J.  K.  Trautman,  Science  257,  189  (1992). 

,3C.  Obermuller  and  K.  Karrai,  Appl.  Phys.  Lett.  67,  3408  (1995). 

I4C.  Durkan,  I.  V.  Shvets,  and  J.  C.  Lodder,  Appl.  Phys.  Lett.  70,  1323 
(1997). 

15  K.  P.  Muller  and  H.  C.  Petzold,  Proc.  SPIE  1263,  12  (1990). 

I6W.  A.  Alia  and  C.  C.  Davis,  Appl.  Phys.  Lett.  70,  405  (1997). 

17 B.  Hecht,  H.  Bielefeldt,  Y.  Inouye,  D.  W.  Pohl,  and  L.  Novotny,  J.  Appl. 
Phys.  81,  2492  (1997). 


1346  OPTICS  LETTERS  /  Vol.  23,  No.  17  /  September  1, 1998 


Apparent  superresolution  in  near-field  optical  imaging 

of  periodic  gratings 


Igor  I.  Smolyaninov  and  Christopher  C.  Davis 

Department  of  Electrical  Engineering,  University  of  Maryland,  College  Park ,  Maryland  20742 

Received  June  10,  1998 

A  superresolution  phenomenon  in  images  of  diffraction  gratings  obtained  with  a  near-field  optical  microscope 
probe  positioned  at  distances  of  the  order  of  1  fim  from  the  surface  was  recently  reported  [Appl.  Opt.  33,  876 
(1994)].  This  experiment  was  interpreted  in  terms  of  evanescent  waves  observed  in  the  far  field.  We  show 
that  the  effect  observed  in  that  experiment  is  a  very  well  known  far-field  optical  phenomenon  called  the  Talbot 
effect.  The  existence  of  self-imaging  Talbot  planes  for  gratings  with  periods  of  the  order  of  a  few  micrometers 
appears  promising  for  many  practical  applications.  ©  1998  Optical  Society  of  America 
OCIS  codes:  180.5810,  070.6760. 


Near-field  scanning  optical  microscopy  (NSOM)  has  be¬ 
come  an  important  tool  in  surface  studies.1  However, 
in  many  situations  the  problems  of  the  contrast  mecha¬ 
nism  and  resolution  of  NSOM  still  require  discussion. 
A  recent  overview  of  these  questions  is  available.2 

An  interesting  NSOM  experiment  was  reported3  that 
has  given  rise  to  ongoing  discussion  and  controversy. 
A  study  of  the  imaging  properties  of  the  external- 
reflection  NSOM  setup  was  conducted  with  differ¬ 
ent  diffraction  gratings  as  samples.  Surprisingly,  the 
best  profiles  of  a  1.6-/xm-period  rectangular  diffrac¬ 
tion  grating  were  measured  at  a  tip -surface  distance 
of  0.7  pm.  This  experimental  observation  was  used 
to  justify  a  general  theoretical  claim  that  evanescent 
waves  contribute  to  the  radiation  fields  of  sources  and 
to  the  far  fields  of  scatterers.4  Recently  Wolf  and  Foley 
demonstrated,5  by  considering  a  spherical  scalar  wave 
and  a  linear  electric  dipole  field,  that  this  claim  is  un¬ 
true.  The  only  question  that  remains  unanswered  is: 
what  was  seen  in  the  experiment?  In  this  Letter  we 
show  that  the  superresolution  observed  in  the  experi¬ 
ment  reported  in  Ref.  3  is  a  manifestation  of  a  well- 
known  far-field  optical  phenomenon  called  the  Talbot 
effect.6 

The  Talbot  effect  exhibits  itself  when  a  grating  is  il¬ 
luminated  with  a  highly  spatially  coherent  plane  wave. 
Self-images  of  the  grating  are  produced  in  certain 
planes  behind  and  in  front  of  the  grating.  These  im¬ 
ages  are  produced  solely  by  free-space  propagation  of  a 
diffracted  optical  field.  Applications  of  the  Talbot  ef¬ 
fect  include  optical  image  processing  and  production  of 
optical  elements  (see,  for  example,  Refs.  7  and  8). 

The  theory  of  the  Talbot  effect  is  most  easily  derived 
in  the  one-dimensional  case.  The  scalar  field  behind 
an  infinite  grating  with  a  period  a  positioned  at  the 
origin  (z  =  0),  when  the  grating  is  illuminated  by  a 
plane  wave,  can  be  described  by  Fourier  series: 

E(x,z)  =  ^Cn  exp(27rinx/a)exp(iknz), 

n 

kn  =  27r(l/A2  -  n2/a2)m.  (1) 

The  usual  integer  Talbot  effect  happens  when  a  »X 
and  knz  ~  2'irz/ \  -  2irn2(\z/2a2).  Since  the  square 
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of  an  integer  n  is  itself  an  integer,  at  Zt  =  2 a2/ A  all 
Fourier  components  of  the  optical  field  have  exactly  the 
same  phase  differences  with  respect  to  one  another  as 
they  have  at  plane  z  =  0.  Thus  the  field  distribution 
of  a  grating  reproduces  itself.  The  same  self-imaging 
is  observed  at  distances  that  are  integer  multiples  of 
the  Talbot  distance:  z  —  nZt.  More-complex  field  pat¬ 
terns  occur  at  fractional  values  of  the  Talbot  distance. 
These  patterns  are  located  at  z  ~  ( p/q)Zt ,  where  p 
and  q  are  integers.  In  the  case  of  the  fractional  Tal¬ 
bot  planes  many  Fourier  components  of  the  field  (but 
not  all  of  them)  have  the  same  phase  differences  with 
respect  to  one  another,  so  the  field  pattern  resembles 
(but  does  not  exactly  reproduce)  the  field  distribution 
at  z  =  0.7 

The  experiment  discussed  in  Ref.  3  and  the  Talbot  ef¬ 
fect  observations  have  the  following  common  features: 

(a)  field  distributions  at  some  plane  near  gratings  illu¬ 
minated  by  spatially  coherent  sources  are  studied,  and 

(b)  sudden  contrast  enhancement  or  self-imaging  oc¬ 
curs  at  a  certain  distance  from  the  grating.  The  differ¬ 
ence  between  these  experiments  is  that  in  Ref.  3a>  A. 
Since  the  experimental  geometry  does  not  depend  on 
the  y  coordinate,  the  optical  field  remains  scalar,  and 
Eq.  (1)  is  still  true  in  this  case.  But  only  a  few  terms 
in  Eq.  (1)  are  now  propagating  waves,  and  the  rest  are 
evanescent.  Therefore  we  cannot  use  an  expansion  of 
knz  for  small  X/a  anymore.  On  the  other  hand,  it  is 
possible  that  at  some  modified  Talbot  distance  Ztm  the 
phase  differences  among  a  small  number  of  propagat¬ 
ing  Fourier  terms  will  be  exactly  the  same  as  in  the 
plane  2  =  0.  For  this  to  occur  the  following  condition 
must  be  satisfied: 

(l/A2  -  1/ a2)V2Ztm  =  (1/A2  -  2 2/a2)mZtm  +  n2 

=  (1/A2  —  32/  a2)mZtm  +  n.3  =  , 

(2) 

where  ri2,  etc.,  are  different  integer  numbers. 
Equation  (2)  can  have  a  number  of  solutions  for  Ztm 
that  are  no  longer  periodically  spaced.  Nevertheless, 
the  field  distribution  at  the  planes  2  =  Ztm  will  possess 
the  Talbot  effect  property  of  being  self-focused. 
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Fig.  1.  Field  distributions  of  a  grating  at  different  dis¬ 
tances  near  the  Talbot  plane  at  Ztm  =  3.36  yzm. 


In  the  following  example  we  show  that  Talbot  planes 
can  indeed  appear  near  gratings  with  periods  a  ~  A. 
To  make  the  numerical  analysis  as  simple  as  possible 
we  consider  a  monochromatic  (A  =  0.633  ^m)  plane 
wave  propagating  through  a  periodic  set  of  narrow 
openings  made  in  a  plane  opaque  screen.  To  keep  our 
considerations  simple  we  assume  that  the  period  of  the 
grating  is  slightly  larger  than  3A  and  that  the  size  of 
the  openings  is  d  <<A,  so  there  are  just  three  propa¬ 
gating  Fourier  components  of  the  optical  field  behind 
the  grating  with  approximately  the  same  amplitude 
( Ci  =  C2  =  C3).  Simple  numerical  calculations  show 
that  such  a  grating  with  a  period  a  =  1.936  /zm  ex¬ 
hibits  self-focusing  at  a  distance  Ztm  =  3.36  /zm.  In 
this  case  the  integer  numbers  in  Eq.  (2)  are  /i2  =  1 
and  ns  =  4.  Optical  field  intensities  calculated  with 
Eq.  (1)  for  different  planes  near  the  Talbot  plane  at 
z  =  Ztm  —  3.36  Jim  are  shown  in  Fig.  1.  The  planes 
are  selected  far  enough  from  the  grating  that  there 
is  no  comparable  contribution  from  the  evanescent 
Fourier  components  of  the  diffracted  field.  Figure  1 
shows  a  sudden  self-focusing  effect.  The  full  widths 
at  half-maximum  of  the  periodic  peaks  of  the  opti¬ 
cal  field  are  equal  to  0.23  /zm  =  0.36  A.  Here  we 
should  recall  that  the  theoretical  limit  to  the  reso¬ 
lution  of  a  near-field  scanning  microscope  with  un¬ 


coated  tapered  optical  fiber  tips  is  A/3.9  Thus  the 
experimental  results  obtained  in  the  study  with  un¬ 
coated  fiber  tips  reported  in  Ref.  3  no  longer  look  sur¬ 
prising.  While  imaging  gratings,  one  can  expect  to 
achieve  optical  resolution  of  the  order  of  A/3  in  some 
Talbot  planes  in  the  far-field  that  are  many  microme¬ 
ters  away  from  the  sample.  Such  results  do  not  re¬ 
quire  any  nontraditional  explanation. 

On  the  other  hand,  the  existence  of  self-imaging 
Talbot  planes  for  gratings  with  periods  of  the  order  of  a 
few  micrometers  is  itself  an  interesting  phenomenon, 
since  self-imaging  can  be  used  in  many  practical 
applications.  For  example,  self-imaging  can  be  used 
in  the  production  of  optical  elements  with  much  smaller 
features  than  those  obtained  with  the  usual  integer  or 
fractional  Talbot  effect. 

In  summary,  we  have  shown  that  experimentally 
detected  superresolution  in  near-field  scanning  mi¬ 
croscopes  working  at  micrometer  distances  from  pe¬ 
riodic  gratings  used  as  samples  is  a  manifestation  of 
the  Talbot  effect.  Superresolution  does  not  require  an 
explanation  involving  unusual  physics  such  as  eva¬ 
nescent-wave  detection  in  the  far  fields  of  scatterers. 
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Abstract 


The  physical  mechanism  of  shear-force  interaction  between  a  fiber  tip  and  a  sample  in  near-field  optical  microscopy 
has  been  studied  using  tunneling  current  measurements  between  the  tip  and  a  sample.  Optical  and  non-optical  (tuning 
fork)  shear-force  distance  controls  have  been  used.  In  both  cases  a  tunneling  current  between  the  tip  and  the  sample  has 
been  detected  at  the  very  beginning  of  shear-force  approach  curve  when  the  amplitude  of  tip  vibration  was  just  slightly 
smaller  than  its  initial  value.  These  results  indicate  that  fiber-tips  generally  mechanically  touch  samples  under  shear-force 
distance  feed-back  control.  ©  1998  Elsevier  Science  B.V.  All  rights  reserved. 

PACS:  07.79.Fc 
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merits.  Several  possible  interaction  mechanisms 
have  been  suggested  to  contribute  to  SFDC,  such 
as  capillary  and  van  der  Waals  forces  [1].  More 
recently  it  has  been  suggested  that  an  effective  non¬ 
linear  bending  force  (NBF)  connected  with  one¬ 
sided  mechanical  probe-surface  contact  is  a  main 
source  of  the  shear-force  [2].  The  latter  mechanism 
could  be  the  “universal”  one  that  is  common  to  all 
environments. 

It  is  possible  to  detect  a  mechanical  contact  be¬ 
tween  a  conductive  probe  and  a  conductive  sample 
by  measuring  the  probe-sample  electric  resistance 
[3].  For  relatively  large  tip-surface  separations  of 
several  nanometers,  the  current  increases  exponen¬ 
tially  with  decreasing  separation  (at  a  fixed  tunnel¬ 
ing  voltage).  On  further  approach,  a  discontinuous 
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I.  Introduction 

Shear-force  distance  control  (SFDC)  is  com¬ 
monly  used  in  near-field  optical  microscopy.  It  is 
based  on  a  decrease  in  the  oscillation  amplitude  of 
a  resonantly  dithered  tapered  fiber  during  its  ap¬ 
proach  to  a  surface.  The  nature  of  SFDC  has  been 
the  subject  of  ongoing  discussion,  since  it  is  used 
with  equal  success  in  air,  in  ultra  high  vacuum,  and 
in  superfluid  helium.  This  “universal”  nature  of  the 
shear  force  suggests  some  common  physical  mecha¬ 
nism  of  fiber  tip  damping  in  these  different  environ- 
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current  jump  is  observed  [4],  which  is  attributed  to 
mechanical  instability  caused  by  significant  atomic 
motion  under  the  influence  of  adhesive  forces,  in¬ 
itiating  point  contact  between  tip  and  sample  [5]. 
The  resistance  of  a  point  contact  is  determined  by 
the  Sharvin  formula:  Rs  =  4pl/3na2  [3],  where  a  is 
the  radius  of  the  contact  area  and  /  is  the  electron 
mean  free  path.  The  Sharvin  formula  is  valid  when 
a«l  For  a  contact  of  atomic  dimensions  one  could 
expect  resistance  to  be  of  the  order  of  Rs  ~  35  K£1 
Thus,  the  value  of  tip-sample  resistance  indicates 
approximately  the  tip-sample  distance  and/or  an 
area  of  tip-sample  contact. 

We  have  used  regular  STM  tips  (cut  gold  wires 
and  electrochemically  etched  tungsten  tips)  and 
metal  coated  fiber-tips  that  are  conventionally  used 
in  near-field  scanning  optical  microscopy  in  the 
SFDC  of  our  near-field  microscope.  Typical  shear- 
force  approach  curves  were  observed  together  with 
the  corresponding  tunneling  current  measured  be¬ 
tween  these  tips  and  a  freshly  prepared  thermally 
evaporated  gold  film  used  as  a  sample.  Depending 
on  the  mechanical  stiffness  of  the  STM  tip  we 
observed  both  contact  and  non-contact  van  der 
Waals  (or  capillary)  modes  of  SFDC  operation.  All 
the  data  indicate  that  metal  coated  fiber-tips  mech¬ 
anically  touch  samples  under  shear-force  distance 
control. 


2.  Experimental  setup 

Two  types  of  SFDC  have  been  used  in  our  ex¬ 
periments.  In  optical  SFDC  (Fig.  la)  a  regular 
STM  tip  or  metal  coated  fiber  tip  is  attached  to 
a  piezotube  that  provides  both  XYZ  motion  in  the 
nanometer  range  and  dithers  the  tip  at  its  resonant 
frequency.  The  vibration  amplitude  is  detected  with 
the  help  of  a  laser  beam  passing  by  the  fiber  near  its 
tip  in  a  direction  perpendicular  to  the  axis  of  the 
fiber.  This  beam  is  directed  onto  a  pinhole  covering 
a  photodiode.  The  component  of  optical  signal 
from  the  photodiode  at  the  frequency  of  tip  dither¬ 
ing  is  extracted  using  a  lock-in  amplifier  and 
directed  into  the  feedback  loop.  Typical  eigen- 
frequencies  of  fiber  tips  are  -3  kHz.  This  value  is 
limited  by  the  lowest  eigenfrequency  of  a  piezotube, 
which  is  4.5  kHz.  The  typical  amplitude  of  tip 


Fig.  1.  Schematic  diagrams  of  the  optical  (a)  and  non-optical 
(tuning  fork)  (b)  SFDC  setups. 

vibration  has  been  estimated  from  the  dithering 
voltage  and  the  Q-factor  of  the  fiber  tip  to  be 
5-10  nm. 

Alternatively,  a  non-optical  SFDC  [6,7]  utilizing 
a  quartz  crystal  tuning  fork  acting  as  a  resonant 
vibration  sensor  has  been  implemented.  A  fiber  tip 
is  mounted  on  one  of  the  prongs  of  the  fork 
(Fig.  lb)  to  serve  as  both  a  shear-force  pickup  and 
a  near-field  optical  probe.  The  crystal  is  driven  at 
its  resonance  frequency  (around  33  kHz)  by  an¬ 
other  piezo  crystal  mounted  next  to  the  tuning  fork. 
The  current  through  the  tuning  fork  is  converted  to 
a  voltage  and  is  fed  back  into  the  tip-sample  dis¬ 
tance  control  loop.  Typical  Q-factors  of  the  fork 
with  an  attached  fiber  tip  are  on  the  order  of 
1000-2000.  The  tip  vibrates  with  an  amplitude  of 
0.5  nm  during  approach  to  the  sample. 
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Fig.  3.  Time  scans  of  the  optical  shear-force  signal  (a)  and  the  tip-sample  current  (b)  measured  simultaneously  using  a  tungsten  STM  tip 
in  the  optical  SFDC  setup.  Electrochemically  etched  tungsten  tip  (eigenfrequency  v  =  4.288  kHz)  and  a  gold  film  sample  have  been  used. 
SFDC  level  is  0.90.  Period  of  oscillations  in  (a)  is  equal  to  v“  \  (c)  and  (d)  are  the  time  scans  of  a  current  between  a  chromium  coated  fiber 
tip  and  a  gold  sample  at  0.90  and  0.82  level  of  SFDC.  Period  of  oscillations  in  tunneling  current  is  again  equal  to  the  period  of  fiber  tip 
oscillations. 


steady.  Fig.  3  shows  time  scans  of  current  measured 
with  a  tungsten  STM  tip  and  a  chromium  coated 
fiber  tip  under  optical  SFDC.  The  spikes  of  current 
indicate  that  tips  are  actually  knocking  on  the 
sample  during  SFDC  feed  back.  The  tip-sample 
resistance  in  the  peaks  of  current  in  Fig.  3b  and 
Fig.  3d  is  measured  to  be  around  30  KQ.  This  value 
is  of  the  order  of  magnitude  of  the  resistance  of  an 
“atomic”  point  contact  according  to  Sharvin’s  for¬ 
mula.  However,  the  current  may  flow  through 
a  thicker  layer  of  contaminations  on  the  surface  of 
the  metal  film.  In  any  case,  all  these  data  support 
the  nonlinear  bending  force  hypothesis  [2]  for  the 
shear-force.  Yet  another  paper  supporting  the  idea 
of  mechanical  tip  contact  has  been  presented  in  this 
conference  by  Wei  and  Fann  [9]. 


The  drawback  of  our  optical  SFDC,  which  is 
common  to  most  other  types  of  optical  SFDC  is  the 
rather  large  amplitude  of  tip  vibration  in  free  space. 
Non-optical  SFDC  using  a  tuning  fork  as  a  sensi¬ 
tive  element  offers  an  improvement  in  this  respect. 
The  vibration  amplitude  of  a  tuning  fork  (x  ~ 
0.5  nm)  is  an  order  of  magnitude  smaller.  Despite 
its  stiff  spring  constant  (typically  20  000  N/m)  keep¬ 
ing  the  Q-factor  of  the  fork  high,  on  the  order  of 
1000-2000,  makes  it  dynamically  soft  (k/Q  - 
10  N/m),  although  it  is  still  stiffer  than  a  typical 
fiber  tip  probe  with  Q  =  50  and  k/Q  ~  2  N/m. 
Nevertheless,  an  actual  force  acting  on  the  lip 
(~kx/Q)  is  about  the  same  in  both  SFDC  setups 
The  approach  curve  of  a  tuning  fork  SFDC  setup  is 
shown  in  Fig.  4. 
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We  have  used  tapered  single  mode  optical  fibers 
coated  “head-on”  with  200  nm  of  gold-palladium 
in  our  tunneling  experiment.  We  were  able  to  run 
SFDC  feed  back  at  the  level  of  99.9%  of  the 
initial  undamped  shear-force  signal.  Moreover  the 
sensitivity  of  the  SFDC  has  been  increased  by  mov¬ 
ing  the  driving  signal  slightly  off  resonance  (by 
slightly  modifying  our  experimental  setup  in 
Fig.  lb).  Nevertheless,  we  always  see  spikes  in  the 
tunneling  current  of  the  kind  shown  in  Fig.  3.  This 
indicates  that  the  tip  is  knocking  on  the  sample  in 
the  tuning  fork  SFDC  setups.  One  has  to  be  careful 
running  this  experiment.  A  gold  or  gold-palladium 
coating  does  not  stick  really  well  to  the  fiber.  It  is 
easy  to  rip  off  the  coating  from  the  very  end  of  the 
fiber  tip.  On  a  number  of  occasions  we  did  not 
observe  a  tip-sample  current  under  SFDC  or  only 
observed  it  for  larger  tuning  fork  damping  (for 
example  at  70%  of  the  initial  shear-force  signal).  In 
all  these  cases  careful  study  of  the  tip  coating  re¬ 
vealed  tip  damage  as  shown  in  Fig.  5a.  This  dam¬ 
age  was  not  due  to  any  drawback  in  the  feed-back 
control  system.  The  vertical  fluctuations  of  the 

piezo  when  the  feed-back  control  was  on  were  less 
than  1  nm. 

It  should  be  noted  that  an  aluminum  coating 
does  stick  well  to  the  fiber.  It  is  possible  to  run 
tuning  fork  SFDC  with  much  larger  damping 
without  causing  physical  damage  to  the  coating. 
Fig.  5b  shows  the  very  end  of  an  aluminum  coated 
fiber  tip  after  scanning.  This  tip  was  used  in  the 
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“!!erd  /t]hat  did  not  reveal  any  tunneling  current  under 
UC  and  (b)  an  aluminum  coated  fiber  tip  after  scanning. 
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same  conditions  as  the  tip  in  Fig.  5a.  No  damage  is 
apparent  in  this  photo. 

All  our  data  indicate  that  the  tip  of  a  near  field 
microscope  is  knocking  on  the  sample  under  SFDC 
feed  back.  This  is  the  “universal”  physical  mecha¬ 
nism  behind  the  shear-force  in  all  environments. 
There  is  nothing  surprising  in  this  fact,  since  for 
example  AFM  is  a  contact  scanning  probe  micro¬ 
scope  as  well.  However,  we  should  point  out  the 
following  important  consequences  of  our  observa¬ 
tions: 

(1)  One  should  not  believe  shear-force  approach 
curves.  They  cannot  give  the  actual  distance  be¬ 
tween  the  tip  and  a  sample. 

(2)  It  is  much  safer  to  work  at  the  smallest  pos¬ 
sible  shear-force  damping,  since  in  this  case  there  is 
much  less  probability  of  damage  to  a  tip  aperture 
or  of  scratching  a  sample. 


if#! 


(3)  It  is  important  to  look  at  the  tip  aperture  not 
only  before  scanning  but  after  scanning  as  well. 
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Abstract 

We  present  conclusive  experimental  quantitative  evidence  of  the  resolution  limit  of  uncoated  optical  fiber  probes  in  the 
internal  reflection  mode.  Additionally,  we  present  a  new  technique  for  unambiguously  determining  the  resolution  of 
a  near-field  scanning  optical  microscope  without  topographical  influences.  A  sample  with  nearly  no  topography  but  with 
a  large  dielectric  step  junction  was  created  by  evaporating  a  thin  chromium  layer  on  a  silicon  wafer  and  subsequently 
cleaving  the  wafer.  The  cleaved  edge  is  then  scanned  over  the  step  junction,  allowing  a  quantitative  determination  of  the 
lateral  resolution  without  topographical  influences.  ©  1998  Elsevier  Science  B.V.  All  rights  reserved. 


1.  Introduction 

In  most  near-field  microscopes,  the  sides  of 
a  tapered  optical  fiber  are  coated  with  metal  to 
prevent  side  leakage  and  thus  confine  the  beam  to 
only  the  aperture’s  size,  since  side  leakage  would 
severely  degrade  the  lateral  spatial  resolution  in 
collection  or  reflection  near-field  scanning  optical 
microscopy  (NSOM).  Indeed,  when  the  fiber  dia¬ 
meter  goes  below  half  a  wavelength,  the  HEU 
mode  is  no  longer  confined  by  the  waveguide  and  is 
essentially  ejected  from  the  core  [1].  However,  if 
the  fiber  tip  is  used  both  as  source  and  detector  (the 
so-called  internal  reflection  mode),  the  double  pass 
through  the  aperture  effectively  produces  a  system 
where  the  point-spread  function  of  the  uncoated 
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aperture  is  traversed  twice,  eliminating  a  large  part 
of  the  unwanted  background  radiation. 

Several  claims  have  been  made  regarding  the 
resolution  of  uncoated  optical  fibers  in  the  internal 
reflection  configuration  [2-8],  but  these  are  at  best 
qualitative  estimates  that  may  be  seriously  exag¬ 
gerated  as  a  result  of  the  topographical  contrast 
present  in  the  samples.  Indeed,  a  recent  investiga¬ 
tion  [9]  of  the  internal  reflection  mode  revealed 
that  the  same  group’s  claim  of  better  than  70  nm 
(A/7)  resolution  [4,5]  with  uncoated  fibers  was 
simply  a  topographical  artifact.  A  fluorescence 
measurement  using  uncoated  fiber  probes  in  the 
internal-reflection  mode  estimates  the  lateral  res¬ 
olution  to  be  about  200  nm  (A/3)  by  measuring  the 
intensity  profile  at  a  fluorescent  film  edge  [10],  but 
this  estimate  may  also  be  affected  by  the  topogra¬ 
phy  of  the  fluorescent  film.  The  most  popular  tech¬ 
nique  involving  the  use  of  a  shadow  mask  [11] 
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suffers  from  a  large  topographical  element,  and  as 
such  claims  of  resolution  using  this  sample  remain 
questionable. 

The  aim  of  the  present  work  was  to  devise  an 
experimental  technique  of  determining  the  optical 
resolution  of  a  near-field  scanning  optical  micro¬ 
scope  without  the  influence  of  topographically  in¬ 
duced  intensity  changes  that  can  give  rise  to  a  false 
high  optical  resolution.  The  described  experimental 
technique  can  also  be  used  to  quantify  the  resolu¬ 
tion  of  other  near-field  microscopes,  including 
those  using  coated  probes  [12]  and  apertureless 
probes  [13]. 


2.  Experiment 

A  sample  of  virtually  no  topography  but  with 
a  well-defined  refractive  index  variation  was  con¬ 
structed  in  order  to  quantify  the  lateral  resolution 
of  the  uncoated  fiber  probe  in  the  internal-reflec¬ 
tion  configuration.  The  sample  we  have  used  was 
a  silicon  wafer  on  which  we  have  evaporated 
500  nm  of  chromium  and  then  cleaved  perpendicu¬ 
lar  to  the  (10  0)  plane  in  the  [Oil]  direction.  We 
chose  a  thickness  of  500  nm  so  that  the  chromium 
would  not  peel  after  cleaving.  The  cleaved  edge 


presents  a  nearly  atomically  flat  junction  between 
the  silicon  and  the  chromium  that  can  be  used  as 
a  step  junction  to  determine  the  lateral  resolution 
of  the  microscope.  The  sample  is  mounted  at 
a  slight  tilt  in  order  that  the  tip  glide  smoothly  near 
the  edge  of  the  chromium  region  rather  than  cata¬ 
strophically  fall  off.  The  measured  length  at  which 
the  optical  signal  rises  to  its  half  its  maximum  value 
gives  a  quantitative  measure  of  the  optical  resolu¬ 
tion.  A  shear-force  mechanism  using  tuning-fork 
distance  regulation  [14]  was  used  to  insure  that  the 
tip  was  always  in  the  near  field  of  the  sample.  The 
vibration  amplitude  was  set  at  less  than  1  nm  to 
insure  that  it  did  not  contribute  to  any  loss  of 
resolution.  Our  tip  diameter  used  was  less  than 
100  nm  determined  through  SEM  imaging  before 
and  after  scanning  the  junction. 


3.  Results  and  conclusions 

In  Fig.  1,  we  show  a  2  x  2  pm  scan  of  the  cleaved 
Si-wafer  edge  where  a  region  of  chromium  lies 
directly  at  the  edge.  The  surface  is  representative  of 
the  shear-force  topographical  image  while  the 
shading  is  rendered  by  the  near-field  signal.  A  sud¬ 
den  change  in  slope  occurs  when  the  tip  begins  to 


Fig.  1.  2  X  2  Jim  scan  of  the  silicon-chromium  interface.  The  surface  topography  is  given  by  the  shear  h me  itgit*).  wM* 
rendered  by  the  optical  signal.  The  500  nm  chromium  region  is  clearly  resolved. 
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Fig.  2.  Line  scan  of  silicon-chromium  cleaved-edge  interface.  The  full-width  half-maximum  reveals  a  resolution  of  300  nm. 


fall  off  the  edge  of  the  chromium  region;  otherwise, 
there  is  no  observable  topography  within  the  noise 
limit  ( <  3  nm)  of  our  microscope.  The  entire  height 
variation  in  Fig.  1  is  only  30  nm,  due  to  the  tilt  of 
the  sample.  As  shown  in  the  figure,  the  500  nm 
region  of  the  chromium  is  clearly  resolved.  A  line 
scan  of  the  optical  signal  through  the  region  is 
shown  in  Fig.  2.  From  this  scan,  we  can  quantitat¬ 
ively  deduce  the  resolution  of  the  uncoated  fiber 
probe  in  the  internal-reflection  mode  to  be  about 
300  nm  (2/2.6)  (from  the  full-width  at  half-max¬ 
imum  taken  from  the  baseline  of  the  step  transition 
before  the  edge),  which  is  significantly  better  than 
the  diffraction  limit  of  390  nm  for  our  780  nm 
wavelength  source.  We  should  also  note  that  our 
best  conventional  microscope  was  not  able  to  suc¬ 
cessfully  image  the  chromium  region.  The  vari¬ 
ations  before  and  after  the  chromium  region  are 
likely  due  to  interference  effects  at  the  step  junction 
and  resemble  a  far-field  microscope  diffraction  pat¬ 
tern  generated  at  a  step  junction.  A  large  back¬ 
ground  field  generated  by  the  uncoated  fiber  tip 
often  gives  rise  to  such  imaging  artifacts  at  a  sharp 
junction.  These  effects  may  be  eliminated  by  using 


a  coated  fiber  tip  since  such  a  source  has  a  well- 
defined  aperture. 

In  conclusion,  we  present  quantitatively  the  res¬ 
olution  limit  of  uncoated  fiber  tips  used  in  the 
internal-reflection  configuration.  The  results  agree 
well  with  recent  experiments  [9]  that  suggest  the 
uncoated  fiber  configuration  is  nearly  diffraction 
limited.  Furthermore,  we  describe  a  new  experi¬ 
mental  technique  for  characterizing  the  optical  res¬ 
olution  in  a  near-field  microscope  without  the 
topographically-induced  artifacts  present  in  cur¬ 
rent  resolution  standards. 
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Near-field  optical  images  of  electromagnetic  mode  structures  in  two-dimensional  photonic  and  plasmonic 
crystals  are  presented.  Interference  patterns  of  plasmon  Bloch  waves  in  a  plasmonic  crystal  have  been  ob¬ 
served.  Fourier  analysis  of  the  photonic  crystal  images  allows  the  reconstruction  of  the  refractive  index  spatial 
distribution  in  the  crystal.  A  similar  procedure  performed  on  images  of  a  plasmonic  crystal  allows  the  recovery 
of  parameters  of  the  surface  plasmon  interaction  with  the  periodic  grooves  that  constitute  the  plasmonic  crystal 
structure.  [S0163-1829(99)13203-X] 


INTRODUCTION 

The  optical  properties  of  spatially  periodic  dielectric 
structures  called  photonic  crystals  have  attracted  much  recent 
interest.  The  analogy  between  electromagnetic  wave  propa¬ 
gation  in  multidimensionally  periodic  structures  and  electron 
wave  propagation  in  real  crystals  has  proven  to  be  a  valuable 
one.  Electron  waves  traveling  in  the  periodic  potential  of  a 
crystal  are  arranged  into  energy  bands  separated  by  gaps  in 
which  propagating  states  are  prohibited.  Analogous  band 
gaps  exist  when  electromagnetic  waves  propagate  in  a  peri¬ 
odic  dielectric  structure  with  a  period  comparable  to  the 
wavelength.  These  frequency  gaps  are  referred  to  as  “pho¬ 
tonic  band  gaps.”  A  demonstration  of  band  structure  effects 
in  photonic  crystals  was  in  the  microwave  regime.2  Recently, 
these  effects  have  been  measured  in  the  near-infrared  and 
visible  regions.3  Photonic  crystals  can  have  a  profound  im¬ 
pact  in  many  areas  of  pure  and  applied  physics.  For  example, 
inside  a  photonic  band  gap,  optical  modes,  spontaneous 
emission,  and  zero  point  fluctuations  are  all  absent.1  From  a 
practical  point  of  view  photonic  crystals  can  be  designed  to 
transmit  or  reflect  light  within  a  specific  range  of  frequen¬ 
cies,  and  their  properties  may  be  tunable  by  modification  of 
their  periodicity  or  refraction  index. 

Until  recently  photonic  crystals  have  been  studied  by 
measuring  their  bulk  optical  properties.  With  the  invention  of 
the  near-field  scanning  optical  microscope4  (NSOM)  local 
microscopic  properties  of  photonic  crystals  become  acces¬ 
sible  to  experimental  measurement.  An  application  of 
NSOM  techniques  to  characterize  the  local  density  of  photon 
states  and  the  electromagnetic  mode  structure  of  a  two- 
dimensional  photonic  crystal  has  been  reported  in  Ref.  5.  In 
this  work  a  photonic  crystal  was  locally  excited  by  light 
leaving  a  subwavelength  aperture  at  the  end  of  a  metal- 
coated  fiber  tip.  The  light  transmitted  by  the  crystal  was  col¬ 
lected  in  the  far  field  with  a  microscope  objective.  The  nu¬ 
merical  aperture  (NA)  of  the  objective  was  varied. 
Depending  on  the  value  of  the  NA,  the  transmitted  light  with 
different  spatial  momenta  was  integrated  and  recorded  as  a 
function  of  the  position  of  the  fiber  tip.  The  contrast  in  the 
images  obtained  depended  strongly  on  the  region  in  the  pho- 
ton  momentum  space  probed. 

Theoretically  this  situation  can  be  understood  as  follows: 
the  photonic  crystal  was  excited  by  a  photon  source  that  may 


be  idealized  as  a  8  function  in  coordinate  space.  The  trans¬ 
mitted  light  was  analyzed  as  a  function  of  the  absolute  value 
of  its  momentum  and  the  spatial  position  of  the  light  source. 
A  complicated  theoretical  model  has  to  be  used  in  order  to 
analyze  this  experiment. 

In  the  present  work  we  study  an  alternative  approach  to 
the  near-field  optical  characterization  of  photonic  crystals. 
We  suggest  exciting  the  crystal  with  an  incident  laser  beam. 
This  may  be  idealized  as  an  excitation  by  a  S  function  in 
momentum  space.  We  use  NSOM  in  the  transmission  mode 
to  measure  directly  the  local  electromagnetic  mode  structure 
on  the  surface  of  a  photonic  crystal  in  coordinate  space.  The 
NSOM  images  are  Fourier  analyzed  in  order  to  extract  infor¬ 
mation  on  the  band-gap  structure  and  refractive  index  distri¬ 
bution  in  the  photonic  crystal.  We  show  this  approach  to  be 
straightforward  and  simple. 

We  also  apply  our  technique  to  the  study  of  a  two- 
dimensional  “plasmonic  crystal”:  a  gold-coated  rectangular 
bigrating,  which  possesses  properties  similar  to  the  proper¬ 
ties  of  a  photonic  crystal.  Far-field  optical  properties  of  simi¬ 
lar  plasmonic  crystals  have  been  studied  very  recently6  but 
no  near-field  studies  of  such  crystals  have  been  reported  to 
date. 

The  surface  plasmon  (SP)  is  a  fundamental  electromag¬ 
netic  excitation  mode  of  a  metal-dielectric  interface.7  The  SP 
is  free  to  propagate  along  the  metal  surface,  but  its  field 
decays  exponentially  in  both  media  in  the  direction  perpen¬ 
dicular  to  the  interface.  Thus  the  spatial  distribution  of  the 
SP  field  is  inaccessible  to  far-field  optical  techniques.  The 
development  of  NSOM  has  created  ways  to  probe  the  SP 
field  directly.  SP  propagation  along  the  interface  has  been 
imaged  in  Ref.  8.  SP  scattering  by  in  situ  created  individual 
surface  defects  has  been  studied  and  some  prototype  two- 
dimensional  optical  elements  able  to  control  SP  propagation 
have  been  demonstrated.9,10 

A  periodic  array  of  defects  in  the  metal-dielectric  inter¬ 
face  should  exhibit  the  properties  of  a  two-dimensional 
“plasmonic  crystal”  when  the  periodicity  of  the  structure  is 
comparable  with  the  wavelength  of  the  SP.  A  complicated 
system  of  band  gaps  has  indeed  been  measured  in  Ref.  6 
using  different  bigratings.  In  the  present  paper  we  report  our 
measurements  of  the  local  electromagnetic  mode  structure  of 
a  plasmonic  crystal  and  derive  some  parameters  of  SP  inter¬ 
action  with  defects  that  constitute  the  bigrating.  We  present  a 
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FIG.  1.  Far-field  transmission  image  of  the  photonic  crystal  structure  taken  with  an  optical  microscope  using  bottom  face  illumination 
with  white  light. 


similar  description  of  the  results  obtained  with  two- 
dimensional  photonic  and  plasmonic  crystals. 

EXPERIMENTAL  SETUP  AND  SAMPLES 

We  have  studied  two  different  kinds  of  photonic  band-gap 
structures.  The  first  sample  studied  is  similar  to  the  photonic 
crystal  used  in  Ref.  5.  It  is  a  fiber  bundle  obtained  from 
Schott  Glass  Technologies,  Inc.  made  of  two  glasses  with 
similar  indices  of  refraction.  The  bundle  is  a  two- 
dimensional  rectangular  lattice  of  glass  with  a  higher  refrac¬ 
tive  index  embedded  in  a  matrix  glass  with  a  lower  refractive 
index.  The  period  of  the  structure  is  approximately  3  /im. 
Figure  1  shows  an  image  of  a  top  face  of  the  sample  taken 
with  a  conventional  optical  microscope  using  bottom  face 
illumination  with  a  white  light  source.  The  image  shows 
good  periodicity  of  the  structure  although  some  dislocations 
and  blemishes  are  present.  Both  faces  of  the  sample  are  pol¬ 
ished.  The  thickness  of  the  sample  is  about  100  mm. 

The  second  sample  is  a  gold-coated  two-dimensional  rect¬ 
angular  calibration  grating  with  a  period  of  about  1  /im  ob¬ 
tained  from  Digital  Instruments,  Inc.  An  atomic  force  micro¬ 
scope  image  of  the  sample  surface  is  shown  in  Fig.  2. 
According  to  the  manufacturer,  a  30-nm  layer  of  gold  has 
been  deposited  by  thermal  evaporation  on  top  of  a  glass  grat¬ 
ing  coated  with  800  nm  of  aluminum  and  4  nm  of  chromium. 
In  the  absence  of  grooves  such  a  “sandwich”  structure  sup¬ 
ports  surface  plasmon  propagation  with  a  spectrum  close  to 
the  spectrum  of  SPs  on  the  surface  of  an  infinitely  thick  flat 
gold  sample.  (The  decay  length  of  the  SP  field  inside  gold  is 
about  /  =  \/27re)^.7  At  a  wavelength  of  633  nm  this  gives 
/=36nm.)  The  grooves  scatter  the  SP  field  in  the  plane  of 
propagation.  Also  they  couple  the  SP  field  to  photons  in  free 


space.  Since  the  period  of  the  structure  is  comparable  with 
the  wavelength  of  plasmons:  Xp  =  \[(e  Au+ l)/€Au]1/2 
=  590  nm,  where  X  =  633  nm  is  the  wavelength  of  light  in 
air,  this  sample  must  behave  as  a  “plasmonic  crystal.”  There 
should  be  forbidden  gaps  in  the  SP  momentum  space.  The 
spatial  topology  of  such  SP  forbidden  gaps  in  a  similar  rect¬ 
angular  bigrating  has  been  studied  in  Ref.  6. 

Our  NSOM  setup  is  described  in  detail  in  Ref.  10.  It  al¬ 
lows  us  to  measure  simultaneously  surface  topography  using 
optical  shear-force  feedback,4  and  the  spatial  distribution  of 
the  optical  signal  in  near-field  proximity  to  the  sample  sur¬ 
face.  In  order  to  study  the  samples  described,  our  experimen¬ 
tal  setup  has  been  modified  as  follows  [Figs.  3(a)  and  3(b)]. 
The  photonic  crystal  from  Fig.  1  has  been  studied  using 
NSOM  working  in  a  transmission  mode  as  shown  in  Fig. 


FIG.  2.  Atomic  force  microscopic  image  of  the  plasmonic  crys¬ 
tal  surface  topography. 
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FIG.  3.  Experimental  setups  for  studying  optical  properties  of  (a)  the  photonic  crystal  and  (b)  the  plasmonic  crystal. 


3(a).  The  angle  of  the  bottom  face  illumination  was  varied 
and  the  field  distribution  on  the  top  face  was  studied  as  a 
function  of  this  angle. 

The  SPs  on  the  surface  of  the  plasmonic  crystal  from  Fig. 
2  were  excited  in  reflection  [Fig.  3(b)].  It  is  impossible  to 
excite  a  SP  wave  on  the  flat  gold-air  interface  by  illumina¬ 
tion  of  the  surface  from  air.  At  a  given  frequency  the  mo¬ 
mentum  of  photons  in  air  is  smaller  than  the  momentum  of 
plasmons.  The  presence  of  the  periodically  spaced  grooves 
leads  to  a  modification  of  the  momentum  conservation  law: 
integer  multiples  of  the  reciprocal  lattice  vectors  of  the  big¬ 
rating  can  be  added  to  the  photon  and  plasmon  momenta. 
This  leads  to  the  possibility  of  SP  excitation  by  photons  il¬ 
luminating  the  bigrating  at  certain  resonant  angles.  Accord¬ 
ingly,  the  angle  of  sample  illumination  was  varied  in  the 
experimental  geometry  shown  in  Fig.  3(b)  in  order  to 
achieve  SP  excitation.  P-polarized  light  from  a  He-Ne  laser 
has  been  used.  The  SP  field  distribution  has  been  measured 
using  an  uncoated  tapered  fiber  tip  in  our  NSOM.  Although 
the  resolution  of  a  NSOM  is  higher  when  metal-coated  tips 
are  used,4  an  uncoated  fiber  tip  introduces  a  much  smaller 
perturbation  in  the  SP  field  distribution.10  This  results  from 
the  much  smaller  difference  between  the  dielectric  constants 
of  the  glass  tip  and  air  (erip-l)  in  comparison  with  the 
difference  between  the  dielectric  constants  of  air  and  metal. 

PHOTONIC  CRYSTAL  IMAGING 

Near-field  optical  images  of  the  photonic  crystal  sample 
obtained  at  different  angles  of  excitation  are  shown  in  Figs. 
4(a)-4(c).  The  topographical  images  taken  simultaneously 
show  no  features  higher  than  10  nm  and  do  not  correlate  with 
the  optical  ones.  The  optical  images  in  Fig.  4  exhibit  general 
periodicity  of  the  refractive  index  structure  (a  3-/zm  period  in 
both  directions).  Some  distortions  in  the  images  may  be  at¬ 
tributed  to  the  dislocations  observed  in  Fig.  1. 

The  main  feature  of  these  images  is  an  appearance  of 
additional  periodicity  and  a  richer  spatial  frequency  spectrum 
at  larger  excitation  angles.  This  feature  is  quite  clear  from 
the  Fourier  spectra  presented  in  Figs.  5(a)  and  5(b).  Such 
behavior  is  easy  to  understand  when  the  general  theory  of  a 
photonic  band-gap  material  is  applied  to  our  particular  case 
(an  extensive  theoretical  description  of  two-dimensional 
photonic  crystals  can  be  found  in  Ref.  11). 


When  a  monochromatic  electromagnetic  wave  of  fre¬ 
quency  o)  propagates  in  a  dielectric  medium  with  a  periodic 
dielectric  constant  e(r)  =  60-f-61(r)  (e0  is  an  average  dielec¬ 
tric  constant  of  a  medium),  the  classical  wave  equation  may 
be  written  in  a  form  resembling  the  Schrodinger  equation: 

-  V2E+  V(  VE)  -  fc)2€,(r)E/c2=  €0<o2E/c2.  (1) 

In  this  equation  fj(r)  plays  a  role  analogous  to  the  peri¬ 
odic  potential  [/(r)  in  the  Schrodinger  equation  for  an  elec¬ 
tron  in  a  crystal.  In  the  two-dimensional  crystal  under  con¬ 
sideration  the  photon  motion  may  be  decomposed  into  the 
motion  in  a  periodic  potential  (perpendicular  to  the  fiber 
bundle  axis)  and  the  free  motion  along  the  structure  with 
momenta  kx  and  so  that  *#+*£  =  €0<o2Ic2.  Thus  larger 
excitation  angles  correspond  to  a  larger  momentum  of  mo¬ 
tion  in  a  periodic  potential.  Consequently,  higher  harmonics 
may  be  expected  in  the  near-field  image.  To  some  extent  our 
data  are  similar  to  the  images  of  metal  or  semiconductor 
crystals  obtained  in  scanning  tunneling  microscopy  with  dif¬ 
ferent  tunneling  voltages  applied.  In  both  cases  one  can  vary 
the  particle  (electron  or  photon)  energy  in  the  periodic  po¬ 
tential  as  a  free  parameter. 

The  intensities  of  different  spectral  peaks  in  Fig.  5  contain 
information  about  the  periodic  part  of  the  dielectric  constant 
^i(r)  of  the  photonic  crystal.  The  Bloch  theorem  can  be  used 
in  order  to  extract  this  information.  Exact  theory11  gives 
slightly  different  equations  for  different  polarization  states  of 
the  electromagnetic  field  in  the  two-dimensional  photonic 
crystal.  On  the  other  hand,  regular  individual  optical  fibers 
(not  polarization  maintaining  ones)  similar  in  design  to  the 
fibers  that  constitute  our  photonic  crystal  (Fig.  I)  scramble 
the  polarization  of  incoming  light.  Generally,  NSOM  fiber 
tips  also  do  not  preserve  the  polarization  of  the  light,12  so  we 
can  disregard  this  difference  and  consider  the  field  to  be 
scalar.  In  this  case,  the  wave  equation  looks  exactly  like  the 
Schrodinger  equation  for  an  electron  in  a  periodic  potential 
U=(o2€ i  / c 2  and  we  can  directly  use  the  standard  represen¬ 
tation  of  the  Schrodinger  equation  in  momentum  space11: 

(tf+k21-e0a>2/c2)Ak±  +  'Z  -r.  =  0.  (2) 
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FIG.  4.  Transmission  near-field  optical  images  of  the  photonic 
crystal  taken  at  different  angles  of  bottom  face  illumination-  (a)  0° 
(b)  7 °,  (c)  15°. 

where  f/(r)  =  EGt/G^/G  r  is  the  periodic  potential,  G  are  the 
reciprocal  lattice  vectors,  and 


E=eik^ |  Akx_G*-G  '  (3) 

is  the  Bloch  wave  solution  of  the  equation. 

Usually  it  is  reasonable  to  assume  that  the  coefficients  Ak 
are  very  small  for  large  absolute  values  of  k.  We  also  see  this 
from  our  experimental  measurements  in  Fig.  5.  This  reduces 
the  Schrodinger  equation  (2)  to  a  small  number  of  linear 
equations  relating  Ak  and  UG . 


Spectrum  2D 


5.  Images  (a)  and  (b)  are  the  Fourier  spectra  of  the  images 
from  Figs.  4(a)  and  4(b),  respectively.  Additional  periodicity  and  a 
richer  spatial  frequency  spectrum  appear  at  larger  excitation  angles. 
Spectral  peaks  near  the  first-  and  second-order  reciprocal  lattice 
vectors  are  shown  by  arrows  in  (b). 

We  have  performed  a  one-dimensional  reconstruction  of 
the  photonic  crystal  refractive  index  as  follows.  The  physical 
value  measured  in  the  experiment  is  the  local  electromag¬ 
netic  field  intensity  on  the  face  of  the  photonic  crystal.  If  a 
one-dimensional  Bloch  wave  propagates  in  the  crystal  along 
the  x  direction,  its  field  distribution  can  be  written  as  E 
~  ^n^ne  *  ^')e‘2trnx/a ,  whgre  a  js  period  of  the  crystal 
in  the  direction  of  propagation  and  n  is  an  integer.  An  aver¬ 
aged  field  intensity  is  equal  to 


7-1/22  AmAn  cos[2'7t(/i  tn)xl dj.  (4) 

mn 

Thus  measurements  of  the  intensities  of  spectral  peaks  at 
iTTnla  in  Fig.  5  allow  us  to  recover  the  coefficients  A„  of  the 
Bloch  wave.  The  first  three  coefficients  have  been  measured. 
Substitution  of  these  A„  in  Eq.  (2)  allows  us  to  obtain  three 
coefficients  Un  of  the  Fourier  spectrum  of  the  crystal  poten¬ 
tial  in  this  direction:  U(x)  =  lnUnei2mx,a.  The  results  of  our 
measurements  are  represented  in  Fig.  6,  which  shows  varia¬ 
tions  of  the  refractive  index  of  the  photonic  crystal  along  one 
direction.  The  amplitude  of  these  variations  agrees  well  with 
the  indices  of  the  bulk  unprocessed  glasses  used  to  make  the 
crystal,  which  were  measured  by  Schott  Glass  Technologies, 
Inc.  at  590  nm  to  be  1.657  and  1.676.  Thus  the  simple  ap- 


2458 


SMOLYANINOV,  ATIA,  AND  DAVIS 


PRB  59 


FIG.  6.  Reconstructed  refractive  index  variations  in  the  photonic 
crystal  along  one  direction. 

proach  introduced  above  allows  a  reasonable  description  of 
the  photonic  crystal  properties. 

PLASMONIC  CRYSTAL  STUDY 

A  surface  plasmon  can  be  excited  on  the  surface  of  the 
rectangular  bigrating  (Fig.  2)  at  a  number  of  resonant  angles. 
These  angles  a„  are  determined  by  the  quasimomentum  con¬ 
servation  law:  k  sin  an=kp+2mi/a,  where  kp  is  the  SP  wave 
vector  on  the  flat  metal  surface,  a  is  the  period  of  the  struc¬ 
ture  in  the  direction  of  propagation,  and  n  is  an  integer.  The 
angular  dependence  of  the  optical  signal  measured  at  the 
output  of  the  fiber  tip  of  our  NSOM  near  one  of  the  resonant 
angles  is  shown  in  Fig.  7.  Resonant  excitation  of  SP  has  been 
clearly  detected.  The  angular  width  of  the  resonant  curve  is 
narrower  than  in  Ref.  6.  This  may  be  attributed  to  the  better 
surface  quality  of  the  bigrating  used  in  our  experiments. 

We  simultaneously  measured  topography  of  the  sample 
and  the  optical  field  distribution  with  resonant  SP  excitation. 
Some  examples  of  these  images  are  shown  in  Figs.  8(a)- 
8(c).  The  periodicity  properties  of  the  optical  images  are 
much  more  complicated  than  those  of  a  topographical  image. 
This  is  evident  from  the  Fourier  image  of  Fig.  8(c)  shown  in 
Fig.  9.  The  brightest  spectral  peak  in  this  image  corresponds 
to  a  spatial  period  of  3.5  /xm.  Our  calculations  show  this 
value  to  be  approximately  equal  to  half  of  the  largest  effec- 


FIG.  7.  Angular  dependence  of  the  NSOM  optical  signal  near 
one  of  the  SP  resonances. 


FIG.  8.  Topography  (a)  and  optical  near-field  distributions  in  the 
plasmon  crystal  measured  on  different  scales  (b),  (c)  with  surface 
pJasmon  excitation. 

tive  wavelength  of  SP  propagating  along  the  surface  of  the 
plasmonic  crystal:  X  *  =  2 ir/k*  =  2 7r/(kp  +  2  )  =  7.26 

fixn  at  n  —  —  2.  This  indicates  that  we  observe  interference 
pattern  of  the  plasmon  Bloch  waves  in  the  image  in  Fig.  8(c). 
Complicated  interference  patterns  formed  by  SPs  on  the  flat 
gold  surface  had  been  observed  in  Ref.  14.  They  were  attrib¬ 
uted  to  the  presence  of  surface  defects.  The  surface  quality  of 
our  plasmonic  crystal  structure  is  very  good  (Fig.  2),  never¬ 
theless  some  surface  defects  can  be  seen.  So  the  appearance 
of  SP  Bloch  wave  interference  patterns  in  the  images  is  natu¬ 
ral. 
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FIG.  9.  Fourier  spectrum  of  the  image  in  Fig.  8(b).  The  spectral 
peak  at  a  spatial  period  of  3.5  fim  is  shown  by  an  arrow. 

NSOM  allows  us  to  modify  a  surface  under  investigation. 
Direct-write  lithography  may  be  performed  by  coupling  UV 
light  into  a  fiber  tip9  or  by  delicate  mechanical  touching  of 
the  sample  surface.15  Mechanical  surface  modification  is  es¬ 
pecially  easy  in  the  case  of  a  soft  sample  such  as  a  gold 
surface.  We  have  verified  experimentally  that  artificially  cre¬ 
ated  defects  lead  to  the  appearance  of  very  pronounced  long- 
period  interference  patterns  in  the  vicinity  of  defects.  The 
results  of  such  an  experiment  are  shown  in  Fig.  10.  A  defect 
has  been  created  in  the  middle  of  the  NSOM  field  of  view  on 
top  of  the  periodic  plasmon  crystal  structure,  which  has  also 
been  slightly  modified.  A  shadowlike  phenomenon  is  seen 
just  behind  the  defect  (the  direction  of  SP  propagation  is 
shown  by  an  arrow).  A  set  of  stripes  with  a  period  on  the 
order  of  3  fxm  has  become  very  pronounced.  These  results 
indicate  that  we  are  indeed  dealing  with  the  interference  of 
SP  Bloch  waves. 

The  SP  field  distribution  and  band-gap  structure  of  a  plas- 
monic  crystal  can  be  calculated  analytically  in  the  limit  of 
small  roughness  if  the  surface  topography  is  known.16  In  our 
case  the  surface  roughness  cannot  be  considered  to  be  small 
since  the  depth  of  the  grooves  in  Fig.  2  is  comparable  with 
the  thickness  of  the  gold  film.  Rather  complicated  numerical 
calculations  are  required  to  compare  the  results  of  our  mea¬ 
surements  with  an  exact  theory  also  presented  in  Ref.  16. 
Such  a  comparison  is  desirable  and  could  give  a  good  test  of 
the  electrodynamics  of  metals  on  a  submicrometer  scale. 

A  simpler  approach  is  to  use  the  concept  of  effective  plas¬ 
mon  refractive  index  introduced  in  Ref.  10.  SP  phase  veloc¬ 
ity  depends  on  the  metal  film  thickness.  The  effective  refrac¬ 
tive  index  of  a  surface  defect  may  be  introduced  as  a  ratio  of 
SP  phase  velocities  on  the  defect  and  on  the  background 
metal  film  area.  Use  of  this  concept  makes  the  description  of 
photonic  and  plasmonic  crystals  similar.  Effective  SP  refrac¬ 
tive  index  distribution  over  the  crystal  surface  may  be  de¬ 
rived  using  the  procedure  described  above  for  the  case  of  a 
photonic  crystal.  Unfortunately,  SP  Bloch  wave  interference 
makes  the  situation  more  complicated  since  they  contribute 
substantially  to  the  Fourier  spectrum  of  the  SP  field  distribu- 


FIG.  10.  Topography  (a)  and  optical  near-field  distribution  (b) 
around  artificially  created  defect  on  the  surface  of  the  plasmon  crys¬ 
tal  structure.  The  direction  of  SP  propagation  is  shown  by  an  arrow. 

tion  and  the  amplitude  of  the  reflected  waves  is  generally 
unknown.  Only  an  order  of  magnitude  estimate  of  the  effec¬ 
tive  refractive  index  can  be  derived  from  Fig.  9.  It  follows 
from  Eq.  (2)  that  o)2ex  i c2 ~  (2tt/ a)2A„+1  / An .  This  leads  to 
Art~(\/<2)2An  +  1  /An~0.05  as  an  estimate  of  the  amplitude 
of  the  effective  SP  refractive  index  variations.  This  number 
is  in  general  agreement  with  the  values  of  SP  refractive  in¬ 
dices  of  different  surface  defects  derived  in  Ref.  10. 

CONCLUSION 

The  main  results  of  this  paper  may  be  summarized  as 
follows. 

An  alternative  method  for  studying  photonic  crystals  us¬ 
ing  NSOM  has  been  introduced.  In  a  previous  study  5  the 
photonic  crystals  were  excited  by  light  from  a  NSOM  tip. 
This  may  be  idealized  as  an  excitation  with  a  8  function  in 
coordinate  space.  The  transmitted  light  distribution  was  ana¬ 
lyzed  as  a  function  of  the  absolute  value  of  its  momentum. 

In  the  present  work  photonic  crystals  were  excited  with  an 
incident  laser  beam.  This  may  be  idealized  as  an  excitation 
by  a  8  function  in  momentum  space.  NSOM  images  are  Fou¬ 
rier  analyzed  in  order  to  extract  information  on  refractive 
index  distributions  in  the  photonic  crystals. 

Interference  patterns  produced  by  SP  Bloch  waves  in  the 
plasmonic  crystal  have  been  directly  observed. 
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A  technique  to  characterize  the  local  nonlinearity  of  lead  zirconate  titanate  piezoceramic  (PZT)  has 
been  developed.  Near-field  optical  microscopy  has  been  used  to  image  variations  in  local  optical 
second  harmonic  generation  from  the  PZT  surface.  Individual  crystalline  grains  and  grain 
boundaries  are  the  main  features  visible  in  the  images.  The  technique  allows  us  to  determine  the 
local  poling  direction  of  individual  submicrometer  size  crystalline  grains  of  ceramic  by  near-field 
second  harmonic  imaging  at  different  angles  of  incidence  and  polarization  states  of  fundamental 
excitation  light.  ©  1999  American  Institute  of  Physics.  [S0003-695 1(99)005 14-8] 


The  ability  of  ferroelectric  materials  such  as 
Pb(ZrtTi]_x)03  (PZT)  piezoceramics  to  switch  from  one 
stable  polarization  state  to  another  forms  the  basis  of  a  new 
thin  film  technology  for  data  storage.1  Thin  PZT  films  are 
used  in  prototype  nonvolatile  ferroelectric  random-access 
memory  (NVFRAM)  and  dynamic  random-access  memory 
(DRAM)  devices.  PZT  materials  also  have  numerous  appli¬ 
cations  in  actuators,  transducers,  resonators,  and  sensors. 
Crucial  parameters  of  piezoceramic  performance  in  different 
applications  are  hysteresis,  nonlinearity,  polarization  reten¬ 
tion  (or  the  loss  of  polarization),  creep,  etc.  For  proper  ap¬ 
plication  of  the  material,  it  is  imperative  to  understand  this 
nonlinear  behavior  on  a  microscopic  level.  Optical  micros¬ 
copy  of  second  harmonic  generation  (SHG)  in  PZT  is  a  natu¬ 
ral  tool  to  study  this  local  nonlinearity. 

Far  field  observations  of  SHG  in  PZT  ceramic  near  the 
ferroelectric  phase  transition2  demonstrated  a  sharp  drop  in 
SHG  in  the  paraelectric  state.  The  second  harmonic  intensity 
versus  temperature  curve  corresponds  qualitatively  to  the 
temperature  dependence  of  the  spontaneous  polarization, 
which  is  an  order  parameter  in  the  transition  to  the  polar 
phase.  PZT  ceramic  is  a  strongly  scattering  medium  consist¬ 
ing  of  individual  submicrometer-size  crystallites.  The  elastic 
scattering  lengths  of  fundamental  and  SH  light  in  PZT  is  of 
the  order  of  the  characteristic  size  of  a  crystallite.  Theory  of 
SHG  in  strongly  scattering  media  shows3  that  generation  of 
SH  light  in  PZT  should  occur  in  a  thin  layer  (of  the  order  of 
an  elastic  scattering  length)  near  the  illuminated  interface. 
Thus,  near-field  second  harmonic  microscopy4  is  ideally 
suited  for  studies  of  PZT  nonlinearity  and  poling  quality  at 
the  microscopic  level  (the  level  of  individual  crystallites  and 
crystallite  boundaries).  Its  main  advantages  in  comparison 
with  the  other  scanning  probe  techniques  such  as  recently 
developed  piezoresponse  atomic  force  microscopy5  are  the 
possibility  of  fast  time-resolved  measurements,  and  substan¬ 
tially  smaller  perturbation  of  the  sample  under  investigation 
caused  by  the  optical  probe. 

We  report  second  harmonic  imaging  of  the  surface  of 
PZT  ceramic  obtained  in  a  near-field  microscopy  setup  (Fig. 

1)  using  a  Ti:sapphire  laser  system  consisting  of  an  oscillator 
and  a  regenerative  amplifier  operating  at  810  nm  (repetition 
rate  up  to  250  kHz,  100  fs  pulse  duration,  and  up  to  10  /xJ 
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pulse  energy).  The  SHG  from  the  surface  of  a  PZT  sample 
was  excited  at  an  angle  of  incidence  around  60°  by  a  weakly 
focused  beam  (illuminated  spot  diameter  on  the  order  of  50- 
100  jxm)  of  the  Tirsapphire  laser.  The  excitation  power  at  the 
sample  surface  was  kept  well  below  the  ablation  threshold. 
The  local  SHG  was  collected  using  an  uncoated  adiabatically 
tapered  fiber  tip,  which  is  drawn  at  the  end  of  a  single  mode 
fiber  by  the  standard  heating  and  pulling  procedure.  The  fiber 
tip  was  scanned  over  the  sample  surface  with  a  constant  tip- 
surface  distance  of  a  few  nanometers  using  shear  force  feed¬ 
back  control  described  elsewhere.4  Therefore,  surface  topog¬ 
raphy  images  were  obtained  with  a  resolution  on  a 
nanometer  scale,  while  simultaneously  recording  the  SH 
near-field  image.  The  SH  signal  was  measured  with  a  photo¬ 
multiplier  and  gated  electronics.  The  SH  signal  at  every 
point  of  the  image  was  averaged  over  30-100  ms.  The  char¬ 
acteristic  SH  photon  counting  rate  was  on  the  order  of  one 
photon  count  per  30  laser  pulses.  Thus,  the  image  acquisition 
time  was  on  the  order  of  20-30  min.  This  is  approximately 
an  order  of  magnitude  improvement  with  respect  to  our  pre¬ 
viously  reported  work4  where  a  Nd:yttrium-aluminum- 
gamet  (YAG)  laser  was  used  as  a  source  of  fundamental 
light. 
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FIG.  1.  Schematic  view  of  the  near-field  optical  microscope  for  second 
harmonic  generation  studies. 
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FIG.  2.  Simultaneously  measured  (a)  shear-force  topographical  and  (b)  Ml 
near-field  optical  images  of  a  PZT  sample  surface.  The  size  of  the  images  is 
1.2  /xm  hy  3  fim.  Height  variation  within  the  topographical  image  is  5H0 
mu.  The  nent-licid  image  has  been  measured  using  /’'polarized  excitation 
light.  Crystallite  boundaries  are  indicated  by  arrows. 

We  have  used  a  commercial  radially  poled  PZT  lube  as  a 
sample.  Measurements  have  been  performed  at  a  point  A  (see 
I  ip,  I)  where  the  electric  field  of  /'-polarized  fundamental 
light  Ims  a  component  along  the  macroscopic  poling  direc¬ 


tion  of  the  piezotube.  Electric  field  of  S-polarized  light  is 
perpendicular  to  the  poling  direction.  Simultaneously  mea¬ 
sured  topographical  and  near-field  SH  images  of  the  surface 
of  the  PZT  sample  obtained  using  P-polarized  fundamental 
light  are  shown  in  Fig.  2.  The  polycrystalline  structure  of  the 
ceramic  exhibits  itself  in  these  images.  Different  crystallites 
show  different  brightness  of  SH  emission  that  corresponds  to 
variations  in  local  poling.  These  variations  should  contribute 
to  the  effects  of  hysteresis  and  creep  of  piezoceramic.  Some 
crystallite  boundaries  also  show  up  in  the  SH  image  (as  is 
indicated  by  arrows).  They  are  believed  to  contribute  to  the 
nonlinear  behavior  of  PZT  ceramic.6 

Symmetry  considerations  show  that  the  second  harmonic 
polarization  tensor  d)fk  that  relates  SH  and  the  fundamental 
excitation  light  (  P{2m=  €t)d^jkEfEf)  has  the  same  form  as 
the  piezoelectric  tensor  dijk ,  which  gives  the  electric  polar¬ 
ization  Pj  at  a  given  stress  (rJk :  Pt-  diJkaJk  J  Thus,  it  should 
be  possible  to  study  local  variations  of  the  piezoelectric  ten¬ 
sor  of  the  PZT  ceramic  dijk  by  looking  at  the  local  SH  re¬ 
sponse  of  the  material  lor  different  polarization  states  of  in¬ 
coming  fundamental  light.  Indeed,  switching  between  P-  and 
.V-poInriml  excitation  light  leads  to  substantial  changes  in 
the  near-lick!  SH  images  of  PZT  ceramic,  as  it  is  evident 
from  l;ig.  3.  Differences  in  polarization  directions  of  indi¬ 
vidual  crystallites  exhibit  themselves  as  variations  of  SH 
brightness  for  different  polarizations  of  excitation  light.  A 
complete  recovery  of  the  local  poling  direction  should  be 
possible  if  near-field  SH  measurements  are  made  at  different 
polarization  stales  and  different  angles  of  incidence  of  the 


MG,  3,  Simultaneously  measured  (a)  shear-force  topographical  and  SH  near-field  optical  images  of  a  PZT  sample  surface  measured  using  (b)  P~  and  (c) 
.V-polari/cd  exctlahon  light.  Size  of  the  images  is  1.7  ^mX3  gm.  Height  variation  within  the  topographical  image  is  350  nm.  Image  (d)  represents  a  schematic 
of  tlifh'frnl  vulblc  in  the  images  («)~(e)  with  the  in-plane  components  of  the  local  poling  vector  shown  by  the  arrows. 
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fundamental  light  together  with  the  control  of  SH  light  po¬ 
larization. 

Individual  crystallites  of  PZT  have  a  perovskite  crystal 
structure  and  belong  to  the  tetragonal  symmetry  class  4mm 
—  C4u.  BaTi03  is  a  very  similar  piezoelectric  material  with 
the  measured  nonlinear  susceptibilities8  of  d2^.  —  d2^ 
=  =4"  =  "  17.7X  10” 12  m/V,  d:^~d2z^=~\S.S 

X  10“ 12  m/V,  and  d2zzz  =  —  7.1  X  10“  12  m/V,  where  z  is  the 
poling  direction.  A  three  times  smaller  value  of  d2zz  with 
respect  to  the  other  components  of  the  SHG  tensor  results 
from  the  much  smaller  value  of  the  permittivity  ezz  of  a 
piezoelectric  single  domain  crystal  in  the  poling  direction. 
Thus,  local  SH  intensity  should  have  a  pronounced  minimum 
(equal  to  approximately  one  tenth  of  its  average  value)  when 
the  electric  field  of  fundamental  excitation  light  is  polarized 
along  the  local  z  direction.  This  is  one  of  the  recipes  for  local 
poling  direction  recovery,  although  it  may  not  be  the  most 
convenient  way  since  it  requires  angle  of  incidence  adjust¬ 
ment  of  the  fundamental  excitation  light. 

It  is  much  easier  to  analyze  data  in  Fig.  3  if  an  assump¬ 
tion  is  made  that  the  local  poling  direction  deviates  slightly 
from  the  macroscopic  one  (this  assumption  should  be  true  in 
the  case  of  a  poled  PZT  sample).  If  the  local  poling  direction 
is  rotated  by  a  small  angle  6  around  the  X  axis  (see  Fig.  1) 
the  electric  field  of  the  S-polarized  fundamental  excitation 
light  remains  perpendicular  to  the  poling  direction.  No 
changes  in  local  SHG  should  be  expected.  On  the  other  hand, 
if  the  local  poling  direction  is  rotated  by  a  small  angle  <f> 
around  the  Y  axis,  the  components  of  second  harmonic  po¬ 
larization  change  as  follows: 

P[2)=e0d2z“E2cos2<i> 

and 

P[2)  =  e0 d]“.E2  sin2  (1) 


where  E  is  the  fundamental  electric  field.  Thus,  the  SH  im¬ 
age  in  Fig.  3(c)  reflects  variations  in  cos4<t>.  Analysis  of  Fig. 
3(b)  is  more  complicated  since  variations  in  0  and  are 
mixed  in  the  case  of  P-polarized  excitation.  However,  ex¬ 
pressions  similar  to  Eq.  (1)  are  easy  to  derive  and  possible  to 
analyze  completely  if  the  data  for  cos43>  is  known.  This 
analysis  is  especially  easy  when  there  is  no  variation  in  local 
S-polarized  SHG  as  in  the  upper  right  part  of  Fig.  3(c).  In 
this  case  variations  of  local  SHG  with  P-polarized  excitation 
correspond  to  variations  in  ©.  Figure  3(d)  represents  a  sche¬ 
matic  view  of  different  grains  visible  in  the  images  in  Figs. 
3(a)-3(c)  with  the  in-plane  components  of  the  local  poling 
vector  shown  by  the  arrows.  The  directions  of  local  poling 
vectors  were  derived  using  the  procedure  described  above. 

In  conclusion,  we  have  developed  a  microscopy  tech¬ 
nique  to  study  local  nonlinear  properties  of  PZT  piezocer¬ 
amic.  The  technique  allows  us  to  determine  the  local  poling 
direction  of  individual  submicrometer  size  crystalline  grains 
of  ceramic  by  near-field  second  harmonic  imaging  for  differ¬ 
ent  angles  of  incidence  and  polarization  states  of  fundamen¬ 
tal  excitation  light. 
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We  report  observation  of  strongly  enhanced  surface  second  harmonic  generation  (SHG)  in  BaTi03 
due  to  excitation  of  a  photorefractive  surface  electromagnetic  wave.  Surface  SH  intensity  may  reach 
10  2  of  the  incident  fundamental  light  intensity.  Angular,  crystal  orientation,  and  polarization  depen¬ 
dencies  of  this  SHG  are  presented.  Possible  applications  of  this  effect  in  nonlinear  surface  spectroscopy 
are  discussed. 
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Linear  surface  electromagnetic  waves  (SEW)  such  as 
surface  plasmons  or  surface  polaritons  [1]  play  a  veiy  im¬ 
portant  role  in  such  surface  optical  phenomena  as  surface 
enhanced  Raman  scattering,  surface  second  harmonic  gen¬ 
eration  (SHG),  etc.  They  are  extremely  useful  in  appli¬ 
cations  such  as  chemical  and  biological  sensing  since  the 
electromagnetic  field  of  a  SEW  is  strongly  enhanced  near 
the  interface.  A  linear  SEW  may  be  excited  at  the  inter¬ 
face  between  media  with  opposite  signs  for  their  dielectric 
constants  e,  such  as  at  a  metal-vacuum  interface.  Another 
example  is  the  interface  between  a  vacuum  and  a  dielec¬ 
tric  that  has  a  sharp  absorption  line.  Such  a  dielectric  has 
e(a>)  <  0  for  frequencies  just  above  the  absorption  line. 
In  both  cases  the  SEW  free  propagation  length  in  the  vis¬ 
ible  range  does  not  exceed  a  few  micrometers  because  of 
high  losses  [1],  This  limits  the  potential  advantages  of  us¬ 
ing  SEW  in  surface  enhanced  nonlinear  optical  studies  and 
sensor  applications. 

Recently,  a  new  kind  of  nonlinear  SEW  called  a  pho¬ 
torefractive  surface  wave  has  been  predicted  [2]  and  ob¬ 
served  experimentally  [3]  in  BaTi03.  This  phenomenon  is 
closely  related  to  self-trapped  optical  beams  (also  known  as 
spatial  solitons)  and  self-bending  beams  observed  in  pho¬ 
torefractive  crystals  [4],  A  photorefractive  SEW  occurs 
when  a  beam  self-bent  towards  the  positive  direction  of 
the  optical  axis  (the  poling  direction)  undergoes  a  cycle 
of  deflections  towards  the  face  of  the  crystal  and  total  in¬ 
ternal  reflections.  The  resulting  nonlinear  SEW  is  local¬ 
ized  near  the  crystal-air  interface  with  a  penetration  depth 
as  small  as  10  ptm  into  the  photorefractive  crystal  [2], 
This  leads  to  a  strong  enhancement  of  the  optical  field 
near  the  interface  that  is  common  for  all  SEWs.  On  the 
other  hand,  since  BaTi03  is  transparent  in  the  visible  range, 
the  free  propagation  length  of  the  photorefractive  SEW 
along  the  surface  is  limited  only  by  the  size  of  the  crys¬ 
tal.  As  a  result,  a  very  strong  enhancement  of  all  nonlin¬ 
ear  surface  optical  phenomena  (such  as  surface  adsorbed 
molecular  luminescence,  Raman  scattering,  surface  SHG, 
etc.)  may  be  expected  due  to  photorefractive  SEW  excita¬ 
tion.  This  effect  may  also  be  used  in  combination  with 
further  field  enhancements  produced  by  surface  topo- 
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graphical  defects  or  by  the  probe  tip  of  a  scanning  probe 
microscope  [5], 

In  this  Letter  we  report  the  first  observation  of  strongly 
enhanced  surface  SHG  due  to  the  photorefractive  SEW 
excitation  in  BaTi03.  Surface  SHG  in  BaTi03  is  a  very 
suitable  phenomenon  for  demonstrating  the  potential 
of  photorefractive  SEWs  in  nonlinear  surface  optics. 
Phase-matched  optical  SHG  (collinear  or  noncollinear) 
is  forbidden  in  the  bulk  of  BaTi03  in  the  visible  range 
because  of  strong  dispersion:  The  refractive  indices 
for  ordinary  and  extraordinary  waves  are  2.67  and  2.57, 
respectively,  at  400  nm,  and  2.36  and  2.32,  respectively, 
at  800  nm  light  wavelength  [6],  The  momentum  conser¬ 
vation  law  cannot  be  satisfied  in  the  volume  SHG  process. 
Thus,  in  contrast  to  experiments  done  with  KTP  crystals 
[7],  observation  of  the  fundamental  and  phase-matched 
SHG  fields  mutually  trapped  in  the  volume  spatial  solitary 
wave  is  impossible  in  BaTi03.  On  the  other  hand,  phase¬ 
matching  conditions  are  modified  for  the  surface  SHG. 
Only  the  momentum  component  parallel  to  the  interface 
must  be  conserved.  This  leads  to  an  extremely  strong 
phase-matched  surface  SHG  when  a  photorefractive  SEW 
is  excited. 

The  surface  nature  of  this  SHG,  in  combination  with 
the  strong  photorefractivity  of  BaTi03,  leads  to  quite  pe¬ 
culiar  angular  and  orientational  behavior  of  the  SH  light 
intensity.  We  believe  that  our  observations  performed  in 
a  well  controlled  geometry  may  clarify  a  lot  of  questions 
concerning  still  unclear  “anomalous”  SHG  in  BaTi03  re¬ 
ported  recently  by  a  number  of  groups  [8,9].  This  is 
especially  important  for  the  rapidly  developing  field  of 
ferroelectric  oxide  film  growth  and  characterization. 
Much  recent  effort  in  this  field  is  caused  by  the  ap¬ 
plications  of  these  films  in  nonvolatile  ferroelectric 
random-access  memory  and  dynamic  random-access 
memory  devices  [10],  SHG  has  been  used  to  determine 
the  crystallographic  orientation  and  the  degree  of  poling 
of  these  films  (in  particular,  BaTi03  films  in  [9]). 

Our  experimental  geometry  is  shown  in  Fig.  1.  The 
iron  doped  crystal  of  BaTi03  from  Sanders  Inc.  used  in 
the  experiments  is  a  8  mm  X  8  mm  X  8  mm  cube  which 
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was  cut  with  two  opposing  faces  perpendicular  to  the  c 
axis  and  poled  along  the  c  axis.  The  crystal  was  mounted 
on  a  three-axis  translational  and  rotational  stage.  Linearly 
polarized  light  from  a  Ti:sapphire  laser  operating  at  a 
wavelength  of  810  nm  (repetition  rate  76  MHz,  100  fs 
pulse  duration,  30  mW  average  power)  was  focused  onto 
the  top  edge  of  the  crystal.  The  direction  of  the  beam  was 
parallel  to  the  top  face  of  the  crystal. 

As  the  vertical  position  of  the  crystal  was  scanned,  very 
bright  second  harmonic  light  emission  coining  from  the 
top  face  of  the  crystal  was  observed.  The  brightness  of 
SH  emission  was  very  sensitive  to  the  vertical  position  of 
the  crystal  as  is  evident  from  Fig.  2(a),  which  was  obtained 
using  a  focusing  lens  with  a  focal  length  of  60  mm.  No 
SH  emission  had  been  detected  when  the  fundamental  light 
passed  through  the  crystal  in  any  direction  far  from  the  top 
face.  This  is  consistent  with  the  fact  that  phase-matched 
SHG  is  prohibited  in  the  bulk  of  BaTiC>3  crystal.  Also, 
no  comparable  SHG  has  been  detected  from  the  other  five 
faces  of  the  same  crystal.  The  full  width  at  half  maximum 
of  the  SH  peak  in  Fig.  2(a)  is  equal  to  80  /im. 

The  SH  emission  appeared  to  be  localized  in  the  xy 
plane  ,  of  the  top  face  of  the  crystal  coming  out  of  the 
top  face  within  a  wide  (almost  180°)  angle.  This  peculiar 
spatial  distribution  of  SHG  is  illustrated  in  Fig.  3  which 
shows  a  pattern  of  SH  illumination  on  the  screen  placed 
just  behind  the  crystal  in  the  x z  plane.  The  SH  light 
is  vertically  polarized.  Its  intensity  is  proportional  to 
the  square  of  the  fundamental  light  intensity  and  depends 
strongly  on  the  polarization  state  of  the  linearly  polarized 
fundamental  light  [Fig.  2(b)].  The  SH  intensity  is  approxi¬ 
mately  6  times  stronger  in  the  case  of  fundamental  light 
linearly  polarized  perpendicular  to  the  optical  axis  of  the 
BaTiC>3  crystal. 

All  features  of  the  observed  SH  emission  are  consistent 
with  the  proposed  surface  SHG  enhancement  due  to  the 
photorefractive  surface  wave  excitation  in  BaTi03.  The 
crystal  cut  allows  photorefractive  SEW  propagation  only 
on  the  top  face  of  the  crystal  (only  near  this  face  can 
photorefractive  self-bending  and  total  internal  reflection 
compete  with  each  other).  This  is  consistent  with  the 
observation  of  SHG  only  from  the  top  face  of  the  crystal. 
At  the  same  time,  the  full  width  at  half  maximum  of 
the  vertical  position  dependence  of  SHG  in  Fig.  2(a) 
corresponds  to  the  observed  range  of  photorefractive  SEW 
launching  in  BaTiC>3  [3]. 


BaTiOj 

FIG.  1.  Schematic  view  of  our  experimental  geometry. 


An  explanation  of  the  angular  distribution  of  the  ob¬ 
served  SHG  needs  a  more  detailed  analysis.  Nonlinear 
SEW  solutions  of  the  Maxwell  equations  in  a  photorefrac¬ 
tive  medium  may  be  found  as  follows  [2].  It  is  assumed 
that  a  SEW  written  as  E(z,  y)  =  E(z)  exp  (—iky)  induces  a 
dielectric  constant  change  of  the  form  e(z)  =  e  +  5e(z) 
due  to  the  photorefractive  effect  [here  E{z)  is  supposed 
to  be  real,  e  is  the  real  dielectric  constant  of  the  medium 
in  the  absence  of  the  wave,  the  photorefractive  medium  is 
supposed  to  be  optically  isotropic,  and  the  field  distribu¬ 
tion  is  assumed  to  be  independent  of  x.  These  simplifying 
assumptions  may  not  be  true  in  the  most  general  case,  but 
they  allow  us  to  illustrate  the  basic  physics  of  the  phenome¬ 
non].  Substituting  this  into  Maxwell’s  equations  results  in 
the  following  equation  for  E(z): 

[d2/dz2  +  (*o  ~  k2)  +  klde{z)/e']E{z)  =  0,  (1) 

where  ko  —  (o(e€o^o)^2  is  the  wave  number  of  the 
light  in  a  linear  medium  with  the  same  unperturbed 
dielectric  constant  e.  Assuming  the  diffusion  mechanism 
for  nonlinearity  [11],  the  distribution  Se(z)  may  be  related 


Polarization  direction  of  the  fundamental  light  (degrees) 


FIG.  2.  (a)  SH  intensity  as  a  function  of  the  vertical  position  z 

of  the  crystal.  Large  positive  z  corresponds  to  the  fundamental 
light  passing  through  the  crystal.  Approximate  position  of 
the  focal  spot  of  the  lens  is  shown  by  the  arrow,  (b)  SH 
intensity  as  a  function  of  the  polarization  direction  of  the 
linearly  polarized  fundamental  light,  90°  corresponds  to  the 
polarization  direction  parallel  to  the  optical  axis  of  the  BaTiCb 
crystal  ( z  direction). 
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FIG.  3.  The  pattern  of  SH  illumination  produced  on  a  screen 
placed  just  behind  the  crystal.  A  bandpass  filter  is  used  to  cut 
off  the  illumination  by  fundamental  light. 


to  E(z)  through  the  space-charge  electric  field  Esc: 

Esc  =  -(kBT/e)[dI(z)/dz]/I(z),  (2) 

where  I(z)  is  the  intensity  of  light.  The  space-charge 
electric  field  induces  refractive  index  changes  via  the 
electro-optic  effect  [11]: 

8e(z)  =  2nAr(kBT/e)  [dE{z)  /  dz]/ E{z) ,  (3) 

where  r  is  the  linear  electro-optic  coefficient.  Thus,  we 
obtain  the  wave  equation  in  the  form 

[{d2/dz2)/kl  +  (2yd/dz)/k0  - 

2(k  -  k0)/k0]E(z)  =  0,  (4) 

where  y  =  kon2r(kBT/e).  This  equation  has  a  solution 
exponentially  decaying  into  the  photorefractive  medium; 

E(z)  =  exp(— y£0z)cos{[2(£  -  k0)k0]x/2z  +  <t>}.  (5) 

In  the  case  of  BaTi03  the  photorefractive  SEW  field 
penetration  depth  is  dz  =  (y£o)_1  ~  10  /im  [2],  Since 
the  resulting  equation  (4)  is  linear,  this  penetration  depth 
does  not  depend  on  the  field  intensity. 

Real  laser  beams  have  finite  width  in  x,  but  if  the  width 
of  the  beam  is  much  larger  than  dz  the  simplified  theory 
described  above  should  be  applicable.  Thus,  a  real  life 
photorefractive  SEW  must  have  a  localization  length  in 
the  x  direction  (dx)  much  larger  than  dz.  Indeed,  this  is 
evident  from  the  profile  of  the  photorefractive  SEW  in 
BaTi03  measured  in  [3],  where  dx  ~  500  /im.  Unlike 
dz,  dx  may  depend  on  the  intensity  of  SEW. 

A  very  fruitful  approach  to  all  phenomena  related  to 
spatial  solitons  is  the  representation  of  solitons  as  linear 
waves  propagating  in  the  self-induced  optical  waveguides 
[12,13].  Let  us  follow  this  way  of  thinking  and  consider 
fundamental  light  propagating  in  a  surface  waveguide 
with  a  rectangular  dz  X  dx  profile  such  as  dx  »  dz  ~ 
10  /im  [Fig.  4(a)]  (this  is  an  approximation  since  such 
a  self-induced  waveguide  does  not  have  sharp  physical 
boundaries,  but  this  approximation  works  fine  for  spatial 


5  mm 


FIG.  4.  (a)  Phase-matched  SHG  in  the  self-induced  surface 
waveguide,  (b)  A  picture  of  the  top  face  of  the  crystal  taken 
through  the  bandpass  filter  which  cuts  off  the  fundamental 
illumination  light.  Simultaneous  weak  illumination  with  the 
white  light  from  the  flashlight  allows  the  edges  of  the  BaTi03 
crystal  to  show.  The  focused  fundamental  beam  was  coming 
from  the  top  right  comer  of  the  image.  The  surface  beam  of 
SH  light  trapped  by  the  propagating  photorefractive  SEW  is 
clearly  visible. 

solitons).  The  phase  matching  conditions  for  SHG  in 
the  surface  waveguide  are  further  modified.  Only  the 
momentum  component  parallel  to  the  waveguide  direction 
must  be  conserved.  At  the  same  time,  such  a  waveguide  is 
highly  multimode  with  many  optical  modes  corresponding 
to  geometrical  optics  rays  propagating  in  a  zigzag  manner 
parallel  to  the  surface.  Thus,  there  will  always  be  a  SH 
waveguide  mode  phase  matched  with  the  fundamental 
light.  Upon  resonant  excitation,  such  a  mode  will  be 
coupled  to  many  other  modes  existing  in  the  waveguide  at 
the  same  frequency.  As  a  result,  when  the  SH  light  reaches 
the  edge  of  the  crystal,  it  leaves  the  surface  waveguide  in 
a  ray  fan  parallel  to  the  face  of  the  crystal.  This  results  in 
the  pattern  of  SH  illumination  detected  in  Fig.  3. 

We  have  performed  further  experiments  to  confirm 
the  physical  picture  of  the  phenomenon  described  above. 
Using  a  3.6-mm  focal  length  40X  microscope  objective  as 
a  focusing  lens  in  Fig.  1,  we  have  observed  much  brighter 
surface  SH  emission  with  an  average  power  on  the  order 
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of  some  tens  of  microwatts.  No  damage  to  the  surface 
of  the  BaTi03  crystal  was  apparent  after  this  experiment. 
Thus,  as  much  as  10-2  of  the  fundamental  light  power 
had  been  converted  into  SH  emission.  This  is  a  much 
higher  conversion  efficiency  than  is  usually  observed  in 
surface  SHG  experiments:  simple  calculations  show  that 
as  much  as  107  second  harmonic  photons  per  1  nJ  laser 
pulse  have  been  generated.  This  is  3  orders  of  magnitude 
higher  than  the  conversion  efficiency  observed  in  the  case 
of  surface  SHG  enhanced  by  a  surface  plasma  wave  [14]. 
It  is  comparable  (just  an  order  of  magnitude  smaller) 
with  the  best  SH  conversion  efficiencies  obtained  in  any 
other  experimental  geometry.  At  such  a  high  power  self- 
pumped  phase  conjugation  builds  up  with  the  laser  power 
dependent  time  constant  on  the  order  of  a  few  hundreds  of 
milliseconds,  causing  instability  of  the  Ti:sapphire  laser. 
So,  optical  isolation  is  necessary  for  SHG  to  be  stable 
(which  itself  exhibits  a  similar  characteristic  buildup  time 
following  sudden  changes  in  optical  alignment). 

A  picture  of  the  top  face  of  the  crystal  taken  through  the 
bandpass  filter  (which  cuts  off  the  fundamental  excitation 
light)  is  shown  in  Fig.  4(b).  Simultaneous  weak  illumi¬ 
nation  with  the  white  light  from  a  flashlight  allows  us  to 
show  the  edges  of  the  BaTi03  crystal.  The  focused  funda¬ 
mental  beam  was  coming  from  the  top  right  comer  of  the 
image.  The  surface  beam  of  SH  light  trapped  by  the  propa¬ 
gating  photorefractive  SEW  is  clearly  visible.  The  beam 
starts  from  the  focal  point  of  the  microscope  objective.  It 
is  possible  to  see  the  SHG  because  of  surface  scattering. 
The  width  of  the  self-induced  surface  waveguide  is  esti¬ 
mated  to  be  on  the  order  of  400  fim. 

The  polarization  properties  of  SH  emission  mostly 

reflect  the  properties  of  the  second  harmonic  polarization 

(2) 

tensor  dijk  of  BaTi03,  which  relates  SH  and  the  funda¬ 
mental  excitation  light  (pf*  =  e0 df^E^E^).  For  a 
surface  SH  process  in  the  geometry  of  our  experiment 

only  two  nonzero  components  of  dijk  are  available:  dlxx  = 
-18.8  X  10"12  m/V  and  dg\  =  -7.1  X  10”12  m/V 
[15].  Thus,  the  ratio  of  SH  intensity  for  excitation  with 
fundamental  light  linearly  polarized  in  the  x  (horizontal) 
and  z  (vertical)  directions  should  be  (d^/d^)2  =  7.0 
which  is  very  close  to  the  ratio  observed  in  the  experiment 
[Fig.  2(b)]. 

At  the  same  time,  the  photorefractive  coupling  constant 
r  is  smaller  in  the  case  of  horizontally  (ordinary)  polarized 
light.  This  means  that  the  photorefractive  SEW  field  pene¬ 
tration  depth  dz  is  bigger  in  this  case  and  the  fundamental 
excitation  light  spends  more  time  away  from  the  surface. 
According  to  the  simple  self-induced  waveguide  model  de¬ 
scribed  above,  we  should  expect  stronger  attenuation  of 
SHG  under  these  circumstances.  In  order  to  account  for 
this  discrepancy  a  complimentary  point  of  view  on  the  na¬ 
ture  of  SHG  enhancement  due  to  the  photorefractive  SEW 
excitation  may  be  suggested.  In  the  picture  of  a  SEW  as  a 
beam  undergoing  a  cycle  of  photorefractive  deflections  to¬ 
wards  the  face  of  the  crystal  and  total  internal  reflections, 


a  periodic  self-induced  modulation  of  the  refractive  index 
near  the  surface  of  BaTi03  crystal  may  be  expected.  Such 
a  situation  would  closely  resemble  the  geometry  of  SHG 
experiments  with  periodically  poled  nonlinear  crystals 
such  as  lithium  niobate,  which  show  substantial  enhance¬ 
ment  of  SHG  efficiency  (periodic  poling  creates  periodic 
modulation  of  refractive  index  near  the  interface). 

In  conclusion,  we  have  observed  strongly  enhanced  sur¬ 
face  SHG  in  BaTi03  due  to  the  excitation  of  a  pho¬ 
torefractive  surface  electromagnetic  wave.  The  surface 
SH  intensity  may  reach  10“2  of  the  incident  fundamen¬ 
tal  light  intensity.  A  physical  picture  of  this  SHG  has 
been  introduced  that  assumes  phase-matched  SHG  in  the 
self-induced  surface  waveguide.  Peculiar  angular,  crys¬ 
tal  orientation,  and  polarization  dependencies  of  this  SHG 
are  presented.  They  may  account  for  anomalous  SHG 
in  BaTi03  recently  reported  in  the  literature.  The  ob¬ 
served  phenomenon  may  have  many  potential  applications 
in  nonlinear  surface  spectroscopy.  Also,  it  may  be  con¬ 
sidered  as  an  example  of  guiding  of  light  with  light  in  a 
self-induced  surface  waveguide,  which  may  find  applica¬ 
tions  in  such  emerging  soliton  related  techniques  as  writ¬ 
ing  virtual  photonic  circuits  [16],  etc. 
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Recently  a  new  kind  of  nonlinear  surface  elec¬ 
tromagnetic  wave  (SEW)  called  a  photorefractive 
surface  wave  was  predicted1  and  observed  experi¬ 
mentally2  in  BaTiCV  This  phenomenon  is  closely 
related  to  self-trapped  optical  beams  (also  known  as 
spatial  solitons)  and  self-bending  beams  observed  in 
photorefractive  crystals.3  A  photorefractive  SEW 
occurs  when  a  beam  that  is  self-bent  toward  the  posi¬ 
tive  direction  of  the  optical  axis  (the  poling  direction) 
undergoes  a  cycle  of  deflections  toward  the  face  of 
the  crystal  and  subsequent  total  internal  reflections. 
The  resultant  nonlinear  SEW  is  localized  near  the 
crystal -air  interface  with  a  penetration  depth  as 
small  as  10  /xm  into  the  photorefractive  crystal.1  This 
localization  leads  to  strong  enhancement  of  the  op¬ 
tical  field  near  the  interface,  which  should  cause  a 
significant  enhancement  of  all  nonlinear  surface  op¬ 
tical  phenomena  [such  as  surface-adsorbed  molecular 
luminescence,  Raman  scattering,  and  surface  second- 
harmonic  generation  (SHG)].  Giant  enhancement  of 
surface  SHG  in  BaTiOs  as  the  result  of  photorefractive 
SEW  excitation  was  observed  recently.4 

All  SEW’s  are  sensitive  to  the  optical  properties  of 
the  guiding  interface,  which  are  different  from  the 
properties  of  the  bulk  medium.  The  concept  of  the 
transition  layer  is  usually  used  to  describe  the  behavior 
of  the  dielectric  constant  (and  other  optical  parameters) 
near  the  interface  in  a  macroscopic  electrodynamics  ap¬ 
proach.  Detailed  measurements  of  the  SEW  field  dis¬ 
tribution  and  comparison  of  these  data  with  theoretical 
models  may  provide  valuable  information  on  the  prop¬ 
erties  of  transition  layers. 

In  this  Letter  we  report  near-field  optical  measure¬ 
ments  performed  on  photorefractive  surface  waves  ex¬ 
cited  on  the  surface  of  a  BaTiOs  crystal.  The  field 
distribution  of  the  photorefractive  surface  wave  near 
the  crystal -air  interface  was  measured  and  com¬ 
pared  with  theory.  Experimental  data  indicate  that 
a  micrometer-wide  transition  layer  with  dielectric  and 
photorefractive  properties  that  are  different  from  the 
properties  of  bulk  BaTiC>3  exists  near  the  interface. 

Our  experimental  geometry  is  shown  in  Fig.  1.  The 
iron-doped  crystal  of  BaTi03  (from  Sanders,  Inc.)  used 
in  the  experiments  is  a  8  mm  X  8  mm  X  8  mm  cube 

0146-9592/99/191367-03$15.00/0 


that  was  cut  with  two  opposing  faces  perpendicular  to 
the  c  axis  and  poled  along  the  c  axis.  This  crystal  cut 
permits  photorefractive  SEW  propagation  only  on  the 
top  face  of  the  crystal  (only  near  this  face  can  photore¬ 
fractive  self-bending  and  total  internal  reflection  com¬ 
pete  with  each  other).  The  crystal  was  mounted  upon 
the  three-axis  translation  sample  stage  of  a  near-field 
optical  microscope  described  elsewhere.6  Horizontally 
(ordinarily)  polarized  light  from  an  argon-ion  laser  op¬ 
erating  at  488  nm  was  focused  with  a  40  X  microscope 
objective  onto  the  top  edge  of  the  crystal  within  a  few 
micrometers  of  the  top  face.  The  direction  of  the  beam 
was  parallel  to  the  top  face  of  the  crystal.  This  ge¬ 
ometry  corresponds  to  the  geometry  of  photorefractive 
SEW  excitation  reported  in  Ref.  2.  Near-field  mea¬ 
surements  were  performed  on  the  top  face  and  on  the 
side  (exit)  face  of  the  BaTiOa  crystal. 

Figure  2  shows  the  data  for  field  distribution  mea¬ 
surements  performed  near  the  edge  of  the  exit  face  of 
the  crystal  where  the  photorefractive  SEW  was  leav¬ 
ing  the  surface  of  the  crystal.  The  position  of  the  edge 
is  determined  from  the  topographical  data  (shown  in 
the  same  figure  and  marked  by  the  arrow).  The  flat 
part  of  the  topographical  data  (at  distances  larger  than 
8  fi m)  corresponds  to  feedback  saturation  where  the 
tip  of  the  microscope  falls  off  the  edge  of  the  crys¬ 
tal.  The  field  distribution  calculated  according  to  the 
theory  given  in  Ref.  1  is  shown  in  the  same  figure 


NSOM  tip  position 
for  the  data  shown 
in  Fig.3 


Fig.  1.  Schematic  view  of  our  experimental  geome¬ 
try:  NSOM,  near-field  scanning  optical  microscopy. 
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Fig.  2.  Comparison  of  measured  photorefraetive  SEW 
field  distribution  near  the  edge  of  the  exit  face  of  the  crys¬ 
tal  with  theory.  The  position  of  the  edge  is  determined 
from  the  topographical  data  and  is  marked  by  an  arrow. 
The  flat  part  of  the  topographical  data  (at  distances  larger 
than  8  ^m)  corresponds  to  feedback  saturation  when  the 
tip  of  the  microscope  falls  off  the  edge  of  the  crystal. 

(labeled  theory).  Although  our  data  follow  theoretical 
predictions  quite  well,  there  is  a  substantial  differ¬ 
ence  in  the  SEW  field  behavior  near  the  interface. 
To  understand  what  parameters  of  the  model  must  be 
changed  to  account  for  this  discrepancy,  let  us  recall 
the  basic  points  of  the  theory  of  photorefraetive  SEW. 

Nonlinear  SEW  solutions  of  the  Maxwell  equations 
in  a  photorefraetive  medium  can  be  found  as  fol¬ 
lows1:  It  is  assumed  that  a  SEW  written  as  E(z ,  y)  = 
E(z)exp(-iky)  induces  a  dielectric  constant  changes  of 
the  form  e(z)  =  e  +  8  e{z)  because  of  the  photorefraetive 
effect.  [Here  E(z)  is  assumed  to  be  real,  e  is  the  real 
dielectric  constant  of  the  medium  in  the  absence  of  the 
wave,  the  photorefraetive  medium  is  assumed  to  be  op¬ 
tically  isotropic,  and  the  field  distribution  is  assumed 
to  be  independent  of  x .  These  simplifying  assump¬ 
tions  may  not  be  true  in  the  most  general  case,  but 
they  allow  us  to  illustrate  the  basic  physics  of  the  phe¬ 
nomenon.]  Substituting  this  equation  into  Maxwell’s 
equations  results  in  the  following  equation  for  E(z): 

[d2/dz2  +  (kl  -  k 2)  +  kl8e{z)/e]E{z)  =0,  (1) 

where  ko  =  w(efoMo)1/2  is  the  wave  number  of  the 
light  in  a  linear  medium  with  the  same  unperturbed 
dielectric  constant  e.  Assuming  the  diffusion  mecha¬ 
nism  for  nonlinearity,6  the  distribution  of  8e(z)  can 
be  related  to  E{z)  through  the  space-charge  electric 
field  Esc: 

Esc  =  ~(kBT/e)  [dl(z)/dz]/l(z) ,  (2) 

where  I{z)  is  the  intensity  of  light.  The  space-charge 
electric  field  induces  refractive-index  changes  by 
means  of  the  electro-optic  effect6: 

8e{z)  =  2nir(kBT/e)  [dE(zydz]/E(z) ,  (3) 

where  r  is  the  linear  electro-optic  coefficient.  Thus  we 
obtain  the  wave  equation  in  the  form 

\(d2/dz2)/k%  +  (2yd/dz)/ko  ~  2(k  -  k0 )/k0]E(z)  =  0, 

(4) 


where  y  =  k0n2r(kBT/e).  Equation  (4)  has  a  solu- 
tion  that  decays  exponentially  into  the  photorefraetive 
medium: 

E{z)  =  exp(-y£02)cos{[2(£  -  k0 )ko\V2z  +  <f>}.  (5) 

In  the  case  of  BaTi03  the  photorefraetive  SEW  field 
penetration  depth  is  d2  =  (y&0)_1  ~  10  /xm.1  Because 
the  resultant  equation  (4)  is  linear,  this  penetration 
depth  does  not  depend  on  the  field  intensity.  The 
phase  (f>  is  easy  to  find  from  the  boundary  conditions. 
The  theoretical  field  distribution  found  in  this  way  is 
shown  in  Fig.  2  (labeled  theory).  Similar  results  can 
be  found  in  Ref.  1. 

The  value  of  field  penetration  depth  dz  observed  in 
the  experiment  is  close  to  the  theoretical  prediction. 
Only  the  behavior  near  the  interface  is  unaccounted  for. 
The  simplest  way  to  explain  this  difference  is  to  assume 
the  existence  of  a  transition  layer  with  dielectric  and 
photorefraetive  properties  that  are  different  from  the 
properties  of  bulk  BaTi03.  Thus  Maxwell’s  equations 
have  to  be  solved  for  a  three-layer  medium.  The 
curve  labeled  modified  theory  in  Fig.  2  represents 
the  best  fit  obtained  in  this  way.  It  corresponds 
to  a  1-yLtm-wide  transition  layer  with  a  refractive 
index  0.01  higher  than  the  refractive  index  of  the 
bulk  BaTiC>3  ( — 2.6  at  488-nm  light  wavelength)  and 
much  weaker  photorefraetive  properties  (much  smaller 
y).  The  exact  value  of  the  refractive-index  change 
is  (to  some  extent)  a  function  of  the  transition  layer 
width.  But  it  is  clear  from  the  shape  of  the  measured 
optical  signal  that  a  transition  layer  of  the  order  of 
1  fim  does  exist  near  the  interface.  Attenuation  of 
photorefraetive  properties  near  the  interface  may  be 
expected  in  a  material  with  a  diffusion  mechanism  for 
nonlinearity  (electric  charges  cannot  diffuse  into  the 
air  side  of  the  interface).  Another  likely  reason  for  the 
existence  of  the  transition  layer  is  surface  charging. 
It  is  an  interesting  theoretical  problem  to  derive  the 
properties  of  this  transition  layer  from  microscopic 
theory. 

Characteristic  near-field  images  of  the  photorefrac- 
tive  SEW  propagation  along  the  top  face  of  the  crys¬ 
tal  are  shown  in  Fig.  3.  The  optical  image  in  Fig.  3(b) 
exhibits  a  set  of  periodic  higher-  and  lower-intensity 
regions  in  the  direction  of  SEW  propagation.  This 
behavior  may  be  expected  in  a  model  of  the  photore- 
fractive  SEW  as  a  beam  undergoing  a  cycle  of  photore- 
fractive  deflections  toward  the  face  of  the  crystal  and 
subsequent  total  internal  reflections.  Thus  a  periodic 
self-induced  modulation  of  the  refractive  index  near  the 
surface  of  BaTi03  crystal  may  be  expected.  Such  a 
situation  would  closely  resemble  the  geometry  of  SHG 
experiments  with  periodically  poled  nonlinear  crystals 
such  as  lithium  niobate,  which  show  substantial  en¬ 
hancement  of  SHG  efficiency  (periodic  poling  creates 
periodic  modulation  of  the  refractive  index  near  the  in¬ 
terface).  This  analogy  may  explain  some  features  of 
strongly  enhanced  surface  SHG  in  BaTi03  owing  to  re¬ 
cently  observed  photorefraetive  SEW  excitation.4 

In  conclusion,  we  have  performed,  for  the  first 
time  to  our  knowledge,  a  near-field  optical  study  of 
photorefraetive  SEW’s  in  BaTi03.  A  micrometer- wide 
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(a) 


transition  layer  with  dielectric  and  photorefractive 
properties  that  are  different  from  the  properties  of  the 
bulk  BaTi03  was  observed.  Near-field  images  of  SEW 
propagation  along  the  interface  show  periodic  bouncing 
of  light  off  the  interface  as  a  result  of  the  competition 
between  photorefractive  self-bending  and  total  internal 
reflection. 

I.  I.  Smolyaninov’s  e-mail  address  is  smoly@ 
eng.umd.edu. 
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Fig.  3.  (a)  Topography  and  (b)  near-field  optical  image  of 
photorefractive  SEW  propagation  along  the  top  face  of  the 
crystal. 


